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Abstract 


A description is given of the flow of two superimposed layers of fluid over a barrier. This 
represents a partial experimental investigation of a problem considered theoretically in Part I. 

In general three regimes of motion are possible: If the velocities of the fluids are sufficiently 
small the interface is little disturbed except for a slight depression over the barrier. If the velocities 
are sufficiently high the interface swells symmetrically over the obstacle. At intermediate speeds 
a hydraulic jump occurs in the lee of the barrier and the lower layer increases in depth upstream. 

Two occurrences do not fit into the above description: If the obstacle is small compared to 
the depth of the lower layer, weak lee waves appear at low speeds, increasing in amplitude as 
the approach velocity of the fluid is increased. This seems to be the only case in which pertur- 
bation theory provides an adequate prediction of the flow. The second anomalous occurrence 
is the appearance of a “jump down’ or hydraulic “‘drop’”’ in the lee when the speed of the 
fluid is moderately high, the obstacle large, and the upper fluid relatively thin. 

The description of the experiments is supplemented by a theoretical discussion, employing 
the assumption of a hydrostatic pressure distribution. In general this theory provides a satis- 
factory explanation of the observed behavior. The paper concludes with a discussion of mete- 


orological implications. 


I. Introduction 


Part I of this series of papers (LONG, 1953 a) 
contained a rather general theoretical analysis 
of the two-dimensional flow of a stratified 
liquid in a gravity field. Assuming only that 
the liquid was frictionless and the flow steady, 
the equations of motion and continuity were 
integrated once to yield a second-order, partial 
differential equation. This equation, valid for 
any arbitrary basic density stratification and 
velocity distribution, was examined from 
several viewpoints: In the first place it was 
assumed that a solution existed for given 
finite configurations of the boundary. A 
sufficient condition was then obtained for this 
solution to be unique. By analogy with 


1 This research was supported by the U.S. Weather 
Bureau. 
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hydraulics a regime of flow satisfying the 
uniqueness criterion was called supercritical. 
If only infinitesimal perturbations are con- 
sidered, uniqueness implies that the flow is in 
such a state that no free wave can maintain its 
position with respect to the fixed bottom, but 
must instead be swept downstream. This means 
that the wave energy generated by an obstacle 
in the stream cannot be propagated upstream 
since the maximum group velocity in such 
systems is the phase velocity of long waves. 
In the flow of water in a channel the occur- 
rence of a supercritical regime at some section 
is a necessary condition for a hydraulic jump. 
In the problems of water flow with a free 
surface, and high speed flow of a compressible 
fluid, this shock-wave phenomenon is by far 
the most important characteristic of the motion. 
In the case of a stratified fluid such as the 
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atmosphere there seems little doubt not only 
that shocks exist but that they are of very 
fundamental importance to meteorology. TEP- 
PER (1952), for example, assumes that squall 
lines are really surges (moving hydraulic 
jumps) on an inversion surface in the tropo- 
sphere. He then applies a simplified mathe- 
matical theory of such shock waves to actual 
surface pressure observations and predicts 
motion and development. Despite some 
objections (Furks, 1951), it is likely that 
Tepper has given a correct physical inter- 
pretation of most if not all squall lines. 

Further evidence of the existence of hydraulic 
jumps in the atmosphere is provided by 
observations of the very violent type of 
motion that occurs in the lee of the Sierra 
Nevada Range near Bishop, California (Cor- 
SON, 1952). Resemblances between this pheno- 
menon and a crude laboratory model (Long, 
1953 b) are numerous. At Froude numbers 
roughly equivalent in model and prototype 
similar violent upward motions occur in both 
cases in the lee. 

The unquestioned importance of stability in 
the atmosphere has inspired a number of 
studies of internal gravity waves set up by 
mountain barriers and other sources of pertur- 
bations (Lyra, 1943; QUENEY, 1947; SCORER, 
1949; Markus and STERN, 1953). In these 
studies, however, it is first assumed that the 
amplitudes of the perturbations are infinitely 
small in order to linearize the differential 
equations. Streamlines are then drawn in which 
the amplitudes of the waves are finite. This 
inconsistent approach raises a number of 
questions. Such a procedure yields a shock-free 
flow of a fluid past a barrier when, in fact, it 
is not known whether such a solution really 
exists. The experiments described in this paper 
indicate that hydraulic jumps occur at very 
moderate Froude numbers from a meteoro- 
logical viewpoint, provided the obstacle is not 
too small. The experiments further indicate 
that the jumps are the only phenomenon of 
any great importance unless the height of the 
obstacle is small compared to the depth of the 
lower fluid. In such cases there are small lee 
waves which increase in amplitude with the 
fluid speed. For an obstacle of appreciable size, 
however, these waves do not simply change 
in wave length and amplitude but completely 
disappear at low speeds or are replaced by 


finite jumps at the higher speeds. It is evident 
from the experiments that a finite obstacle 
cannot be considered simply as a “big” infini- 
tesimal barrier. 

In Part I mention was made of some pre- 
liminary experiments with a three-fluid system 
in which flow over an obstacle produced a 
hydraulic jump at each interface. A more basic 
approach to the experimental problem of 
investigating stratified flows, however, would 
be to set up the most simple model possible 
and attempt to understand it fully before 
constructing more complex models. Obviously 
a system of two fluids is simpler than a three- 
fluid system and the former is the subject of 
this paper. The latter will be investigated at a 
later time, along with other models of interest. 

Some mention should be made of the reason 
for the choice of multiple-layer systems instead 
of those with continuous density gradient. To 
some extent the choice is one of expediency. 
Layers of immiscible fluids can be set up with 
ease and the resulting step-wise density distri- 
bution can be maintained indefinitely. After a 
given run the fluids settle back quickly to their 
original undisturbed state and the experiments 
continue without interruption. A fluid with a 
continuous density gradient may be obtained 
by a mixture of water and salt, for example, 
(GORTLER, 1943), and some preliminary ex- 
periments with such a set-up have already been 
undertaken. Yet multiple-layer systems have 
considerable practical interest in meteorology. 
Rather strong temperature inversions are of 
frequent occurrence in the atmosphere and, 
dynamically, they resemble closely interfaces 
between two homogeneous fluids (adiabatic 
atmospheric layers!). 

Moreover, even the two-fluid case probably 
has considerable application to atmospheric 
problems. In the preliminary experimental 
investigation with three fluids it was suggested 
that the lower two fluids were analogous to 
the troposphere, and the upper fluid to the 
stratosphere. When the density difference was 
the same across each interface, small but 
significant effects (hydraulic jumps) were 
noticed at the “tropopause” for atmospheric 
values of the Froude number and obstacle sizes. 
The great stability of the stratosphere, how- 


1 The effect of compressibility will be discussed in a 
forthcoming paper. 
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ever, would indicate that a closer approxi- 


mation to the atmosphere would be obtained 
by increasing the density difference between 
the middle and upper fluid. When the latter 
density difference was made four times the 
difference across the first interface, the “tropo- 


pause” was almost undisturbed for reasonable 


Froude numbers and the fluid system behaved 
very much like a two-fluid system with an 
upper rigid or free surface. At the present stage 
of the research this should be regarded only 
as a suggestion. It is possible to arrange various 
layers of immiscible fluids to obtain density 
distributions more closely resembling those in 
the atmosphere. Such experiments will help to 
resolve some of the questions that arise in 
applications of simple stratified systems to 
atmospheric problems. 

Section 2 of this paper is devoted to a 
discussion of experiments with a two-fluid 
system of stratified fluids. A qualitative de- 
scription is given of the most important 
phenomena together with a number of photo- 
graphs. Section 3 contains an introduction to 
the theory of simple stratified systems. In this 
section the flow of a single fluid over a barrier 
is treated exhaustively, since a full under- 
standing of this model leads to a fuller under- 
standing of the more complicated two-fluid 
problem. 

Section 4 presents the theory of the general 
two-fluid problem of flow over a barrier and 
a comparison with the experimental observa- 
tions. The last section discusses the meteoro- 
logical implications of this work. 


2. Experimental observations 


The apparatus for the experimental work is 
a channel about 10 ft long, 5 in wide, and 20 
in high (fig. 1). The fluids are mixtures of 
water and salt and of carbon tetrachloride and 
a commercial cleaning fluid. An obstacle of 
“easy shape” is drawn by a motor drive at a 
uniform speed along the bottom.! This obstacle 
extends from one glass wall to the other. The 
density difference between layers is 0.025 gm 
cm-? and is not varied. Theoretical considera- 
tions indicate that this involves little or no loss 


1 In the early experiments an obstacle of semi-circular 
cross-section was used. The boundary-layer separation 
effects were very marked, however, and the resulting 
large energy losses to turbulence were undesirable. 
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Fig. 1. Experimental channel. 


of generality if interest is confined to models 
and prototypes in which density differences 
are small. 

It follows from dimensional considerations 
that the phenomena appearing in this experi- 
ment depend on a minimum of three non- 
dimensional numbers 


U ho b 

Vera 
Q 

where F; is the internal Froude number, hy is 
the initial height of the lower layer upstream, 
H is the total depth of the two fluids, and b 
is the maximum height of the obstacle. In the 
following description we will identify each 
experiment with the values of these three 
numbers. 

(a). Ro < 1/2. For any R, in this range and 
for moderate values of ß, a so-called absolutely 
subcritical state of flow exists provided the 
Froude number is small enough. In this state 
the fluids move over the barrier with little or 
no turbulence. The interface is level except for 
a slight symmetrical dip over the barrier. Fig. 
2 is a typical photograph of this regime of 
motion. Small values of F; and B, and large 
values of Ry are favorable for its existence. 
During the steady motion of an obstacle in the 
absolutely subcritical state, the upstream and 
downstream levels do not change sensibly and 
no deformation of the free surface can be seen. 
The only cases for which the above description 
of the absolutely subcritical flow does not hold 
are those in which ß is small compared to Ry. 
In such cases the absolutely subcritical regime, 
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Fig. 5. Undular jump phe- 
nomenon. F, = .284, Ro = 


— 3330 Br 6067 BE xcept 
for the higher Froude 
number the conditions are 
the same as in fig. 


Fig. 6. Strong hydraulic 
jump. F; = .272, Ro = .33, 


Dr — 2205: 


Fig. 7. Absolutely supercri- 
tical flow. F, = .583, Ro= 
= 6, 0 = AOS. 
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as indicated theoretically in section 4, may 
involve wave motions in the lee as in fig. 3. 
The Froude number, .220, is absolutely sub- 
critical for the values Ry = .33, P = .067. 


If the fluid levels and the obstacles are kept 
constant while the Froude number is increased, 
the absolutely subcritical regime changes, at a 
certain Froude number, to one which we may 
call critical. The interface is no longer symme- 
tric about the crest of the barrier, but continues 
to lower on the lee. Instead of returning 
smoothly to the undisturbed downstream level, 
the interface jumps more or less abruptly to 
the downstream level as in fig. 4. The interface 
has a rough appearance but assumes a relatively 
constant level downstream. With small ob- 
stacles the rise to the downstream level takes 
place over a series of downstream waves, 
resembling the undular jump in open channel 
flow. This is illustrated in fig. 5. 


If the Froude number is lowered the turbu- 
lent jump moves toward the crest, diminishing 
in intensity. The value of the Froude number 
when the jump just disappears is taken as the 
lower limit of the critical regime. 


If the Froude number is increased the jump 
moves downstream and increases in intensity. 
An example of a stronger jump is given in 
fig. 6. Unless the jump is very weak, the 
obstacle acts to block the flow of the lower 
fluid upstream and the upstream depth in- 
creases markedly. This is very evident in fig. 6. 
The depth ratio was originally .33; at the time 
of the picture it has been increased to nearly 
.50 by this blocking action. The increase is 
accomplished by a gravity wave of elevation 
which moves further and further ahead of the 
barrier. If the channel were much longer than 
our present one this wave would eventually 
be able to move so far ahead of the obstacle 
that we could neglect it and use the height of 
the lower layer just upstream to compute a new 
value of Ry. Experimental comparison with 
the theory of section 4 indicates that the channel 
is too short for this, and a true unsteady state 
exists whenever a strong jump occurs. 


If the Froude number is sufficiently high 
initially, a jump, forming in the lee of the 
obstacle, is left behind as the barrier moves on. 
In this state the blocking action does not 
appear; the fluid level upstream remains at its 
initial height. The interface swells up over the 
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barrier and is symmetric about the crest. This 
is called absolutely supercritical and is shown 
in fig. 7. 

Depending on the experimental conditions 
the absolutely supercritical regime occurs for 
Froude numbers of the order of .so and 
greater. At such speeds the free surface is 
sensibly undisturbed. If F; is now raised to 
2 or 3, the free surface becomes disturbed and 
a lee jump may even occur on it (fig. 8). This 
regime is considered to be outside the area of 
interest of this paper. 

When the initial height of the lower fluid 
is less than that of the obstacle, two types of 
flow are possible: If F; is very high the lower 
fluid moves symmetrically over the barrier in 
a swell, or for the larger obstacles, the interface 
rises over the barrier remaining high on the 
lee for some distance and then appears to 
drop down violently to a lower downstream 
level (fig. 9). This particular phenomenon is 
discussed further below. If the Froude number 
is smaller, however, the high barrier blocks 
the lower fluid completely and the lower fluid 
builds up eventually until it spills over the 
obstacle, forming a lee jump. This is brought 
out in figs 10 and 11. In fig. 10 the blocking 
action has increased the upstream depth con- 
siderably. Fig. 11, taken a short time later, 
shows the spilling over the barrier and the lee 
ump. 

(b). Ro > 1/2. When the lower fluid is 
deeper than the upper, the absolutely sub- 
critical regime still exists for sufficiently small 
Froude numbers, and a very small obstacle will 
again produce lee waves in this regime. If the 
Froude number is increased, a small jump 
appears. As in the case of Ry < 1/2, the jump 
increases in size and moves downstream with 
an increase in F;. When the obstacle is large 
or moderately large and the Froude number 
is high, instead of an ordinary lee jump, or a 
symmetrical, absolutely supercritical type of 
motion, the interface rises on the upwind side, 
reaches a maximum at the crest and then 
appears to “drop” suddenly on the lee. An 
example of the hydraulic drop is shown in 
fig. 12. An increase in F; merely causes the 
hydraulic drop to move downstream. In some 
cases, however, a drop is followed by an 
orthodox jump in the interface level, but this 
does not always happen. Although a hydraulic 
drop is not possible in a one-fluid system it is 
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Fig. 8. Jump on free surface. ss : 04 09 25 53 


F; 2:20 0 Ry = -7S D = pat ee 


== .00: 


Fig. 9. Hydraulic drop. 
The obstacle is higher than 
the upstream fluid depth. 
The interface remains high 
on the lee dropping with 
accompanying turbulence 
to the downstream level. 
F. = .577, Ro = -33, B = 
== .50. 


Fig. 10. Complete block- 
ing of lower fluid by a bar- 
rier higher than the inter- 
face: 
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a definite thecretical possibility in multiple 
fluids. In a two-fluid system it is favored by a 
relatively deep lower fluid (see section 4). 


3. Hydraulic analogy 


The theory of the flow of two fluids over an 
obstacle is developed in section 4 of this paper. 
The analogous problem of the flow of water 
over a barrier is much simpler and it seems 
advisable to develop first the theory of the 
latter system. Since this problem has long been 
of interest to hydraulic engineers, the present 
section probably does not contain any original 
material. The procedure is quite different, 
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Fig. 11.' Final spilling of 
the fluid over the barrier 
of fig. 10 and formation of 
a lee jump. 


Fig. 12. Hydraulic drop. 
‚542, Ro = .75, B= 
= .60. 


however, from the standard approach and 
seems to contain a number of advantages.! 
Fig. 13 1s a drawing of a channel in which 
water of basic depth, H, approaches a barrier 
of variable height, b(x), over the otherwise 
horizontal bottom. The water has a vertical 
depth, h(x), over the topography and the 
flow with respect to the obstacle is assumed 
steady. The fluid is considered frictionless and 
the pressure distribution hydrostatic. In the 
experimental work of this paper an obstacle is 
towed through a resting fluid. Relative motion 


* Iam indebted to Professor George S. Benton for a 
number of invaluable suggestions regarding the material 
of this section. 
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Fig. 13. Theoretical model of single-Auid flow over a 
barrier. 


of the fluid and the boundaries occurs as the 
liquid passes over the barrier. It is not likely 
that the resulting energy changes will be 
important, since turbulence is largely avoided 
by using a smooth obstacle of “easy shape”. 
The hydrostatic assumption is certainly satisfied 
upstream where the vertical motion is zero. 
In the vicinity of the barrier the paths of the 
fluid particles are necessarily curved and there 
will be vertical centripetal accelerations. 


Some mention may be made of the effect of 
dissipative forces and centrifugal forces. In the 
absence of a jump there will still be some loss 
of energy, as the fluid passes over the barrier, 
by friction with the moving obstacle. This is 
accompanied in open channel flow by a loss 
of total head (total depth of the fluid). Con- 
sequently, in a steady state the water level on 
the downstream side must be somewhat lower 
than that on the upstream side. If the channel 
is infinite in length, however, sufficiently far 
upstream and downstream the levels must be 
equal if the obstacle, starting from rest, has 
moved a finite length of time. The conclusion 
is inevitable that a steady state can never be 
achieved because of the continuing adjustment 
of upstream and downstream levels. While this 
effect will be very small in laminar motion, a 
jump, occurring upstream or in the lee of the 
obstacle, causes considerable turbulence and an 
enormous increase in energy loss. The relative 
increase of upstream head is accomplished by 
a blocking action of the barrier. A long wave 
of elevation moves upstream, raising the 
approach level sufficiently to permit the jump 
to remain in a steady state. It follows that for 
large jumps, since the approach Froude number 
is large, the elevation wave cannot move very 
far ahead of the barrier in the limited length 
of the experimental channel. A quasi-steady 
state cannot be assumed and the theoretical 
deductions in this paper are not applicable. The 
Tellus VI (1954), 2 
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blocking effect is certainly worthy of separate 
study, however, since it will be operative in 
the flow of air over mountain barriers. 

The neglect of the curvilinear nature of the 
flow by assuming a hydrostatic pressure distri- 
bution is not completely justified unless the 
characteristic radius of curvature of the bound- 
ary form is considerably greater than the depth 
of the fluid. This is not fulfilled in the ex- 
perimental work and some error results. The 
centrifugal forces act upward in the flow over 
a barrier so that the effect is to reduce gravity 
apparently. This should tend to reduce the 
critical Froude number at which lee jumps 
begin to occur in the lee. The experiments do 
not indicate an appreciable error in this regard 
provided a real shock wave occurs for a given 
height of the obstacle and fluid depth. As 
mentioned in section 4 a shock wave predicted 
by hydrostatic theory does not materialize in 
the experiments if the barrier is too small. 
Instead lee waves or so-called undular jumps 
(without accompanying turbulence) are found. 
Since a truly hydrostatic pressure distribution 
permits only infinitely long waves, hydro- 
static theory obviously breaks down in such 
cases. On the whole, however, if this limitation 
is kept in mind, the hydrostatic assumption is 
much more useful and valid than small 
perturbation theory. 

In the model of fig. 13 the Bernoulli equation 
for a steady state is 

2 
48 + ae + 2= constant, (1) 
Cg 2g 
where z is the height of the point in question 
above the horizontal bottom, and q is the speed 
at that point. Using the hydrostatic assumption 
to compute the pressure, and neglecting the 
vertical velocity1, (1) becomes 


— +h + b(x) = constant. (2) 


In hydraulics the sum of the first two terms of 
this expression is called the specific energy, e, 


of the Auid, 
e=— +h. (3) 


1 The ratio v2/u2 = o(h?/A?), where À is the distance 
between two consecutive points at which the wave 
profile meets the undisturbed level. The hydrostatic 
assumption therefore implies that the kinetic energy due 
to vertical motion is negligible. (See LAMB 1932, p. 258). 
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Fig. 14. Specific energy diagram. 


Since the equation of continuity requires a 
constant discharge, Q, at each section, 


Q = uh = constant. (4) 


Substituting for u in (3) we find that 


@z 

tere e(h). (5) 
Thus, for a given discharge, e is a function of 
the depth of the fluid, h. A plot of this relation- 
ship is given in fig. 14. For a given Q, any 
fluid depth is possible but at a certain depth, 
h., the specific energy is a minimum. The 
minimum occurs when u? = gh. Depths less 
than h, are called supercritical and the flow is 
then in such a state that a hydraulic jump is 
possible at some downstream section. The jump 
is characterized by a sudden increase in the 
depth of flow downstream, accompanied by 
turbulence. The flow downstream of the jump 
must have a depth greater than the critical 
depth, h., and so lies in the upper branch of the 
specific energy curve. Some of the energy is 
lost through turbulence as illustrated qualita- 
tively by the arrow in fig. 14. It should be 
noted that in the absence of frictional or 
turbulent losses the specific energy changes 
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only by the virtue of changes in the height of 
the bottom of the channel. As the fluid moves 
over an obstacle, for example, it loses specific 
energy on the upwind side and gains it on 
the downwind side. Quantitatively, from (2) 


and (5), 
de dedh _ -S)E--7 (6) 
dx dhdx j gh? dx dx’ 


u?\ dh db 
Us m 

gh) dx dx 
This equation shows that u? can equal gh (i. e. 
the flow can become critical) only where 
db/dx vanishes. Since the specific energy 
changes only during passage over the barrier, 
if the flow becomes critical anywhere, it must 
do so at the crest where b = b.. Moreover we 
may show that a critical condition at the crest 
necessarily implies a change of regime from 
the upstream to the downstream side of the 
barrier. Thus, differentiating (7) with respect 


COR 
zu? (dh\2 u2\ d?h d2b (8) 
gh? (=) x (: *) dx? acer 
at the top of the barrier. Since it is critical 
there the second term on the left vanishes. 
However, the right-hand side is not zero so 
that dh/dx does not vanish at the crest. Hence, 
if it is subcritical upstream and critical at the 
crest, it will be supercritical just on the lee. 
If it is supercritical upstream, and critical at 
the crest it must be subcritical on the lee side. 
With these considerations in mind we may 
discuss the various possibilities listed in Table 
1, not all of which are possible steady states. 
They may occur instantaneously, however, if 
a given obstacle is accelerated from zero very 
rapidly to a given speed. In (a) the obstacle 
speed is very slow and the flow is everywhere 
subcritical. This is an obvious steady state. If 


or 


Table 1. Regimes of motion in flow of single fluid over an obstacle 


Approach Upstream Jump Crest | Downstream Jump Downstream 
ee eae | VOR EEE 
(a) Subcritical — Subcritical — Subcritical 
(b) Subcritical — Critical yes Subcritical 
(c) Subcritical — Supercritical yes Subcritical 
(d) Supercritical yes Subcritical — Supercritical 
(e) Supercritical yes Critical | — Supercritical 
(f) Supercritical — Supercritical — Supercritical 


Tellus VI (1954), 2 


FLOW OF STRATIFIED FLUIDS 107 


Fig. 15. Subcritical flow of 
single fluid over a barrier. 
F = .044, B = .773. 


Fig. 16. Lee jump in single- 
fluid system. F = .394, 
B = .773. 


Fig. 17. Upstream jump in 
single fluid system. F = 
— 037, » = -773- Experi- 
mentally it was not pos- 
sible to move the obstacle 
fast enough to obtain a 
Froude number greater 
than 1.0. Consequently, 
this figure does not repre- 
sent a condition that quite 
falls in the region (d) of 
fig. 18. Obviously since F 
is slightly less than 1.0 the 
upstream jump could not 
maintain its position with 
respect to the obstacle be- 
cause the fluid would the be 
jumping from subcritical 
to subcritical. Thus the 
jump must be moving 
down the channel faster 
than the obstacle follow- 
ing it. 
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the speed is increased further (b), the flow 
will become critical at the crest of the barrier 
and, necessarily, supercritical somewhere on 
the lee. It is apparent from the specific energy 
diagram that the fluid cannot return to the 
subcritical state required downstream except 
by a jump on the lee. If the obstacle speed is 
increased further (c) the flow will become 
supercritical instantaneously at the crest, 
although remaining subcritical upstream and 
downstream. This is not asteady state, of course, 
and a wave of elevation! will progress upstream 
raising the upstream level sufficiently to permit 
a critical condition at the crest. A steady state 
will ultimately be established with a lee jump. 

If the speed of the barrier is so high that the 
flow is moderately supercritical upstream (d), 
the loss of specific energy in flowing over the 
barrier may reduce the system to subcritical at 
the crest. This is also unsteady and a jump 
upstream must form. We may see this by 
considering the rate of propagation of disturb- 
ance energy originating at the barrier. Since 
the flow in the vicinity of the crest is subcritical, 
the group velocity of the longer waves, tending 
to move upstream, exceeds the fluid velocity 
and energy is propagated upstream. This 
energy cannot move indefinitely upstream, 
however, because of the postulated supercritical 
flow far upstream. The resulting convergence 
of energy is then available to support a jump. 
The fluid will then jump to subcritical some- 
where upstream followed by critical flow at 
the crest, passing smoothly to the supercritical 
condition downstream. 


The obstacle speed may now be increased 
sufficiently to permit a critical condition at the 
crest (e), and therefore, a subcritical condition 
for some distance to the lee. Since a smooth 
transition to the downstream supercritical 
condition is not possible, a discontinuity of 
flow must occur somewhere. We may have 
a jump upstream but it does not seem possible 
to exclude a hydraulic drop in the lee on the 
basis of the tools developed so far for analyzing 
the model. If the principle of conservation of 
momentum is employed, the drop is seen to 
require a negative dissipation of energy and 
is thus impossible. This resort is not available 
for multiple-fluid systems and in the subsequent 


* As in fig. 17 this wave of elevation may break 
forming an upstream bore. 
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sections we must employ experimental evi- 
dence to indicate which of these possibilities 
will be realized. 

The last possibility (f) for the one-fluid 
system will be obtained if the initial speed of 
the obstacle is faster than that required for (e). 
The flow will then be supercritical everywhere. 

It is of interest to consider the form of the 
free surface in the various cases of Table 1. In 
case (a) there will be a slight draw-down over 
the obstacle returning to the original level 
downstream (fig. 15). In (b) and (c) the draw- 
down over the obstacle is followed by a lee 
jump and an abrupt rise to the undisturbed 
downstream level (fig. 16). In cases (d) and 
(e) the upstream jump is followed by a gradual 
lowering of the surface toward the crest (fig. 
17). In (f) the surface will rise smoothly over 
the obstacle. Such a phenomenon is called a 
“swell” in hydraulics. 

Given the undisturbed fluid level, H, the 
height of the obstacle, b,, and the speed of the 
approaching current, U, it should be possible 
to predict which of the possibilities of Table 1 
will be realized. Obviously, the deciding factor 
is the condition that the flow be critical at the 
summit. Using this, plus the energy and con- 
tinuity equations we have, at the crest, 


2 2 
Us u WA (9) 

28 29g 

ue 
gh, I, (10) 
Dh D. (11) 
Eliminating u. and h,, 
: b. 

F?2-3F’s=2 (- 1) ae NO He 


A graph of this equation is shown in fig. 18. 
The regimes are labeled in accordance with 
the system employed in Table 1. 


4. Theory of a two-fluid system 


In the case of a two-fluid system we assume 
that the density difference between the two 
layers, Ao/o, is very small. The approach 
velocities are assumed to be the same in both 
fluids. It is not difficult to see intuitively that 
velocities accompanying marked effects at the 
interface will be so small as to produce little 
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Fig. 18. Critical curve of a one-fluid system. Compare 
Table ı. 


disturbance at the free surface. Kinematically, 
therefore, it should be possible to replace the 
free surface with a rigid surface. If this is done 
it is necessary to consider pressure variations 
(x) on this upper rigid surface, just as the 
dynamic effects of the very slight changes in 
height of the free surface would also be 
important. In the following development, 
therefore, the model of fig. 19 is adopted in 
which the two-fluid system is capped by a 
rigid surface. À parallel development with a 
free surface yields identical results provided 
conditions are such that the flow is well 
within the absolutely subcritical range with 
respect to deformations of the free surface. 
This puts generously high upper limits on the 
size of the obstacle and the internal Froude 
number. These limits may be computed from 
fig. 18. 

With the hydrostatic assumption the energy 
equations in this model are 

2 


x + o'gH+ 0 = constant, ‘ee, 
a —o'g(h+b-— H)+og(h+b)+ a = constant. 
(14) 


The total energy transport through a given 
section is 


if PE re u? 
E,=(Q+ Q') 7+ QVe'gH+ Qe + Qo—+ 


+ Qog(h+b) (15) 
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Fig. 19. Theoretical model of two-fluid flow over a 
barrier. 


- Qo'g(h+b- H)= Es = constant. 


Evaluating E, upstream and substituting for 
rt from equation (13), 


LE (Hh)? (H-h-5)-2]+ 


2 
Ze nee hy #] + (e— 0')g (h— ho + D) =0. 
(16) 
If we define the quantities 
2 _ 20 _ lo af aay 
SACS torre Set Sra er 
(17) 


and assume the same approach velocity, U, in 
both layers, equation (16) becomes 


FF Rÿ CRE 
AE  (1-&-ß)? 


This relationship gives & as a function of ß, 
i.e. the profile of the interface during passage 
over the obstacle. Curves of this equation for 
R, of 0.30 and 0.70 are shown in figs 20 and 21. 

(a) Ry < 1/2. Fig. 20 is typical of flows with 
a lower layer less than one-half the total depth. 
If F; exceeds 0.50, the interface will swell over 
a barrier returning to its original depth on the 
lee. This theoretical behavior is confirmed by 
experiment, as for example in fig. 7. If F; <.50, 
the interface will tend to lower slightly as the 
fluid ascends the barrier. For example, if 
F; = .10 and the obstacle has a maximum 
height less than .155, the interface lowers to 
a minimum at the summit, returning smoothly 
to Ry on the lee. This is the situation in fig. 2. 


|+=+#-Ro=o. (18) 
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Fig. 20. Thickness of lower layer, «(x), during passage of 
two-Auid system over a barrier of height, B(x). The 
upstream depth of the lower fluid is Ry = .30. 


If the barrier is precisely .155, the interface 
will continue to lower on the lee. From the 
one-fluid analogy we would expect a very 
small (theoretically infinitesimal) hydraulic 
jump just on the lee, the interface returning 
to a downstream level slightly less than Ry 
because of the energy loss. If the obstacle is 
appreciably greater than .155, a steady state is 
not possible. This causes a change in R,, i. e. 
a blocking effect takes place that increases Ro 
to a level that permits a steady state between 
some upstream section and the crest of the 
obstacle. Thus, if the obstacle has a maximum 
height .30 and is held at a speed corresponding 
to F; = .10, additional curves of equation (18) 
show that the upstream depth of the lower 
fluid will tend to increase from .30 to .50. The 
system will then have attained a quasi-steady 
state upstream, the interface will lower as the 
fluid ascends the barrier and will continue to 
lower on the lee. Since the blocking effect has 
created a differential in the upstream and down- 
stream levels, the interface can take a finite 
jump on the lee to the lower energy level 
downstream, as in figs 4 and 6. We will see 
later how we can predict the new value of Ro, 
knowing the approach Froude number and the 
height of the obstacle. 

If the obstacle is higher than .30, in the case 
of fig. 20, and F; < .50, the barrier blocks the 
lower fluid completely until it builds up to 
obstacle height, spills over and forms a hydrau- 
lic jump on the lee as in figs to and 11. If 
F; > .50 it may be expected from the curves 
that the interface will swell over the barrier, 
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Fig. 21. Thickness of lower layer, «(x), during passage of 
two-fluid system over a barrier of height, B(x). The 
upstream depth of the lower fluid is Ry = .70. 


returning smoothly to its equilibrium height. 
This is confirmed by experiments, however, 
only if the height of the barrier is not too much 
greater than Ro. In fig. 9, however, the inter- 
face does not return smoothly on the lee but 
appears to drop down with considerab]e tur- 
bulence. A possible explanation is the f„jlow- 
ing: In fig. 9 we observe that the in,.rface 
traces a curve close to the F; = .60 cu ve in 
fig. 20 up to the crest, B = .50. This ist to be 
expected since there should be very little loss 
of energy in this region. As the fluid attempts 
to move down the lee of the obstacle, bounda- 
ry-layer separation and resulting turbulence 
will abstract energy. Fig. 9 shows that the 
interface stays at its maximum height for a 
considerable distance in the lee which means 
that the interface describes a line, « + B= 
const .70. In an intuitive sense we may 
assume that the loss of energy to turbulence, 
throws the fluid system from an energy level, 
corresponding to F; = .60 to a lower energy 
curve close to that of the F; = .40 curve in 
fig. 20. This, however, would require a deep 
lower fluid on the lee (.70), whereas the 
experimental set-up requires the fluid to return 
to the vicinity of .30. This can be, and ap- 
parently is, accomplished by the occurrence of 
a shock wave, in this case a hydraulic drop, 
to form a connection with the downstream 
level. According to this intuitive argument, 
this phenomenon requires a nearby curve of 
lower energy level at the point of the diagram 
corresponding to the crest of the barrier, so 
that a modest loss of energy to turbulence will 
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permit a return of the interface height along 
an entirely different path than the approach. 
This requirement is met for a large -barrier in 
this case. If the barrier has a maximum height 
of, say, only .15, the F; lines at the crest are 
fanned out sufficiently so that there is no nearby 
curve and the interface returns to close to the 
equilibrium height downstream along much 
the same curve as it described on the upwind 
side. It is evident from fig. 20 that the higher 
the barrier the more easily is this shift to a 
different branch of the F; curves accomplished. 

The above argument lacks rigor because the 
various F; curves in fig. 20 do not correspond 
to different energy levels for given approach 
conditions but rather to different approach 
conditions and correspondingly differentenergy 
levels. It does not seem possible to verify the 
argument at this time because the present 
theory is inadequate to permit a calculation of 
alternate energy curves for the given approach 
conditions, Ro and F;. 

(b) R,> 1/2. It has been stated that the 
behavior of a system in which Ry > 1/2, is 
much the same as for shallower lower fluids 
except that for large or moderately large 
obstacles and high Froude numbers the sym- 
metrical regime is missing and a hydraulic drop 
occurs in the lee. The suggestion above for the 
occurrence of a drop may also be applied in 
this case. We see from fig. 21 that the F; 
curves are very tightly packed in the upper 
part of the diagram. Thus, for example, in 
fig. 12 the curve traced by the interface should 
follow the F; = .542 curve until the crest. The 
energy losses to turbulence in the lee will 
permit the interface to describe a curve, 
a + B = .90 approximately. The return to a 
height somewhere near the equilibrium, 
Ry = .75, downstream would require such a 
hydraulic drop as appears in fig. 12. 

It would be desirable to have curves analo- 
gous to fig. 18 for a two-fluid system. Such 
curves are obtained by assuming, as in the 
single-fluid case, that conditions are just 
critical at the top of the barrier. Thus, if we 
differentiate (18) with respect to x, 


A Re RE Nil 
tlt goo le) 
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At the top of the obstacle, therefore, either 
the slope of the interface vanishes or 


F2 Ey a, 


a3 I-a,- B,)8 


C 
where the subscript c indicates that the 
section is at the top of the barrier. If there has 
been no energy loss upstream from the barrier, 
then at the top we have also from (18) 


“(3 (1 — Ro)? 
2% oc ee 


] = I, (20) 


] act Be Ro=o. 
(21) 


Between equations (20) and (21) we may 
eliminate «, and, for a given Rg, obtain a 
relationship between F; and £, when the con- 
ditions at the crest of the barrier are just 
critical. It is obvious from the considerations 
of the previous section that this relationship is 
of fundamental interest. For a given Ry the 
resulting curve will be analogous to the curve 
of fig. 18 for a single fluid. These curves are 
given in fig. 22 for steps of 1/10 in Ry. They 
resemble in general the curve of fig. 18 except 
that the upper branches differ. It follows from 
the preceding discussion that, for a given 
approach depth, Ro, of the lower fluid, points 
on or below the appropriate Ry curve in fig. 22 
represent possible quasi-steady states. If the 
given obstacle height and Froude number 
correspond to a point below the curve, the 
conditions are absolutely subcritical and no 
shock wave will occur. If the point falls on the 
curve a lee jump occurs (or lee waves for very 
small obstacles). If the point is somewhat above 
the curve the upstream depth should increase 
to the R, of the curve passing through this 
point. If F; is sufficiently large either the 
symmetrical, absolutely supercritical regime 
will occur or a hydraulic drop will be ex- 
perienced in the lee. 

In the case of a single fluid it is well known 
that a finite amplitude long wave tends to 
“break” with accompanying turbulence (Lams, 
1932). This follows from the consideration that 
the speed of propagation of a given height of 
the wave is proportional to the height. The 
high portions then eventually catch up with the 
troughs and the wave configuration becomes 
unstable. Since this is intimately connected 
with the hydraulic jump phenomenon, it is of 
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considerable interest to investigate the behavior 
of finite long waves in a two-layer system. In 
particular it is of interest to investigate the 
hydraulic drop observed in the experiments 
under certain conditions when the upper fluid 
becomes more shallow than the lower. 

It is easily shown that the speed, co, relative 
to some coordinate system, of an infinitesimal, 
long wave at the interface is given by 


(10 — Co)? | (Ho — Co)” 
"h ae 'h 
KALT) & "20 
where g’ = gAo/o, and, for convenience, we 


denote the lower layer by subscript ı and the 
upper by subscript 2. Solving for cg, 


=f. (22) 


_ Myohgg + Maglio _ 


à = | 
2 
hy ohoo hyoltoo 2 aa se 
= [ ii rm (#10 — U20) ] » 


where we confine our attention to the wave 
tending to move upstream (toward negative 
x). Let us now consider a long wave of finite 
amplitude in which the height of the interface 
at some section is say hy + 7. Since the inter- 
face has an infinitesimal slope in a long wave, 
the interface is nearly horizontal in the vicinity 
of this section and, moreover, nearby portions 
of the interface are at infinitesimally small 


Fig. 22. Critical curves 
of a two-fluid system 


100 


distances above or below hy + 7. The speed 
of propagation of nearby wave heights should 
therefore be given by equation (23) if we 
substitute the heights and velocities in the 
vicinity of this section. Thus 


Ba (hep — 7) + Hy (hy +) 
H 
g (hyo +n)(hao = n) 
| a à \ (24) 
= (hyo + 9) (hoy — 7) (4 — aoe | 
) 


H? 


In order to obtain the speed of propagation in 
terms of the approach speed, U, of the two 
fluids we use the continuity equations at the 
section in question, namely 


d oe 
= (hyp +) + 7 Un (ha#n)]=0,. (2%) 
and 


0 à 
A (hog — 4) + Be [ua (he — n)]= 0. (26) 


v 


In view of the quasi-steady state the time 
derivatives may be expressed as 


0 0 


Ta (27) 
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If we confine our attention to values of n/H, 
small compared to one, then to the order of 
n/H, the c in equation (27) may be replaced 
by the & of equation (23). Using this and 
integrating equations (25) and (26), 

Ua (hyo + 9) — con = Uhr, (28) 
and 

2 (hag — 9) + con = Uhso. (29) 
Substituting for u, and u, in equation (24) and 
retaining terms only of the order of n/H, 


MT g’hyohzo at 3 (1-2R,) n 
c=U (Sree ) I+ 


PRR) He 
(30) 
where R is hjo/H. Therefore, the speed of 
propagation of the elevation 7, relative to the 
equilibrium speed U of the fluids, is 


Shy ohoo ail TITRE 
( H ) MEN ATEN AT ES 


a result analogous to that of Airy for a single 
fluid (LAMB, 1932). This result shows that a 
crest moves faster than a trough provided the 
equilibrium depth of the lower fluid is less 
than one-half the total depth. On such an 
interface one would expect the lower fluid to 
thicken through a hydraulic jump as in a one- 
fluid system. If the upper fluid is deeper than 
the lower, however, a trough will move more 
quickly than a ridge and the resulting shock 
wave should take the form of a hydraulic drop. 
This suggests roughly why hydraulic drops 
occur in the experiments when the upper fluid 
becomes shallow. The value of the analysis has 
two severe limitations, however: First, the 
wave while finite, was considered small; 
second no account was taken of the shear 
between the two fluids arising in passage over 
the barrier. A full discussion of the initial value 
problem of a finite gravity wave in a two-fluid 
system will appear in a fortcoming paper. This 
independent approach yields results in agree- 
ment with equation (31) to the order of 7/H. 


5. Meteorological implications 


It was pointed out in Part I that the eventual 
goal of the investigations of this series of 
papers was an understanding of small-scale 
atmospheric flow over mountains and hills. 
This has been sidetracked somewhat by the 
complexity of the general problem of statified 
fluids and it has become necessary, as in the 
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present paper, to explore first extremely simple 
systems before attempting to simulate the 
vertical distribution of wind and density in the 
atmosphere, complex terrain structures, effect 
of the earth’s rotation, etc. The outstanding 
questions of similarity may be listed as follows: 

(a) What is the effect of neglecting the 
earth’s rotation? This question has been dis- 
cussed briefly by Lone (1953 a, 1953 b). Ob- 
viously rotation may be neglected if the 
horizontal dimensions of the hill or mountain 
are small enough. If we require that the 
Rossby number be large compared to one, 
for example, rotation may be neglected 
provided the horizontal dimension of a ridge in 
the direction of the current is to km or less. 
This appears to be too severe a requirement, 
however, as is apparent from the following 
considerations: If a two-fluid system flows 
over a ridge in a rotating system, the rotation 
will not be important in determining the 
regimes of motion and critical Froude numbers 
if the kinetic energy of the velocity component 
parallel to the range is small compared to the 
kinetic energy of the basic motion, i.e. if 
v?/u? < 1. Considering the lower fluid only, 
at the crest v will be of the order of fA/ho, 
where f is the Coriolis parameter, A is the 
cross-sectional area of the ridge from base to 
crest, and hy is the height of the inversion 
at the upwind base of the ridge. Since, for 
the atmosphere, u? = o(105 cm? sec?) and 
A = 0 (bL/4), where b is the maximum height 
of the ridge and L is the total horizontal length, 
we must have L? < 16 x 10® x 108 h5/b? cm?. 
If we take a fairly typical value, h,/b = 2, for 
example, we have L<8x10’, so that the 
ridge may be as wide as 100 km before the 
effect of rotation is appreciable upwind of the 
crest. If there is a tendency for supercritical 
flow on the lee, the lower layer may thin con- 
siderably and-higher v-velocities generated. 
This will have an influence on the nature of 
lee waves or lee jumps rather than on the 
criteria for the various regimes of flow. 

(b) What is the effect of the compressibility 
of air? Since the model employs incompressible 
liquids, a question arises as to the effect of 
compressibility. In a forthcoming study of this 
problem it appears that the atmosphere may 
be considered incompressible if we substitute 
potential temperature gradients (40/0) for 
density gradients and if the ratio gb/c? < 1 
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where b is the height of the barrier and c is 
the speed of sound. This requires b < 10%, 
which is satisfied by mountains 1 km high or 
less. Again this is probably too severe a 
restriction. A computation employing the 
methods of section 4, but with two compres- 
sible fluids, shows very little effect of com- 
pressibility for barriers several times larger 
than this. 

(c) To what extent do inversion surfaces in 
the atmosphere resemble the interfaces and free 
surface in the models? In most of the past 
investigations of flow over mountains it has 
been assumed that the atmosphere has a more 
or less uniform potential temperature gradient. 
In such a case no vertical height dimension, 
short of the total depth of the atmosphere, 
enters into the analysis. This is at complete 
variance with the consideration of this paper 
in which heights of inversions, particularly the 
tropopause, are used to define Froude numbers 
whose critical values are decisive for the 
occurrence of the laminar flows, jumps, and 
blocking. This important difference cannot be 
resolved at the present time, but future experi- 
ments should be able to provide a rather 
decisive answer. It is possible, for example, by 
the use of water and salt mixtures and/or layers 
of immiscible fluids, to reproduce faithfully 
typical potential temperature distributions in 
the atmosphere up to levels (perhaps 80,000 ft) 
where we may feel confident any effects of 
topography disappear. 

In conclusion, let us assume for the purpose 
of discussion that the model of this paper is 
essentially similar to the atmosphere, and that 
the free surface corresponds to the tropopause 
and the interface to a tropospheric inversion 
surface. What suggestions regarding the flow 
of air over barriers are contained in the results 
of this paper? In the first place we may say 
that, unless the tropospheric inversion is quite 
high and the barrier quite small, the main 
phenomenon of interest is not lee waves of the 
classical type but rather shock waves resembling 
hydraulic jumps, surges or bores. This results 
from the fact that the pressure distribution in 
the model is closely hydrostatic except in the 
vicinity of the jump. The hydrostatic assump- 
tion should be at least as good in the atmos- 
phere, so that the shock phenomenon should 
be typical of air flow over most mountains and 
hills, provided inversions are present. Lee 
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waves of the classical type should be confined 
to hills with maximum heights considerably 
less than the height of the first inversion surface. 

Lee jumps should form in the lee if the 
internal Froude number exceeds a certain 
value, depending on the height of the inversion 
and the height of the obstacle. For lower 
tropospheric inversions the critical value of 
F; is of the order of .20. In the atmosphere we 
usually have a shear of wind with height which 
is not reproduced in the experiment. If, how- 
ever, we use the mean velocity of the tropo- 
spheric air, for the purposes of the Froude 
number, and take (gH A@/O)=:2 104, jump 
phenomena require winds of the order of 20 
m sec-!. On the basis of this computation we 
would expect hydraulic jump phenomena to 
be a common feature of flow over the larger 
surface barriers. Such phenomena as absolutely 
supercritical flow and hydraulic drops probably 
do not occur since they require mean winds 
of the order of so m sec-t. 

It is of some interest to speculate about one 
possible effect of rotation on the flow of air 
over a range as large as the Appalachians. If, 
for example, a quasi-steady jump would form 
on the castern side, by virtue of flow over the 
range with a tropospheric inversion, it is 
conceivable that the resulting convergence on 
the lower layer in the jump could lead to 
sizeable increases of relative vertical vorticity 
and a possible source of cyclogenesis. Such an 
effect would require that columns of air in the 
lower layer would lose to the ground by 
friction much of the relative negative vorticity 
generated by the vertical shrinking of the 
columns in passage over the mountains. It is 
suggestive that the Cape Hatteras region is a 
favorite site for cyclogenesis. An exploratory 
experiment with a barrier moved in a rotating 
system of two fluids, did not show any such 
effect. This should not be regarded as conclu- 
sive, however, because the frictional action 
between the barrier and the lower layer of 
fluid may be entirely different in the experi- 
ment and in the atmosphere. 
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Abstract 


Three sections across the Gulf Stream in which deep temperature and salinity observations 
were made close intervals show that comparably large values of surface current velocity are 
obtained from at dynamic computations, the Geomagnetic-Electrokinetograph (G.E.K.) and 
from navigational data. The dynamic computations indicate that the crosscurrent slope of the 
isobaric surfaces is uneven, resulting theoretically in three or four zones of high velocity separated 
by zones of lower velocity. The navigational and GEK measurements do not clearly reflect 
this picture. It is suggested that the unevenness of the slopes of the isobaric surfaces is related 
to the surface temperature discontinuities observed in the Gulf Stream by von Arx and RıcH- 
ARDSON (1953) from aircraft, and that the Stream is composed of overlapping, discontinuous 


currents at all levels. 


Introduction 


With the advent of Loran and GEK equip- 
ment it became possible to measure the velocity 
of surface currents from a moving ship with 
some accuracy. In the Gulf Stream velocities 
of 4 to 5 knots (200-250 cm/sec) were 
normally found by Loran in the swiftest part 
of the current (ISELIN and FUGLISTER 1948). 
Von Arx (1950) measured a peak velocity of 
274 cm/sec with the GEK and values of more 
than 200 cm/sec are commonly observed with 
this instrument. 

These velocities were in excess of those 
previously found by Dynamic Computations, 
a previous maximum being 140 cm/sec 
(SEIWELL 1939, Fig. 19) in a typical section 
from Montauk Point to Bermuda. ISELIN and 
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FUGLISTER (1948) pointed out that values 
obtained by means of the geostrophic equation 
in the period 1931 — 1940 must be regarded as 
minimum values since the stations were 
usually too widely spaced to observe the true 
slopes of the isobaric surfaces. Also, with the 
navigational equipment available during that 
period it was impossible to tell at what angle 
to the Stream a section was made. If this angle 
were to deviate more than about 10° from the 
perpendicular a noticeable reduction in com- 
puted velocity would result, since the com- 
putations are concerned only with the com- 
ponent at right angles to the sections. 

A cruise was made in the ATLANTIS in 
October and November 1950 in which the 
primary objective was to make sections across 
the Stream in which the slopes of the isobaric 
surfaces would be closely observed. This 
involved making hydrographic stations, con- 
sisting of serial observations of temperature 
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Fig. 1. Temperature at 200 m during preliminary survey and section 1. 


and salinity, at closely spaced intervals across 
the current and at right angles to its axis. 
Without Loran, which enables the observer to 
fix the ship’s position as frequently as he 
desires, this program could not have been 
attempted. It was planned to measure the 
surface current between stations with the GEK 
and Loran and to arrive at a proper basis for 
comparing the values obtained by these in- 
struments with those computed from the 
temperature and salinity data. 


The Sections 


A preliminary survey was made to establish 
the location and curvature of the Stream. The 
results of this survey and the placings of the 
first two sections are shown in Figure 1. The 
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contours of 200 m temperature are usually 
considered an adequate index of the location 
of the Stream (FUGLISTER 1951). The placing 
of Section 1, stations 4842 — 4853, was decided 
on the basis of this survey. Section 2, stations 
4853 — 4867, was made immediately afterwards 
on a reversed course since negligible changes 
in the position of the Stream had taken place 
between the time of the survey and the time 
when Section 1 was completed. 

Four days elapsed before Section 3 was 
made and no time was available for another 
survey. Consequently the placing of Section 
3 was decided on the basis of a single crossing 
of Stream. This resulted in Section 3’s not 
being made at right angles to the current since 
the Stream had meandered during the four 
days in a manner not recognized by the 


118 Le NEW ORTE ION 
ie} 
S a a 
70° = a a = 
DATE + 
39° 1102. 30 OCT. (1950) 8 = = 
Le] © 
o Se MO a os + 
1600 3 NOV mes BS Bass À à = & 8 
sı 49 47 
4853 52 50 48 46 45 44 43 4842 
38° 
no 
c 
Ww 
+ 
w 
= 
= 
= 
+ 
a 
w 
a 
2056 2 NOV sae 
ca 0900| 29 ocr. 
Fig. 2. Temperature at 200 m before, during and 
after section 3. 


observers. Figure 2 shows the location of 
Section 3, the crossing (29-30 Oct.), on the 
basis of which the location was chosen, and 
the contours of 200 m temperature before, 
during and after the section. GEK current 
direction data indicate that Section ı deviated 
2° from the perpendicular to the axis of the 
Stream, Section 2 four degrees and Section 3 
sixteen degrees. (The axis of the Stream was 
considered to be the average direction of those 
GEK current fixes with velocities of more than 
100 cm/sec in each section.) 

The temperature profiles of the three sections 
are shown in Figures 3, 4, and s, together with 
the distribution of the data, each dot represent- 
ing an observation of temperature and salinity. 
Salinity profiles are not presented here since 
the T/S relationship was consistent, except in 
very shallow (0 —-200 m) observations close to 
the inshore edge of the current. 

In general two series of eleven Nansen 
bottles each used, the shallow series from ı to 
900 m and the deep from 1,000 to 3,200 m by 
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Fig. 3. Section 1:Temperature. 


wire. Large wire angles were encountered, 
particularly in Section 1 where no efforts were 
made to reduce them by maneuvering the ship. 
In the second and third sections, where the 
ship was maneuvered, the bottles reached a 
depth of more than 2,000 m on all of the 
stations but one. Five unprotected thermo- 
meters were usually used with each series and 
much confidence is placed in the accuracy of 
the bottle-depth determinations. 

The station positions given in Figures 1 and 
2, and the times given in Figures 3, 4, and 5 
are those of the shallow series which usually 
reached a depth of 600-800 m. A Loran fix 
was obtained at the time the messenger was 
sent on each series. These sections are the most 
closely spaced available in the Gulf Stream, 
the usual distance between stations in the past 
having been about 20 nautical miles. 
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Fig. 4. Section 2:Temperature. 


The Dynamic Computations 


In all important particulars, the method of 
computation described by La FOND, (1951) 
has been used. The major source of inaccuracy 
in dynamic computations lies in the fact that 
small errors in position may give erroneously 
steep or gentle slopes to the isobaric surfaces, 
resulting in large errors in the computed 
currents. VON ARX (1952) has estimated that 
Loran fixes have an average error of 1 km. 
(about 0.5 nautical miles) in this area. Since 
Loran fixes were usually obtained at half- 
hourly intervals during these sections this error 
should give a scatter of + 1 knot (about 100 
cm/sec) in the values of current velocity 
obtained by this means. The observed scatter 
was considerably smaller than this which leads 
one to suspect that von Arx’s estimate is pessi- 
mistic. However, in presenting the computed 
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Fig. 5. Section 3:Temperature. 


surface currents graphically (Figs. 6, 7, and 
8) the effect of an 0.5 nautical miles error in 
the position of each station has been plotted 
as an envelope around the curves of surface 
velocity derived from the computations. 


The true distances between stations were not 
used in all cases since many stations fall up or 
downstream from a mean line drawn through 
them. Instead it was assumed that all stations 
fall on the mean line. 


The computed velocities are based on the 
assumption of a level of no motion at 2,000 m. 
This level was chosen because nearly all the 
stations extended to that depth, and the errors 
involved in the choice of this level are small. 
In Sections 2 and 3 data were available down 
to 3,000 m. These data show that a 5 % 
increase in computed volume results if the 
3,000 decibar surface is assumed to be motion- 
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Fig. 6. Surtace velocity—section 1. 


less instead of the 2,000 decibar surface. Calcu- 
lated velocities at 2,000 m referred to the 3,000 
decibar surface average, 3.5 cm/sec in Sections 
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Fig. 7. Surface velocity—section 2. 


2 and 3. This figure is approximate since the 
errors of salinity analysis (reckoned at + 0.02 
%) are critical at these depths. 


Current Measurements by GEK 


Current fixes were obtained with the GEK 
as frequently as time allowed when steaming 
between stations. The measured currents have 
been increased by a K factor of 1.2 (VON Arx 
1950). This K factor is an empirical factor 
based on all the comparisons between currents 
measured with Loran and those measured with 
the GEK on this cruise, in which the GEK on 
average gives lower values than Loran. 

More GEK measurements would have been 
desirable but they were sacrificed, perhaps too 
freely, in order that the ideal synoptic section 
might be more closely approached. About 25 
measurements were made with each section. 


Current Measurements by Loran 


Loran fixes were obtained at half-hourly 
intervals when underway and the currents 
calculated according to the ship’s displacement 
from its dead reckoning positions. The effect 
of wind drift was ignored. The errors in the 
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Loran positions have been discussed above. 
About 35 current measurements were made 
with Loran in each section, and the component 
values at right angles to the sections used for 
comparison with the dynamic values. 


The Comparisons 


The three surface velocity profiles show that 
comparable values are found with all three 
methods of measurement, but the dynamic 
calculations indicate an unevenness of structure 
that is not clearly shown by the other two 
methods. To be sure, the Loran gives a more 
uneven picture if all the positions are assumed 
to be exact but the scatter of the Loran points 
is ascribed here, as elsewhere, to unavoidable 
errors in the system. 

In all three sections the geostrophic equation 
consistently argues that the Stream is compose 
of bands of swift moving water separated by 
bands of relatively slow moving water. The 
navigational data are sufficiently accurate to 
ensure that large changes in the cross current 
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slopes of the isobaric surfaces are real. What 
is not clear is whether or not there is a suffi- 
ciently close relationship between the slopes 
of the isobaric surfaces and the surface currents 
to warrant the assumption that the cross- 
current profile presented by the dynamic com- 
putations is more correct than those presented 
by Loran and GEK. 

The dynamic computations ignore the effect 
of frictional forces which would tend to 
smooth out the unevenness required for a 
steady state. It is further possible that local 
winds and waves have a like effect on the 
surface layer. On the other hand the high 
scatter of the Loran measurements may obscure 
unevennesses in the currents that actually exist, 
and recent work by Longuet-Higgins, Stom- 
mel and Stern (not yet published) indicates 
that GEK values are influenced by the motion, 
or lack of motion, of water at some distance 
from the point of observation. The effect of 
this would certainly be toward averaging out 
the high and low velocity zones in these 
sections. 

It is worthy of note that the existence of 
such high and low velocity zones was not 
suspected by the observers in the ship even 
though continuous temperature profiles were 
kept up as the station data came in. The 
changes in the slopes of the isotherms were 
regarded as negligible except in cases near the 
northern edge of the Stream where the slopes 
were actually reversed. 

In retrospect, looking through past bathy- 
thermograph profiles across the stream where 
the 900 foot instrument was used, it was scen 
that changes in the slopes of the isotherms at 
the 200 m level, comparable to the changes 
found in these sections, are the rule rather than 
the exception. These bathythermograph data 
were obtained from ships continuously under 
way, which implies that the changes in slope 
found in these three sections cannot be attri- 
buted solely to the long time lapse between 
stations, and that comparable changes would 
have been found if all the stations in a section 
had been occupied simultaneously. 


Subsurface Currents 


From the computed velocity profiles illus- 
trated in Figures 9, 10, and 11 it can be seen 
that the unevenness of the Stream is not 
confined to the surface layers. If the unevenness 
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Fig. 10. Section 2—velocity cm/sec. 


were due to errors in navigation one would 
expect that a high or low velocity zone would 
be found at all depths between the same pair 
of stations. In Figure 11 dashed lines have been 
drawn through the high and low zones, and 
it can be seen in many cases that a high zone 
at the surface is over a low zone at mid-depths. 


In Section 3 (Fig. II) the Stream appears to 
be divided into three parts, the two major parts 
being separated by a deep countercurrent. It is 
tempting to conclude that here the Gulf 
Stream is moving laterally across the counter- 
current by re-forming to the north, and that 
the low velocity zones found further towards 
the southern edge of the Stream are the rem- 
nants of countercurrents which have been 
absorbed into the main current in the same 
manner, further upstream. Such crossings do 
not imply that Slope water is moving across 
the Stream but that changes in the slopes of 
the isobaric surfaces move slowly across the 
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Fig. 11. Section 3—velocity cm/sec. 


Stream in a manner analogous to 
motion. 


Recent studies using aircraft (VON Arx and 
RICHARDSON 1953) have shown that on the 
surface the Florida Current and Gulf Stream 
fronts are discontinuous. Using the Stommel- 
Parson airborne radiation thermometer they 
found that surface temperature gradients 
associated with the front diminished in a 
downstream direction to the point where they 
could no longer be detected, but that by 
turning to the left a fresh gradient could be 
found which in its turn would fade out. These 
three cross-sections suggest that these discon- 
tinuities are not merely surface phenomena but 
that the Gulf Stream is discontinuous at all 
levels all the way across its width. More 
accurate direct current measurements are 
necessary to determine to what extent the 
velocity pattern is a part of these discontinuities 
in the pressure field. 


Wave 


Tellus VI (1954), 2 


CROSS-SECTIONS OF THE GULF STREAM 123 


REFERENCES 
FUGLISTER, F. C., 1951: Multiple Currents in the Gulf von Arx, Wırıam S., 1950: An Electromagnetic 
Stream System. Tellus 3 (4), 230 - 233. Method for Measuring the Velocities of Ocean 
Iserın, C. O’D., and F. C. FUGLISTER, 1948: Some Recent Currents from a Ship Under Way. Pap. Phys. Oce- 
Developments in the Study of the Gulf Stream. J. anogr. Meteorol. IL (3), 1-62. 
Mar. Res. 7 (3), 317 — 329. — 1952: Notes on the Surface Velocity Profile and 
La Fonp, E. C., 1951: Processing Oceanographic Data. Horizontal Shear across the Width of the Gulf 
U. S. Navy Hydrographic Office. H. ©. Pub. No. 614. Stream. Tellus 4 (3), 211 - 214. 
Washington, I -II4. — and WırLıam S. RICHARDSON, 1953: The Surface 
SEIWELL, H. R., 1939: The Effect of Short Period Varia- Outcrop of the Gulf Stream Front. W.H.O.I. Ref. 
tions of Temperature and Salinity on Calculations 53 —24, April 1953. 1-4. (Unpublished manuscript.) 


in Dynamic Oceanography. Pap. Phys. Oceanogr. 
Meteorol. 7 (3), I —32. 


Tellus VI (1954), 2 


Observations of Lee Wave Clouds in the Jämtland Mountains, 
Sweden 


By LENNART LARSSON, Weather Service of the Swedish Air Force. 


(Manuscript received 4 February, 1954) 


Abstract 


An account is presented of three years of observations of lee-wave clouds in the Jämtland 
mountains, Sweden, with detailed descriptions of some individual cases. Wind and stability 
conditions in connection with the occurrence of the clouds are given. A brief discussion of the 
agreement between observed conditions and theoretical results is attempted. 


I. Introduction 


An air-current passing over a mountain 
range is subjected to various kinds of deforma- 
tions as described in an account by QUENEY 
(1948). One of the effects often occurring even 
in comparatively small mountains is the devel- 
opment of a definite wave-pattern which 
usually is stationary in relation to the ground. 
Its existence and characteristics are determined 
by the vertical distribution of the wind and 
the temperature, the structure of the obstacle 
and its position relative to the wind. If the 
humidity is favourable -not too high - these 
waves may become visible through the devel- 
opment of lenticularis-type clouds. Under 
certain weather conditions the cloud-pattern 
over many mountain areas is characterized by 
such stationary wave-clouds. The air-flow over 
Riesengebirge and the “Moazagotl-clouds” 
formed there have been described in detail by 
KUTTNER (1939). Other areas where wave- 
clouds are common can be found in Greenland, 
Iceland, Crossfell in England and the Italian 
and Swiss alps. Perhaps the most magnificent 
example of this phenomenon is provided by 
the “Sierra Wave” or the “Bishop Wave” in 
California (DE VER CoLson, 1952). 


Glider-pilots have actively contributed to the 
discovery and study of these wave-upwinds 
and ever since their discovery in about 1930, 
they have utilized them for sportflight. In con- 
nection with the exploration of the wave- 
upwinds of the Sierra Nevada regular routine 
flights were made with gliders to altitudes of 
10,000 — 12,000 metres. 


2. Methods of observation 


Observations have been conducted daily for 

a period of three years, October 1950 — Sep- 
tember 1953 and were mainly concentrated to 
the central part of the Scandinavian mountain 
range shown in fig. 1. The writer has made 
observations partly from Frösön where a 
Swedish Air Force wing is stationed and from 
where a free view is obtained of the greater 
part of the area shown in fig. 1, and partly 
from the air in connection with daily weather 
reconnaisance flights over this area. In this 
manner and through the aid of continuous 
reports from specially briefed pilots it has been 
possible to form a satisfactory picture of when 
and where lee-wave clouds appear even at 
times when low clouds have prevented ob- 
servations from the ground. In most cases it 
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Fig. 1. Topographic map of the eastern side of the central part of the Scandinavian mountain range. The thick broken 
line on the left is the Norwegian Border. À indicates the observation place at Frösön. The elliptical figures drawn 


with thick solid lines indicate the positions and the hor 


izontal dimensions of piles of wave clouds. C, D, E, F and 


G also indicate positions of piles of lee wave clouds on different occasions (cf. text). 


has been possible to measure or estimate the 
altitudes of the clouds, their size, position and 
persistance and to some extent structure. No 
instruments specially designed for these pur- 
poses were available in the air plane. Some 
time-lapse-filming was also undertaken. Careful 
notes of the observations were made. Data on 
wind, temperature and humidity distribution 
were obtained from daily radio-soundings at 
Frösön, 0300 and 1500 GMT. Pibals were 
made at 0300, 0900 and 1500 GMT or at times 
when wave-clouds had appeared. The wind 
velocity measurements at 0300 and 1500 GMT 
were obtained by observing the radio-sonde 
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balloon through theodolites and the wind data 
thus obtained were as accurate as pibal ob- 
servations by two theodolites. The wind data 
were determined in this manner in all cases 
discussed below and in most other cases. 

The topographical conditions can be seen 
on fig. 1 which shows the eastern side of the 
central part of the Scandinavian mountain 
range. It is rather irregular but a fairly con- 
tinuous mountain range extends from the 
mountains at the Norwegian border in the 
west to the Ovik Mountains in the east. The 
form and position of these mountains were 
found to be exceedingly favourable for the 
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development of lee-wave clouds- whenever 
WNW winds occurred. A number of lower 
mountains are found oriented in various 
directions between this range and a similar one 
partially shown on the northern part of the 
map. The Ovik Mountains, where the most 
beautiful trains of wave-clouds were observed, 
have a concave form open to the NW. The 
range is not entirely level but has a few peaks. 
The mean slope on the NW side is approxi- 
mately 1 13 and on the SE side 18... Che 
country is fairly level both on the NW side 
and SE side of the mountain and the height 
difference is about 700 metres. 


3. Occurrence of lee-wave clouds 


a. Form. The lee-wave clouds assume their 
most definite shape when the humidity is low, 
30-60 %. They then resemble a lens, an 
almond or a cigar. Seen from the side the form 
is often even more like a crossection of an 
airfoil with the blunt side turning against the 
wind. The base is normally flat or concave 
whilst the upper side is always convex and 
sharply lineated. Biconvex forms have also 
been observed. The length of the clouds are 
normally the same as the obstacle, usually 
10-30 km. The width is normally 2-10 km 
and the thickness 100 — 600 metres. With the aid 
of a good pair of binoculars one can see how 
the cloud-particles swiftly move from the wind- 
ward to the leeward side of the cloud. At times 
the cloud disappears completely, only to be 
regenerated a few minutes later at the same 
place. In spite of strong winds it has not 
changed its position just above or more often 
to the leeward side of the mountain range. 
They have most commonly been observed at 
altitudes between 4,000 and 7,000 metres, when 
they appear one by one or as a few clouds 
behind each other. At times they are magni- 
ficently coloured in red, blue, green, and 
yellow. A fibrous cloud veil (snow) can some- 
times be observed spreading from a stationary 
wave cloud far out in the direction of the wind 
before evaporation. Photo No. 1 provides an 
example. This photograph was taken on the 
Blahammaren towards Sylarna at 1345 GMT 
on July 26, 1953. Originally the cloud was 
distinctly lens-shaped with an attractive irri- 
descence. After 20 minutes a fibrous cloudveil 
spread out from the cloud. This continued for 
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Photo 1. Standing wave cloud with a veil of snow spread- 
ing out downvind. 
Photograph taken from Blahammaren towards Sylarna at 
1345 GMT, July 26, 1953. 


an hour; then the cloud disappeared entirely. 
The whole process was filmed with a time 
lapse camera which made it possible to establish 
that the wave-cloud remained stationary. The 
weather conditions were such that the pibal 
and the radiosounding from Frösön at 1500 
GMT can be assumed to be representative for 
the entire observation area. Above the friction 
layer the wind was SSW, 10 m/sec through- 
out the troposphere up to 9,000 metres 


onu 


Photo 2. À train of wave clouds formed during prevail- 
ing, WNW windt on the lee side of the Anaris and Ovik 
Mountains. The clouds remained stationary and occured 
beetween 0730 and 1030 GMT. To the left there can be 
seen a vertical pile of wave clouds with its top at 4,200 
m. The horizontal wavelength was measured to 5 km 
with the aid of an airplane. 
Photograph taken at Frösön towards the Ovik Mountains 
in SW at o815 GMT, April 17, 1953. 
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Photo 3. Wave clouds with cumulus below. 
Photograph taken ot Frösön looking west over Storsjön 
at 1215 GMT, May 25, 1953. 


and the relative humidity was 40-50 %. The 
altitude of the cloud was estimated to be 5,000 
metres corresponding to a temperature of 
puto to — 15°C. A’ stable layer which 
prevented cu-activity existed between 3,000 
and 4,000 metres. Above this layer a saturated 
indifferent stratification was found. 

When wave-clouds form “trains”, the wave- 
movement can become clearly visible as in 
photos Nos 2 and 3. 

The most spectacular formations of lee-wave 
clouds occur when wave-clouds such as those 
described above are located above each other 
in enormous piles which can extend throughout 
the whole troposphere (compare the situation 
on 5 March and 24 March). Intervals with 
clear sky occur between the individual wave- 
clouds. Piles of great vertical extension have 
been observed only in cases of strong winds 
aloft (at least 15 m/sec above the friction 
layer). A number of piles may occur after each 
other forming “trains”. These have only been 
observed on the leeward side of the mountains. 
Most of the piles were vertical but sometimes 
they tilted towards or from the obstacle. 

The author has observed stationary vertical 
piles of wave clouds extending from 2,000 m 
to 10,000 m consisting of up to 25 well 
defined cloud sheets. This happened on March 
23, 1953. On that occasion they occurred in 
a west-north-westerly current. It is to be re- 
gretted that no photographs were taken at that 
time. Sometimes these clouds are identified as 
being cumulonimbus clouds. It is easily under- 
stood that this may happen, because of the 
vertical extent of the piles. Furthermore their 
stratified structure can not always be observed 
from the ground. However, the differences 
are very pronounced. The pile of wave clouds 
remains stationary and precipitation has never 
Tellus VI (1954), 2 


127 


been observed. No significant turbulence has 
been observed in the clouds. Instead the up- 
ward and downward motion is very gradual 
but may reach values of 7-8 m/s. The meas- 
urements have been made with the aid of slow 
airplanes but never above 5,000 m. These 
values do probably not represent the maximum 
velocities occurring in such clouds. In a few 
cases measurements have been made in a single 
wave cloud or pile of wave clouds. It is then 
often found that the upward motion is con- 
siderably - sometimes 100 % — greater than the 
downward motion on the lee side of the cloud. 
This may indicate a strongly damped wave 
motion. 

At the Ovik Mountains from where the 
majority and the most well-defined “trains” of 
lee-wave clouds were observed the distances 
(the wavelength) between the piles (wave 
clouds) were in many cases measured. Because 
of simplicity these measurements were made 
at lower levels. Wavelengths between 5 and 
25 km, usually 8- 10 km, were obtained. The 
wavelength seems mainly to depend on the 
wind speed, with which it increases. In some 
cases (e.g. 5 and 21 March) the distance from 
the mountain ridge to the first wave was 
shorter than the distances between the indi- 
vidual piles. 

The structure of the wave clouds has not 
been investigated in detail. They seem, how- 
ever, usually to consist of water droplets even at 
low temperatures. Thus some are occationally 
formed on the, planeratı - 10% C or 15°C. 
In a few cases coronas have been observed 
in clouds at temperatures of - 25° C. 

b. Frequency. Wave clouds were observed 
on 125 days during a period of three years, 
October 1950 - September 1953. Distribution 
per month is shown in fig. 2. In this connection 
only clouds, the nature of which the writer 
was able to classify as stationary lee-wave 
clouds with the equipment available, are 
included. Thus lenticular clouds with no 
obvious relation to the topography or which 
could not be established to be stationary have 
not been included. 

c. The occurrence of lee-wave clouds in different 
weather situations. As will be seen in section d. 
lee-wave clouds were mainly observed during 
SW and NW winds, although some deviations 
occurred. In cases with NE-ESE winds no lee- 
wave clouds were seen over the area of 
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Fig. 2. The monthly distribution of 125 days whith stat- 
ionary lee wave clouds, observed over the area shown in 
fig. 1 during the time October 1950—September 1953. 
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observation, but such clouds then appeared on 
the Norwegian side of the range. It will also 
be seen from section d. that the lee-wave 
clouds occur in currents with practically no 
change in wind direction with height. 

When lee-wave clouds occurred in currents, 
undisturbed by fronts, mainly from NW, N 
or SW, they were usually found at high 
altitudes 4,000 -7,000 metres and nearly 
always in one or a few different layers. Humidi- 
ty was normally low 30-60 %, except possibly 
at lower levels, where condensation often 
occurred on the windward side of the moun- 
tain while no clouds at all formed on the 
leeward side. The winds were normally strong, 
never below 10 m/sec above the friction layer. 

Piles of wave clouds of great vertical extent 
on the other hand occurred practically ex- 
clusively in connection with fronts (high 
humidity) still mainly during NW and SW 
winds of great strength aloft (at least 15 - 20 
m/sec above the frictionlayer). In addition high 
reaching piles have been observed practically 
only in the peripheral parts of the fronts -in 
pre- and post-frontal positions. The whole 
month of March 1953 was marked by a 
weather situation very favourable for the 
formation of lee-wave clouds. An extensive 
anticyclone covered Central Europe, South 
Scandinavia and the British Isles, while disturb- 
ances on the polar front moved on a northerly 
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Photo 4. A pile of wave clouds formed during prevailing 
SW winds in a prefrontal cloud-system on the lee side 
of the mountains at Are. The base of the cloud pile is 
about 2,000 m and the top 4,500 m. 
Photograph taken from 2,800 m west of Storsjön looking 
northwest, at 0730 GMT, April 13, 1953. 


track in a very strong current from Iceland 
over northern Scandinavia to northwestern 
Russia. During 14 days lee-wave clouds 
appeared over the Jämtland mountains. The 
majority of situations to be described here 
occurred during this month. 

When lee-wave clouds are formed ahead of 
an approaching front they often give the 
impression of moving with the front. This is 
not the case as was often established by the 
writer during flights of long duration. Actu- 
ally they are first formed along mountains 
nearest the front and gradually new clouds are 
formed along mountains further away in the 
direction of the wind. 

Photo No. 4 shows a pile of wave clouds in 
a prefrontal position during SW wind on 13 
April 1953. According to upper level charts of 
0300 GMT a practically constant SW wind 
prevailed over the observation area. The wind 
speed was 15 m/sec immediately above the 
friction layer and increased to 35 m/sec at 
9,000 metres. The following observations were 
made at weather-flights by the author over 
the central part over the observation area 
0700-0800 GMT: Clouds in many layers up 
to approximately 2,000 metres. A new layer 
of clouds, altostratus, at about 3,000 metres and 
still another at 5,500 metres. These cloud 
layers were continuous except in lee of the 
mountains where “föhn” gaps occurred. Two 
piles of wave clouds with an approximate 
position and extent shown in fig. I, were 
observed to the NE of Are. The base was at 
about 2,000 metres, the top at approximately 
4,500 metres. Photo No. 4 shows the first pile 
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Fig. 3. Wind direction frequency for every 10 degrees at different hights. The wind data are taken from pibals made 
at Frösön when lee wave clouds were observed. The maximum of lee wave cloud occurence is found for WNW 
winds and a secondary maximum for SW winds. 


as seen from the SE. Many such piles occurred 
on the northeastern side of the mountains 
further to the south. 

The piles had disappeared at 1000 GMT 
according to reports from the pilots. The 
cloud layers were now without any holes and 
between 5,000 and 10,000 metres there was a 
continuous cloudcover. The radiosoundings 
0300 GMT from Frösön showed strong sta- 
bility in the layer close to the surface. The 
stratification was moist adiabatic between 1,000 
and 5,500 metres. A front inversion existed 
between 5,500 and 6,000 metres, with approxi- 
mately moist adiabatic stratification above. The 
humidity was high, 80-05 %. 

d. Wind conditions. The gradient wind was 
calculated for the main directions S, SW, NW 
etc. in the 125 cases when lee-wave clouds were 
observed. In proportion to the frequency thus 
obtained, 50 cases were selected in which 
upper-air wind measurements were made up 
to 9,000 metres. (Upper-air wind measure- 
ments were not carried out in all 125 cases) A 
few supplementary wind data have been 
estimated from upper level charts. Wind 
direction frequency for every 10 degree cal- 
culated for the ground, 1,000, 3,000, 5,000, 
7,000, and 9,000 metres are shown in fig. 3. 

Maximum of lee-wave clouds occurrence is 
found for WNW-winds (majority of cases 
observed at the Ovik Mountains). A secondary 
maximum for SW winds (majority of cases 
observed on the NE side of the following 
mountain ranges: the Ovik Mountains — the 
Anaris Mountains - the Helags Mountains — 
the mountains at Äre). The individual cases 
show that the wind is from approximately 
Tellus VI (1954), 2 
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the same direction at all altitudes when lee- 
wave clouds occur, at least up to the highest 
clouds (compare s and 21 March and 24 
March). This is also suggested by fig. 3 where 
mainly the surface wind, due to friction, has 
directions deviating from those at other altitu- 
des. The wind direction seems to be approxi- 
mately at right angles to the mountains gene- 
rating the wave clouds. 

The average windspeed estimated from these 
50 cases are indicated on fig. 4. In most cases 
the windspeed increased with altitude, but 
there are many cases in which the windspeed is 
constant except in the friction layer. In no 
cases lee-wave clouds have been observed 
when the wind velocity has decreased with 
height. 

e. Stability conditions. Mean temperature 
curves (fig. 4) have been calculated from 
soundings made from Frösön during March 
1953 in order to obtain an idea of the average 
stability conditions when lee-wave clouds 
occur. In these cases the soundings (altogether 
19) made when wave clouds have occurred 
(often in all layers in the troposphere) have 
been compared with those (43) when wave 
clouds have not occurred. A considerably 
higher degree of stability is found in layers 
1,000 — 3,000 metres and a somewhat smaller 
stability over 4,000 metres, than was the case 
when no wave clouds occurred. The values of 
the temperature gradient is seen in table 1. 

It should be noted, however, that these 
conditions do not entirely represent “the 
undisturbed airflow’. In the cases of NW 
winds, which were dominating, the track of 
the radiosonde balloon is outside the area of 
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Fig. 4. Mean sounding curves for the radiosonde station 

at Frösön for March 1953. The soundings (altogether 19) 

made when wave clouds have occured, have been com- 

pared to those (43) when no wave clouds occured. The 

diagram also shows the mean wind speed for different 

levels, estimated from so cases with wave clouds (those 
in fig. 3). 


the Ovik Mountains, but minor influences 
from mountains in NW can not be excluded. 
When the wind comes from SW the radiosonde 
balloon is even more affected by the disturbed 
airstream. (Compare fig. 1.) 

A study of the stability conditions in the 
individual cases (more than 125 available, often 
with wave clouds through the entire tropo- 
sphere) can be summarized as follows: 

1. À layer of great stability from the ground 
to 2,000 — 4,000 metres and aloft low stability 
(nearly moist adiabatic stratification). 
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2. When an adiabatic layer occurs at the 
ground there is always an extremely stable 
layer above it (possibly inversion) and low 
stability aloft. 

3. At times a strong inversion (front in- 
version) occurs aloft but it is more common 
that many small inversions or strongly stable 
layers are present between which the stratifi- 
cation is conditionally unstable or even abso- 
lutely unstable. 

As was pointed out by Ludlam in a note to 
Austin (1952) the undisturbed airflow in a 
wavecrest can be 1—10° colder than in the 
trough at the same altitude. A radiosonde 
passing through a lee-wave area therefore 
registrates fictitious inversions with great or 
even superadiabatic temperature gradient be- 
tween these. It has not been possible to make 
sure whether this is the reason for the stability 
conditions described in paragraph 3 above. 


4. Detailed description of three cases. 
Mars 5, 1953 


The general weather situation was marked 
by an extensive anticyclone over the British 
Isles and Central Europe, and a depression 
between Jan Mayen and Spitzbergen moving 
ESE. Its occluded front passed the Scandinavian 
mountain range early in the morning and had 
reached western Finland at 1200 GMT, fol- 
lowed by a highpressure ridge west of Norway 
bringing NW-winds over central and northern 
Scandinavia. Another low at Iceland moved 
towards the NE with its warmfront extending 
from southern Norway to Iceland at 1500 
GMT. The airflow in the middle and high 
troposphere was determined by the ridge from 
the British Isles over Scandinavia towards Spitz- 
bergen. It remained practically stationary re- 
sulting in a NW high altitude flow over the 
observation area during the entire day. 

The author made two weather flights over 
the area shown on fig. 1 between 0700 and 
1000 GMT. At first only sc clouds were 
observed over the mountains and the inland 


Table 1. Mean temperature gradient (0° C/100 m) within various layers when wave clouds occur or do not 


Waryeras (kina) pce een Oat al ee | NN So los 
Lee-wave clouds ........ 0.96 0.27 0.30 | 0.53 | 0.70 | 0.71 | 0.79 | 0.80 | 0.82 | 0.84 
No lee-wave clouds...... 0.58 0.57 0.50 | 0.61 0.64 | 0.68 | 0.73 |0.67 | 0.49 | 0.43 
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Fig. 5. Cross-section of a train of wave clouds formed on the lee side of the Ovik Mountains, March s, 1953. The 
cross-section is parallel to the average wind flow aloft. Wind data computed from upper level charts are shown to 
the right. 


but very little lee ward of the mountains. Cs 
were disappearing at the southeastern horizon. 
Fine wave clouds were observed over the 
Norwegian mountains as early as 0700 GMT 
at levels estimated between 3,000 and 7,000 
metres. New clouds were gradually formed 
near the mountains further east and at 0930 
GMT the first clouds were formed leewards 
of the Ovik Mountains (fig. ı) at many levels 
between 3,000 and 7,000 metres. At the same 
time cs appeared from the west. The wave 
clouds previously observed in the west had 
dissolved by this time -it should be noted that 
they had not been moved by the wind. Be- 
tween 1000 and 1430 GMT a series of practical 
vertical piles were formed leeward of the Ovik 
Mountains (and the Anaris mountains although 
not measured and having a less regular structure 
than those at the Ovik Mountains). Each one 
of these piles consisted of lens-shaped clouds 
lying above each other and well separated from 
each other. The first pile tilted slightly away 
from the mountain. The system remained 
unchanged but within each pile there was a 
constant forming and dissolving of the indi- 
vidual clouds. The orientation of the system, 
the dimensions and position of the piles were 
measured from aircraft at 1000 and 1400 GMT 
and with the aid of a theodolite from Frösön. 
The distance between cach pile was measured 
Tellus VI (1954), 2 


at the level where the lowest wave clouds 
occurred and was approximately 9 km. The 
individual tops of the Ovik Mountains gave 
rise to small irregularities both in the horizontal - 
and the vertical direction in the system of piles, 
but as a whole the appearence was as shown 
in fig. 5. The northernmost top of the Ovik 
Mountains gave rise to wave clouds which 
were situated on the side of the main piles and 
had a fairly irregular character. Fig. 5 shows a 
schematic cross-section of the main system, as 
observed from the aircraft and from Frösön. 
The dimensions and heights of the individual 
clouds in the first wave were precisely measured 
from an aircraft at 1400 GMT. They are drawn 
to scale and shown in fig. 5. Other piles were not 
measured but the heights were estimated from 
the aircraft and from the ground. 

Now and then a large wave cloud reaching 
above the piles I-III were formed at the 6,000 
metres level as now and then wave clouds were 
formed at about the same level above the ridge 
and at the cirrus level (hatchered in fig. 5). The 
upper side of the lowest wave clouds was 
distinctly lens-shaped while the base was diffuse. 
The upper side coincided with the lower limit 
of the temperature inversion. 12 waves were 
observed at this altitude. At higher levels a 
smaller number of wave clouds were observed. 
In addition to the wave clouds sc and fs, 5 — 8/10 
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Photo 5. A part of a train of lee wave clouds formed on the lee side of the Ovik Mountains. The cloud pile to the 
right is wave I, drawn in fig. 5. 
Photograph taken from Frösön towards the southwest, at 1249 GMT, March 5, 1953. 


occurred at a low level, very little leeward of 
the mountain where the wave clouds had 
formed. Sc clouds occurred at low altitude on 
the windward side of the Ovik Mountains, 
forming a “föhn” wall on the crest. The sc 
clouds made measurement of the upper winds 
impossible at Frösön. The upper winds on fig. 
5 are taken from upper air charts (interpolated 
for 0300 - 1500 GMT). Although the track of 
the radiosonde balloon passed 30 km north of 
the wave cloud area at the Ovik Mountains it 
is possible that the temperature sounding taken 
at Frösön at 1500 GMT is not entirely repre- 
sentative for the undisturbed current. The 
stratification is markedly stable (isothermal) 
between 2,000 and 3,000 metres, and above that 
moist adiabatic. Relative humidity above the 
inversion is 55 % -70 %. Due to a faulty 
camera only one photo was acceptable. Photo 
No. 5 was taken from Frösön towards the waves 
I-III at 1240 GMT. Wave I is best developed, 
but in view of what has been stated above it 
is difficult to identify the individual clouds 
with those drawn in fig. 5. At 1430 GMT the 
system was dissolving and further observations 


were impossible because of the appearance of 
SCHWEINE 


March 21, 1953. 


The weather situation is given by the 500 
mb map (fig. 6) and the surface map (fig. 7) 
at 1500 GMT. The strong WNW flow over 
Central and Northern Scandinavia continued 
all day. At 0600 GMT the frontal system 
located over Finland in the afternoon was 
situated off the Norwegian coast with a pre- 
cipitation area over the central and northern 
parts of the mountain range. During its passage 
over the mountains the front was considerably 
weakened by “föhn”. In the mountain area 
orographic precipitation fell throughout the 
day and up to 18 mm was measured, while no 
precipitation at all was observed on the lee- 
ward side. 

No observations from aircraft were made 
this day. At Frösön, scand as had been reported 
since midnight and nearly 10/10 of these clouds 
remained at 0600 GMT. After 0700 GMT the 
layer of sc was partly broken. There were no 
stratified clouds, sc, as or cs immediately on the 
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Fig. 6. soo mb chart, March 21, 1953, 1500 GMT. Fig. 7. Surface map, March 21, 1953, 1500 GMT. 
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Fig. 8. Radiosonde observations from Frösön for March 21, 1953 plotted on a Väisälä Emagram. The left diagram 
shovs wind data from Frösön at 1500 GMT, March 21. 1953. The full thick lines represent the observed wind, and 
the broken thick lines represent wind data computed from upper level charts. 
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Photo 6. A panoramic view looking towards the southwest from Frösön. The Ovik Mountains are situated just 
outside the right edge of the picture. Three piles of wave clouds can be seen. Their positions over the ground are 
indicated by C, D and E in fig. 1. The top of the middle pile, marked with an »x» and tilted towards the obstacle, 
was about 6,000 m. The outer piles were at times as well developed as the middle one. 
Photograph taken at 1645 GMT, March 21, 1953. 


eeward side of the mountains. Instead very 
high reaching piles of wave clouds had formed 
here. Over areas not in the immediate vicinity 
of mountains such as the Storsjö area the 
stratified clouds remained, even if somewhat 
broken up. The stratified clouds moved away 
after 0900 GMT. At that time the wind turned 
from SW to WNW and the temperature in- 
creased. Thus the front probably passed about 
0900 GMT. Afterwards the strength of the 
surface wind was about 15 m/sec. and very 
gusty. The high piles of wave clouds remained, 
however, and completely dominated the sky 
on the leeward side of the mountains. They 
were particularly beautifully developed SE of 
the Ovik Mountains where they reached the 
ci level, probably above 10,000 metres. 

Not until 1600 GMT was the writer able 
to take photographs and make measurements. 
A theodolite was used at Frösön to measure 
the bearing of three stationary piles of wave 
clouds on the leeward side of the Ovik Moun- 
tains (fig. 1). According to pibals from Frösön 
at 1500 GMT the wind direction was about 
the same, 300°, at all levels up to 5,100 metres. 
The winds at higher altitudes in fig. 8 have 
been estimated from the upper air charts at 
1500 GMT. It was assumed that the wave 
clouds were lying directly in the direction of the 
wind. Thus it was possible to determine the 
position of the piles through the intersection 
of the bearings and the wind direction line, 
the sc positions being indicated as C, D and E 
(see fig. 1). The distance between the piles was 


in this way determined to 22-23 km. The 
distance from the mountain ridge to the first 
wave was only 14 km. On some occasions such 
small distances from the first wave have also 
been measured from aircraft (compare March 5, 
1953). The top of the middle pile, marked 
with an “X” on photo No. 6, was estimated to 
reach an altitude of 6,000 metres at 1645 GMT. 
The wave clouds had reached the ci level 
before 1600 GMT and the outer piles were 
occasionally as well developed as the middle 
one. At this time new cs was being formed to 
the W. The base of the lower wave clouds of 
“rotor type” was calculated by means of a 
theodolite and was found to be approximately 
1,000 metres. 7/10 sc was found at low altitude 
over land, particularly N and E of Frésén, but 
not leewards of the Ovik Mountains. Sc 
appeared also over Frösön after 1730 GMT 
and made further observations impossible. 
When studying the temperature curve at 
1500 GMT a stratification is found which is 
typical when wave clouds appear: a “föhn” 
gradient (approximately 1°/100 metres) in the 
lowest layer, great stability (isothermal) up to 
more than 2,000 metres. Above 3,000 metres 
the stratification is moist adiabatically unstable 
up to the inversion at the 6,000 metres level 
and also above this altitude as far up as to the 
tropopause at 13,200 metres. There were, 
however, a number of minor inversions 
between 9,000 metres and the tropopause 
(the whole temperature curve is not shown in 
fig. 8). 
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Fig. 9. The right diagram shows radiosonde observations from Frösön for March 24, 1953 plotted on a Väisälä 

Emagram. The middle diagram shows wind data from Frösön at 1500 GMT, March 24, 1953. Full thick lines 

represent observed wind and broken thick lines represent wind data computed from upper level charts. The left 
diagram is a schematic cross-section of wave clouds formations on the lee side of the Ovik Mountains. 


March 24, 1953. 

In addition to piles of wave clouds this 
situation permitted some other interesting 
observations in connection with a NW “föhn”. 

According to the soo and 700 mb maps 
for 0300 and 1500 GMT the high pressure 
system remained over Central Europe and the 
British Isles (cf. 21 March) with a ridge over 
the sea between Iceland and Scandinavia. A 
very strong WNW wind prevailed over 
Scandinavia the whole day. The surface maps 
for March 24 showed a series of disturbances 
on the polar front from Iceland over Central 
Scandinavia to Western Russia. A small de- 
pression between Iceland and Central Norway 
at 0300 GMT moved swiftly ESE, immediately 
north of the observation area and had its centre 
over the Gulf of Bothnia (Kvarken) at 1500 
GMT. Its warmfront passed over Frösön at 
about 0730 GMT. Airplane reports at 0700 
GMT indicated that the cloud system con- 
sisted of stratified, fairly continuous clouds, the 
top of which were above 5,000 metres. After 
the front passage a mild strong and gusty 
WNW airflow penetrated into the observation 
area. It lasted all day and increased in strength. 
Tellus VI (1954), 2 


Following the frontal passage some clouds 
disappeared and the cloud cover was fairly 
constant over the observation arca after 0800 
GMT and for the rest of the day. At 10,000 
metres there was 4 — 8/ro cs and ci. Turbulence 
strato-cumulus, $ — 8/10, appeared at low level 
over land, particularly north and cast of Frösön, 
while there were no low clouds above Storsjön 
or immediately on the leeward side of the 
Ovik Mountains. A “wall of clouds” remained 
stationary all day west of a line approximately 
through the Anaris Mountains. Its altitude was 
estimated to approximately 5,000 metres and 
it was very similar to the cloud-system of a 
cold front. NW of this line large quantities 
of precipitation fell during the day (20-45 
millimetres). There was no precipitation in the 
southern and eastern parts of the observation 
area, apart from a few places, where precipita- 
tion totalled 1 millimetre. 

After the clearing the writer observed a small 
pile of wave clouds at 0800-0930 GMT. It 
remained stationary leewards of the Ovik 
Mountains and was later dissolved. Its position 
and altitude was determined by means of a 
theodolite from Frösön (see figs. ı and 9). The 
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Photo 7. Wave clouds at one level and with a short wave 
lenghth (about ı km). 
Photograph taken from Frösön towards the southwest, at 
1200 GMT, March 24, 1953. 


gradient wind at the time was WNW and 
approximately ıs m/sec. (Upper winds were 
not measured. 

Later a pilot reported (competent glider- 
pilot and meteorological observer) that he had 
determined the position and measured a single 
mighty pile of wave clouds at 1030 GMT. It 
was situated at F, see figs 1 and 9, that is about 
18 km from the mountain crest in the wind 
direction. In addition to this pile and cs, and 
at low height sc, there were no other clouds 
in the vicinity. He described it as reproduced in 
the cross section, in fig. 1. The pile consisted 
of 13 wave clouds well separated from each 
other, with the top layer at 7,500 metres. The 
pile was vertical, and its horizontal extension 
at right angles to the wind was equal to the 
width of the Ovik Mountains. The pilot had 
been unable to observe any movement of the 
pile during the 30 minutes he observed it. In 
the afternoon the same pilot reported that the 
pile had “moved” and that it was situated at 
G at 1400 GMT, that is 48 km from the moun- 
tain crest (see figs 1 and 9). The highest wave 
cloud was now at 7,300 metres. The pilot 
described the pile as “being not quite as regular 
in shape as previously” 

Unfortunately the writer was unable to make 
any observations from the air himself, and the 
piles were not visible at Frösön for sc clouds. 
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Photo 8. Clouds rotating around a horizontal axis and 
simultaneously mowing downwind. They appeared in a 
layer of strong wind discontinuity and turbulence (see 
fig. 9) at an altitude of about 6,000 m. 
Photograph taken from Frösön looking towards the east, 
at 1255 GMT, March 24, 1953. 


Some other interesting observations were 
made, however. The wave clouds appeared at 
one level (see photo No. 7) and with a short 
wave length (about 1 km) leeward of the Ovik 
Mountains for a period of half an hour. The 
photo was taken from Frösön towards the SE 
at 1200 GMT. The altitude of the clouds was 
estimated at 6,000-7,000 metres. They ap- 
peared to be parallel to the Ovik Mountains. 
Photo No. 8 shows another type of clouds 
which appeared chiefly over the Storsjö area, 
but not leewards of the Ovik Mountains in the 
afternoon. The photo was taken from Frösön 
towards the east at 1255 GMT. The clouds 
were in a strong rotating motion around a 
horizontal axis, at the same time as they moved 
ESE. According to reports from an aircraft 
they appeared at an altitude of 6,000 metres. 
They were thus situated in a layer of strong 
wind discontinuity (see fig. 9). The pilots found 
a layer of strong turbulence between approxi- 
mately 5,500 and 6,500 metres with accelero- 
meter readings of 6 g. There was likewise 
severe turbulence in the surface layer below the 
sc clouds with accelerometer readings of 7 g 
at speeds of 450 km/hour. As a result of their 
construction accelerometers do not give true 
readings under conditions of turbulence, but it 
was the general opinion of the pilots that the 
flying conditions were most difficult both at 
6,000 metres and in the turbulent layer near 
the ground. Turbulence was also noticed at the 
tropopause level, which was described as light- 
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moderate within the 9,500 - 10,200 metres 
layer. Condensation trails from aircraft oc- 
curred within the layer at altitudes between 
9,200 and 10,200 metres. The trails were inter- 
mittent even when the airplane maintained 
steady altitude and airspeed. This seems to 
indicate differences in horizontal humidity 
distribution even within very short distances. 


A study of the wind distribution at 1500 
GMT shows that, except in the friction layer 
where the wind was very gusty, the wind 
direction was approximately 305° throughout 
the troposphere. The piles of wave clouds were 
oriented exactly in this direction (fig. 1). 
Disregarding the friction layer the wind 
strength was practically constant up to 5,000 
metres. À strong wind discontinuity existed 
between 5,000 and 6,000 metres. The wind 
velocity above 7,000 metres was estimated from 
the 300 and 200 mb maps at 1500 GMT. The 
upper air maps showed no noticable change in 
the direction and strength of the current at 
higher levels between 0300 and 1500 GMT. As 
already pointed out the wind intensity in- 
creased considerably at lower levels particu- 
larly in the friction layer afternoon. There is 
no marked difference between the 0300 and 
1500 GMT stratification curves. A minor 
radiation inversion at 0300 GMT was dissolved 
during the morning. A practically adiabatic 
temperature gradient was found near the 
ground and great stability in the 1,000 — 3,000 
metres layer. Above, the stratification was 
practically moist adiabatic up to 7,500 metres. 
Humidity decreased from 0300 GMT and was 
between 80 and 70 % above the stable layer 
at 1500 GMT. 

Observations from the air were discontinued 
after 1500 GMT. There was no change in 
cloud conditions before nightfall. 


5. Comparison between observed and theo- 
retical results. 


There are many theoretical works dealing 
with the deformation of an air current passing 
over a mountain ridge, and the formation of 
vertical, stationary wavesystems on the lee- 
ward side of the ridge. It is difficult to apply 
these theories to a real situation, however, due 
to the various simplified assumptions which 
must be made. 

It may be of interest, however, to investigate 
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whether the formula for the oscillation of an 
air parcel in a stable stratified atmosphere can 
be applied. 


The time of oscillation + is expressed 


250 I (1) 
Vera 


where Ty is the initial temperature, yg the dry 
adiabatic temperature gradient and y the pre- 
vailing temperature gradient. 

Here, however, we are not considering single 
air parcels in vertical oscillation, but an entire 
air layer, which is in vertical oscillation on the 
leeward side of the mountain ridge, which can 
be approximately regarded as sinusoidal oscil- 
lations. Hence, the average values, applying 
to the entire oscillating layer, for temperature 


(T) and the prevailing temperature gradient 
(y) should be introduced in formula (1). When 
we determine the wave length L we should 
also use the average value for the wind speed 
Us 

\ Applying this to the layer where wave 
clouds appeared in the cases described in this 
article, on s, 21 and 24 March we shall find: 


March 5, 1953 
T=245°K; y=0.0061°/m and U=18 m/s 


that is t=516 sec. and the corresponding 
wave length L=9,288 metres (9.3 km). 

The wave length actually measured was 
approximately 9 km. 


March 21, 1953 
T=265° K; y=0.0068°/m and U=30 m/s 
that is T— 596 sec. and L=17,880 metres (18 


km) 5 


The measured wave length was 22-23 km. 
March 24, 1953 
T=252°K; y=0.0060°/m and U=45 m/s 


that is t=517 sec. and L=23,265 metres (23 
km). 

In this case there was only one pile of wave 
clouds and it is doubtful whether its distance 
from the mountain crest (48 km at 1400 GMT) 
represents the actual wave length. 

In summing up it can thus be said that good 
correlation between the optically measured and 
the theoretical wave length was obtained for 
the case on March 5, 1953 and to some extent 
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also for that on March 21, 1953. In some other 
cases the wave length has been measured and 
also computed from formula (1). Fairly good 
agreement has been obtained. 
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Abstract 


A series of twenty-four 24-hour forecasts have been completed using the barotropic model 
of the atmosphere. The results of these forecasts are better than previous tests with this model 
have indicated. The average correlation coefficient between computed and observed height 
changes of the soo mb surface is 0.77, the error in the forecast height changes is 89 m against 
an average observed change of 121 m. The errors are mainly due to the following factors: 


1) Erroneous assumptions at the boundary of the forecast area. 
2) Inaccurate analysis of the initial data over the Atlantic. 
3) Inaccuracy in forecasting components of the motion of a scale only a few times larger than 


the grid size. 
4) Baroclinic developments. 


The computations indicate that the present results may be improved considerably by trying 


to eliminate errors listed under (I) - (3). 


I. Introduction 


During the last years a number of models of 
the atmosphere have been developed for numer- 
ical forecasting purposes, making use of infor- 
mation from one, two or several levels in the at- 
mosphere. The computational schemes for these 
different models involve very large computing 
programs in order to obtain a 24-hour forecast 
for a significantly large region. The intention 
has been to test and evaluate the proposed 
methods by actual computations with the aid 
of electronic computing machines, but it is 
not until recently, that such machines have been 
available on a wider basis. Hence very few 
systematic and conclusive evaluations based 
on actual computations have been made up to 
the present time. The first more extensive 
study was published recently by CHARNEY and 
Puirrirs (1953). Altogether six barotropic and 
six two-layer forecasts for 24 hours were 
presented. All these forecasts were made 
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using data from a case of extreme baroclinic 
development and the two-layer model gave 
considerably better results than the barotropic 
one. Since then a few additiona forecasts 
hava been made at Princeton but no additional 
results have been published as yet. Therefore 
we do not yet know the value and limitations 
even of the barotropic model. It seems im- 
portant to get a good idea of the potentiality 
of this simple model, since more complicated 
models should be tested with regard to the 
improvement they may represent in compari- 
son with the barotropic one. 


Partly for the reasons mentioned above, 
partly because of the limited capacity of the 
machine in Stockholm at the beginning of the 
present computational program, it was decided 
to investigate the barotropic model more ex- 
tensively. Furthermore, it is reasonable to start 
with the simplest model to obtain experience 
before using more complicated ones. As will 
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be seen from the results presented here of alto- 
gether twenty-four 24-hour forecasts a number 
of basic problems associated with numerical 
forecasting seem to require further considera- 
tion and are more easily studied with the aid 
of the barotropic model than with the more 
complicated ones. 


The results presented in this report represent 
a continuation of work started during the 
autumn of 1951 by a group of meteorologists 
assembled at the University of Stockholm. The 
work of the group has since then been supported 
by grants from the Knut and Alice Wallenberg 
Foundation, the U.S. Office of Naval Research, 
the Swedish National Science Research Council, 
and the U.S. Weather Bureau. Very generously 
the Swedish electronic computer (BESK) has 
been placed at our disposal by the Computer 
Board. The participants in the active work 
have been: G. Arnason, Iceland; H. Bedient, 
U.S.A.; P. Bergthorson, Iceland; B. Bolin, 
Sweden; G. Dahlquist, Sweden, B. Döös, 
Sweden; N. Phillips, U.S.A. Some of the 
maps used in the series tested had previously 
been analysed by W. H. Hubert, U.S.A., Ch. 
Newton, U.S.A., and L. Vuorela, Finland. 


2. Method and organization of the 
computations 


By and large the computations have been 
done in the way described by CHARNEY and 
PHirires (1953) and we shall here restrict our- 
selves to a very brief account of the procedure. 


The integration is done by an iterative 
process in which each time step T to T+I may 
be described by the following four equations 


nis the absolute vorticity, m the magnification 
factor of the map projection used, fthe Coriolis 
—9 


parameter, & = f. m? &, ¢ being the relative 
vorticity, 9 the geopotential and At is the 
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timestep used in the extrapolation. J (7, 9) is 
evaluated by finite differences and the Poisson 
equation (4) for @ is solved by relaxation. The 
history of the motion is thus carried by £, 
the relative vorticity, and the new values of 
p are obtained from (4) at each time step. 
Certain boundary conditions have to be fore 
mulated in order to be able to solve the system 
of the equations (1) - (4) uniquely over a lim- 
ited area. It was shown by CHARNEY, FJORTOFT 
and von NEUMANN (1951) that a knowledge 
of y at all times along all boundaries and of 
the vorticity, ¢, at points of inflow probably 
are necessary and sufficient conditions to deter- 
mine the motion uniquely inside the area. We 
do not know what is going to happen along 
these artificial boundaries and we must make 
certain assumptions regarding q and &. These 
assumptions are necessarily in error. They will 
influence the motion and introduce errors in 
the forecasts. However, it is hoped that these 
errors will not propagate far in from the bound- 
aries over a period of 24-hours. We shall later 
obtain some informations about the validity 
of this statement. It was here assumed that 
dpldt = 0, e.g. 9 = const. on the boundary. 
Furthermore & was set equal to its initial value 
#0 at points of inflow and also at points of out- 
flow. This latter assumption at points of out- 
flow really means an overspecification of the 
problem. It was done in order to simplify the 
code and save storage space in the machine. We 
shall furthermore see that the errors thus intro- 
duced usually are not larger over 24-hours than 
those introduced by the arbitrary assumption 
of & = £0 at points of inflow. Both may, how- 
ever, cause considerable errors quite far from 
the boundaries but they seem to be of quite a 
different nature. The values of &® on the boun- 
daries were obtained by supplying information 
of 9" at one row of points outside the area later 
used for the computations. In that way no ex- 
trapolation procedure had to be used. 
Equation (4) was solved by relaxation 
methods. The rapidity of convergence of this 
method depends to a large extent upon a good 
first guess. The best scheme is undoubtedly to 
extrapolate values of @"+! linearly from the 
knowledge of 9-1 and g'. This requires, 
however, the storage of four quantities at each 
grid point instead of three (EI SEE andren: 
The limited present memory of the computer 
did not permit us to apply this method. g* was 
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Fig. Ia. soo mb map at 0300 GMT 
November 25, 1951. Contour lines 
drawn for every 40 meters and labelled 
in decameters. 


used as the first approximation in the relaxation 
to obtain g**}. 

The computations were made over a square- 
grid consisting of 20 x 20 points and the grid- 
size was chosen to be 300 km at 50° N. The 
map projection was polar stereographic. Thus 
the total area covered by the grid was about 
5,700 x 5,700 km. The location of the area var- 
ied slightly for the different forecasts, but did 
never deviate much from the one shown in 
fig. 1a. The forecasts were only verified over a 
grid 11 x II or 12 x 12, in the centre of the area 
since the effects of the boundaries were sup- 
posed to effect significantly the region outside 
this area (see fig. 1b). The time step was 1 hour. 

The computer used for the computation has 
an internal William-type memory of 512 
words of 40 binary digits, which in the near 
future will be supplemented by two magnetic 
drums, each of which will be capable of hold- 
ing 4,096 words. The adding time is 56 micro- 
seconds and the multiplication time is 364 
microseconds. 

It was already mentioned that three quan- 
tities were stored for each grid point, each of 
them with an accuracy of 13 or 14 binary 
digits making it possible to store all informa- 
tions needed for one point in one word. The re- 
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maining 112 words were used for instructions. 
This was not suflicient for storing all instruc- 
tions needed and thus new instructions had to 
be read in continuously during the computa- 
tions. This fact, in addition to the poor initial 
guess when solving equation (4), increased the 
computation time by a factor of at least 4 or 5. 
With the method used here a 24-hour forecast 
required 40 minutes machine time, of which 6 
minutes were used for printing the final re- 
sults. The increase of the storage capacity of the 
machine will make it possible to decrease the 
actual computing time to about 1/5. 


In a few cases the forecasts were made on an 
operational basis, e.g. the forecast was com- 
pleted as soon as possible after the time of 
observation. This was done in order to gain 
experience for routine forecasting by numeri- 
cal methods. The preparation of data for a 
forecast involves a series of steps. The time re- 
quired for a forecast after a certain skill in the 
routine operations has been aquired may be 
estimated to the following values: 


Checking and plotting of data for a 


s00.emb, mapa (Aypeople) jane. th. 30m. 
Analysis of soo mb for the fore- 
cast area and adjacent regions 1h. 30 m. 


November 25 1951 
O3GMT 


Reading off values of y at 440 grid- 


points and checking (2 people).. 1h. 20 m. 
Punching of input data on tape 40 m. 
Checking and correcting tape ... 20 m. 
Forecasting with the machine ... 40 m. 
Plotting and analysis of the final 

FOLECASEVR ERE eee eee Re ELS 30m. 

6h. 30m. 


Some of these steps may be speeded up, others 
carried out simultaneously, but the time to 
be gained in that way is very limited. In the 
case of the barotropic forecast the time given 
above is not prohibitive if one can extend the 
forecast range to 48 or 72 hours but as soon 
as we want to extend the computations to 
two or several layers in the atmosphere other 
schemes for preparation of the data have to be 
devised. Furthermore certain errors were still 
found in the initial values in spite of the 
checking routine which was adopted here. 
Some work to overcome the practical diffi- 
culties has been attempted: 

1) A method of numerical synoptic analysis 
of the initial data seems very important from 
the point of view both of the time required 
using conventional methods of analysis and 
the probability of introducing errors by man- 
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Fig. 1b. Computed changes (dashed 
lines) and observed changes (solid lines) 
of the height of 500-mb surface from 
November 25, 1951, 0300 GMT to 
November 26, 1951, 0300 GMT. Lines 
are drawn for every 100 meters. 


ual methods. This problem is at present being 
studied by some members of the group. 

2) As long as methods for numerical anal- 
ysis are not available the input program must 
be constructed in the most flexible way. Thus 
it should be possible to read in the data in a 
more or less arbitrary order and to introduce 
corrections to these data afterwards and without 
difficulties. The continuity and smoothness of 
the field should be checked automatically and 
if not fulfilling certain criteria automatic cor- 
rections should be applied. One such checking 
scheme was developed by H. Bedient. It con- 
sisted in computing \72\7?@ initially, which 
field quite clearly shows unreasonable values 
of y. This field was printed with one signif- 
icant sedecimal figure and inspected visually. 
It proved quite helpful on a number of occa- 
sions. 

3) The final results were obtained as a table 
of the values of gt and £24 at all grid points. 
The machine may, however, be instructed to 
interpolate between the actual grid points and 
also print the results in the form of a map in 
the same scale as used for preparing the initial 
data. Such a scheme has been worked out by 
G. Dahlquist and H. Bedient and fig. 2 shows 
an example of such a representation of a map. 
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Fig. 2. The initial soo mb height field on 24 April 1954, 
0300 and its graphical representation by typewriter out put. 

The solid lines show the observed height contours of 
the soo mb pressure surface and the typewriter figures 
fill the space between the lines 5,000 meter—s,040 meter; 
5,080 meter— 5,120 meter; etc. The memory of the 
computer contained height valuesin a coordinate system 
x=iNS,y=jAS where ASis the grid interval on 
the map and i and j are integers. The typewriter prints 
digits located in a coordinate system x’ =i’ AS, y = 
—j AS where ij are not integers but 01 = aK, 
öi’/öj = B. « is the ratio of the original map scale and 
the new scale made by the typewriter and ß and K are 
constants depending on the construction of the type- 
writer. At each point x’ y’ the height is interpolated from 
the four corner values of the grid square and converted 
to a single sedecimal figure (base 16) such that o = 5,000 
meters and the value increases by 80 meters for each unit. 
A decision to print the figure is made if the value lies 
between the integral value and ?/, unit higher. The deci- 
sion is made to print nothing if the value lies between 1/, 
and the next integral value. Some “nervousness’’ is shown 
in the through west of Ireland because of the linear 
interpolation scheme. Also the low in the English channel 
is not well represented because of its small size in relation 

to the grid. 


3. Results 


Forecasts have been made from two series 
of maps: 24 November 03 GMT - 27 Novem- 
ber 15 GMT 1951 consisting of eight 24-hour 
forecasts from maps with 12 hour intervals; 
1—10 January 1954, altogether ten 24-hour 
forecasts from maps 24 hours apart. In addition 
to those a number of individual cases were 
treated, chosen for various particular reasons. 
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The January series included the development 
of two severe storms within the forecast area. 

The results of the forecasts are listed in tablel. 
They represent evaluations over the verifica- 
tion area in the centre of the forecast area. 
r denotes the correlation between y, the com- 
puted 24 hour changes of y, and x, the corre- 
sponding observed changes; o, and o,, are the 
root mean squares of the observed and com- 
puted changes respectively; ¢ denotes the 
errors of the computed changes; x and y 
represent the average observed and computed 
I 
N 

The statistical treatment of the results ob- 
tained here is by no means exhaustive. It has 
often been pointed out that the correlation 
between observed and computed changes is 
not a satisfactory measure of the goodness of 
a forecast. This is certainly true. It is believed, 
however, that the quantities listed in table I 
will give a fairly good idea of the results. 
Furthermore the correlation coefficient often 
has been used in evaluating conventional meth- 
ods of forecasting and a comparison with such 
evaluations is desirable (cf. LONNQUIST, 1952). 
We shall later discuss the forecasts in terms 
of predicted changes in the winds over the 
British Isles. 

The following remarks are of particular 
interest in connection with table I: 

The average correlation coefficient for the 
24 forecasts is 0.77, the error in the forecast 
height changes is 89 m in comparison with an 
average observed change of 121 m. The eight 
first forecasts (24-27 November) were made 
from very carefully analysed maps while the 
remaining 16 forecasts were based on analyses 
probably comparable with routine analyses 
in a weather service. There is a distinct differ- 
ence between the results. The mean corre- 
lation coefficient for the eight forecasts in 
November is 0.90. The weather situation 
during this period was characterized by a 
series of regular long waves moving east- 
wards in a fairly strong westerly current in 
middle latitudes (cf. fig. 1 a). No particularly 
strong baroclinic developments occurred dur- 
ing this period. However, regular cyclones 
formed, occluded and dissolved while moving 
from west to east. The difference between 
this first series and the remaining forecasts 
cannot be explained merely as a result of the 


changes (x =—2x, y= 27): 
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x 

Date of initial map momie | + | 6 | & | 6 | em | & [ee 7 | (a) | EN | (ma) | io) oylox 

24, Nova 51 03e ee 0.83 114 116 70 — 16 — A1 1.02 
24 » BSD ETS PDO 0.91 120 130 54 6 3 1.08 
25 >» HBT MOB tts bec 0.95 155 153 48 26 24 0.99 
25» SS Tl TS snes ater 0.96 173 209 73 33 — I 1.20 
200) 5 LOG Rene 0.79 148 164 99 34 37 Tr 
26 =» 51) 15 eue 0.96 272 208 105 18 98 er 
ZT BIN OFEN cr 0.95 160 185 67 32 44 1.16 
2/7) SN EI CE 0.85 160 158 IOI 57 8 0.99 
2h SEN GH, Of 0 oncide & 0.75 86 62 56 — 20 — 13 0.72 
24 » iP OR oe 0.91 140 93 qe — 75 — 29 0.66 
I Jan SAN O3 eee O7 7e 204 169 45 180 2.83 
Zor ate O DSC ele 0.80 107 107 71 45 83 1.09 
30° SGV OS 0.76 126 104 89 — 20 gar 0.83 
A yy SOS er 0.72 17e 133 180 — 82 43 0.77 
5.» BA, OB sere derer. 0.92 182 209 84 7 — 6 1.14 
On» SANDOR wei een: 0.83 102 174 118 18 — 36 1:72 
Te) HAN OS Mec er 0.38 67 148 124 30 130 2221 
8 » MO io das 0.83 98 112 63 35 68 1.15 
9 » CP ao ood ood 0.73 83 Va. 65 20 II 0.93 
Io» 254, 03cm 0.91 118 159 64 — 68 — 77 1.34 
232 Manch2-54, 032 2. 0.60 89 109 OI 6 5 1.22 
24 » = 54) BOR Wickes peers 0.40 I2 178 152 34 80 1.38 
MA) DIRES EOS ee 07 89 88 65 20 — 5 0.99 
24 » ATOS ARR 0.76 MSN LES nn 6 7e Ma TT Er SNA 83 67 I — 51 1.86 
Wversallsmeansce rer | 0.77 | TOT: | 140 vr [owe | | als 89 | 8 | 24 | 1.16 
24 27 Nov. -51......... 0.90 SC one | 165 77 — — 1.10 


1 These forecasts were done on an operational basis. 


different character of the general weather 
situation. The very best forecast (r = 0.95, 

= 48 m), was made during this November 
period: 25 November 03 GMT, which is 
shown in fig. 1. 

Tendency computations had previously been 
carried out for one of the maps in the Novem- 
ber series: November 25, 15 GMT (cf. STAFF 
MEMBERS, UNIVERSITY OF STOCKHOLM, 1952). It 
is of some interest to note that the complete 
24-hour forecast gave considerably better 
result than the estimated 12 hour changes from 
the tendency computations. 

The average computed 24-hour changes are 
15 % larger than the observed ones. The ratio 
0,/o, varies, however, considerably from case 
to case. One reason for this can be found by 
comparing x and y. It is seen that x and y 
usually are quite different in the cases where 
dy/o, deviates considerably from unity, e.g. 
an erroneous rise or fall over the whole area 
has been obtained from the computations in 
those cases. The errors (e) of the prognostic 


height changes are also largest in these cases. 
This is particularly true on Jan. 1, Jan. 4, 
Jan. 7, March 24 and April 24. These fore- 
casts are some of the least successfull. 

It is quite obvious that we cannot expect 
to be able to forecast the motion and devel- 
opment of systems whose scale is larger than 
the size of the forecast area and it is therefore 
not surprising that x-y may be quite large 
in some cases, even if there is a tendency for 
both x and y to approach zero when the size 
of the area increases. Let us introduce the 
notation 


where B is the periphery of the forecast area. 
The computed value of S is of course zero 
because of our boundary assumption, but this 
is not the case for the observed value of S. 
It may deviate from zero by 20 m or at a 
maximum 25 m, for an area the this size. The 
erroneous assumption concerning the value 
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Fig. 3 a. Observed contours (solid lines) and absolute 
vorticity (dashed lines) on November 24, 1951, I,500 
GMT. The heights are given in decameter as unit 
and the vorticity-lines are labelled in units of 1075 
secu. 


Fig. 3 b. 24-hour forecast from map shown in a. 


of g on the boundary then gives rise to an 
error within the verification area of about 
the same amount as the value of S. 

More important in explaining the difference 
between x and y seems, however, the second 
boundary condition about the vorticity &. 
It was earlier mentioned that we assumed 
&* = £ at all boundary points. This is wrong 
in two respects: 

a) It means an over-specification in that & 
at points of outflow is determined by the 
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November 24,1951 


29 


Forecast from 
November 24, 195] 


values of & inside the boundaries (if Op/dt=o 
at the boundary) and cannot be prescribed. 
The effect is very clearly shown in fig. 3. 
First the vortex lines are advected towards 
the boundary and since £*=&? on the bound- 
ary the gradients increase which means an 
amplification of higher harmonics in the flow 
field. The finite difference formula is quite 
inaccurate for these higher harmonics and 
they are distorted considerably during the 
computations. This gives rise to the for- 
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mation of a system of waves. Gradually the 
effect spreads inwards from the boundary. In 
most cases it was, however, hardly noticeable 
in the o-field inside the verification area. It 
should be remarked that a certain tendency for 
a “blocking” of this kind was present on the 
map shown in fig. 3. The pattern was, how- 
ever, considerably distorted during the com- 
putations. 

b) On the boundary segments with inflow 
the true value of the vorticity in reality is 
determined by the vorticity field outside the 
area. If such a boundary intersects a region of a 
strong gradient of the vorticity field perpen- 
dicular to the boundary the prescribed value 
of the vorticity on the boundary will differ 
considerably from the actual value within a few 
time steps. In the case-of strong flow across the 
boundary the effect of such erroneous values of 
the advection of vorticity across the boundary 
may be noticed far away from it. The total gain 
or loss of vorticity over 24 hours for the whole 
forecast area by advection across the bound- 
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Fig. 4 a. 500-mb 
map at 0300 GMT 
January I, 1954. 
Notations as in 
Here 


ary is casily evaluated. It has been done in a 
few cases and the agreement is poor between 
the observed change of vorticity and the one 
computed on the basis of the boundary condi- 
tions. The forecast from January ı is the most 
striking example of this effect (fig. 4). In this 
case the left boundary was placed to the east 
of a deep trough over western Atlantic. The 
initial values of the vorticity on the boundary 
were kept constant for the whole forecast 
period and gave rise to negative values of the 
vorticity advection over most of the western 
parts of the area. A very strong rise resulted. 
In reality, however, the trough to the west 
moved across the boundary, caused a con- 
siderable deepening and brought large values 
of absolute vorticity into the area. The fore- 
cast was an almost complete failure (cf. table I). 
It should be remarked that the analysis over 
the Atlantic never can be very accurate and 
analysis errors may to some extent have 
contributed to the poor result but there is 
nevertheless no doubt about the importance 
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Fig. 4 b. Computed and observed 24-hour changes 
from map shown in fig. a. (Notations as in fig. I.) 


of the effect described above!. Also in some of 
the other cases of unsatisfactory results the 
same effect has played an important role. 

It is quite evident both from table I and 
fig. 5 that the forecasts are better the more 
pronounced the actual changes of the flow 
pattern are. The same thing has been noticed 
with regard to standard methods of fore- 
casting (LÖNNQVIST, 1952). 

The goodness of the forecast in different 
parts of the verification area has been investi- 
gated in the same Way as in the previous report 
(STAFF MEMBERS, UNIVERSITY OF STOCKHOLM, 
1952). The result is shown in fig. 6 a in 
graphical form. We notice the relatively poor 
results in the western parts of the area, caused 
by the erroneous boundary values. Over north- 
ernmost England the! average correlation is 
above 0.90 in spite of the fact that even the 
least successful forecasts have been included. 
This is certainly a very promising result. Fig. 7 
shows the observed height changes at Liver- 
pool (solid line) and the corresponding com- 

uted changes at a grid point close to this 
station (dashed line). We cannot forecast the 
mean height change (y) with any accuracy. 
This fact influences the results considerably. 
It is conceivable that most of the errors thus 


1 Note added in proof: The map for January 1 has later 
been reanalysed and another 24-hour forecast has been 
performed. The results were slightly improved (r = 0,30) 
but the main features of the forecast changes as shown 
in fig. 4 remained unchanged. 
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introduced would be eliminated if a larger 
arca could be used for the computations. We 
can get an approximate idea of the improve- 
ment it would represent by adding (x-y) to 
the forecast values at every point. Fig. 6b 
shows the geographical distribution of the 
correlation coeflicient after applying this correc- 
tion and the dash-dotted curve in fig. 7 in- 
dicates the result of this correction at a partic- 
ular point. 

Since one of the main objectives of upper air 
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Fig. 5. Correlation coefficient (r) plotted as a function of 
observed changes (x). 
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Fig. 6. Geographical distribution of the correlation between observed and computed 24-hour changes. Fig. a 
shows the correlation between x and y at Is points, fig. b shows the correlation obtained between x and 


(y +x—y). The squares show the two most common locations of the verification area. 


forecasting is wind forecasting it is of some 
interest to express the results obtained here in 
terms of the wind. The computed and ob- 
served height changes at three points (A), 
(B) and (C) (see fig. 6a) have been used for 
evaluating the wind changes. The mean ob- 
served wind speed over this area was 26 m/s. 
The average observed change of the wind 
over 24 hours was 20 m/s against a computed 
value of 21 m/s. The mean error of the wind 


forecast was 13 m/s, which seems to be inde- 
pendent of the intensity of the wind field or 
the rapidity of its changes. All these values 
refer to the total wind. These results clearly 
show that considerable errors in the wind 
forecast exist in spite of the relatively high 
correlation (r=0.90). A verification of the com- 
puted wind changes should therefore be given 
as a complement of ordinary correlation co- 
efficients. 
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Fig. 7. Observed changes (x) of the soo mb surface over Liverpool during November 24—27, 1951 and January 


1—10, 1954 (solid lire). Computed changes (y) (dashed line) and computed and corrected changes (y + x— y) (dash- 
dotted line) at gridpoint close to Liverpool. 
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4. Conclusions 


The results of 24-hour forecasts with the 
barotropic model presented here are better 
than previous test (tendency computations and 
complete forecasts) have indicated. The errors 
seem to depend on the following factors, 
listed in the order of importance indicated by 
the present computations: 


1) Erroneous assumptions about the values of 


gp and in particular € on the boundary. This de- 
fect may most effectively be eliminated by 
extending the forecast area beyond regions 
of strong vorticity gradients, e.g. by en- 
closing the larger portion of middle and high 
latitudes of the northern hemisphere. The 
memory of the machine after incorporating the 
magnetic drums, will easily permit the storage 
of this amount of data. It will then also be 
possible to extend the forecasts to 48 and 72 
hours and investigate the break-down of the 
barotropic model (or of the computational 
scheme!). Such a program is under develop- 
ment. 

2) Inaccuracy of the initial analysis over areas 
with sparse data, e.g. the Atlantic. For the next 
series of forecasts we are planning to utilize 
the result of the previous 24-hour forecast 
as an aid in the analysis over regions with 
sparse data. This method has previously been 
suggested to be one of the basic principles of 
numerical analysis of the initial data. 

3) Inaccuracy in forecasting components of the 
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motion of a scale only a few times larger than 
the grid size. For example, the speed of a wave 
of wave length 1,200 km is about 30% in 
error. This probably necessitates an intitial 
smoothing of the vorticity field to remove 
the small components of motion, which cannot 
be computed accurately. A continuous activa- 
tion of the small scale motion takes place 
because of the non-linearity of the forecast 
equation. This tendency is noticeable in some 
of the forecasts and the wave pattern in fig. 3 
may partly be looked upon as a result of this 
effect. When extending the forecasting to 
longer periods a smoothing even during the 
computations may therefore be necessary. 


4) Baroclinic developments. The errors in- 
troduced by using the barotropic model are 
often obscured by the other sources of errors 
mentioned above. In extreme baroclinic de- 
velopments as during the period studied by 
CHARNEY and Puirurps (l.c.) the barotropic 
model is not satisfactory. However, the re- 
sults obtained by these writers do not seem to 
be representative for the applicability of the 
barotropic model in general. It is a common 
impression within our group that the results 
presented here by no means represent the ulti- 
mate limit of the potentialities of the baro- 
tropic model. On the contrary, a number of 
modifications may be carried through which 
will probably still improve these results sig- 
nificantly. 
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Abstract 


By making use of the mathematical properties of certain special types of weighted means — 
taken either areawise or timewise —it has been shown that a variable’s deviation from its mean 
can be expressed in terms of the mean itself. This result has been used to derive a closed system 
of equations for the mean motion of an incompressible autobarotropic fluid. Comparison 
of the equations for area-averaged and time-averaged motion reveals that the area mean has 
certain advantages, inherent in the form of the equation for the unaveraged motion. The eddy 
transport for the case of “strong averaging” is contrasted with that derived from the classical 
“mixing-length’”’ theory of eddy diffusion. The usual restrictions on the statistics of eddy fluctua- 
tions are found to be inconsistent with the unrestricted form of the eddy transport. Finally, 
it has been shown that these methods can be extended to the derivation of closed systems of 
equations for the mean motions of simple baroclinic fluids. 


1. Introduction and General Remarks 


In the view of most synoptic meteorologists, 
one of the marks of a good analysis of pressure 
or streamline patterns in the upper troposphere 
is its ‘smoothness’ — i. e., “smooth” in the 
sense that it does not reflect all the minor 
irregularities weakly indicated by the raw pres- 
sure or wind data. To cite an example, the 
direction of the analyst’s streamlines in the 
immediate neighborhood of an observing sta- 
tion frequently does not coincide with the 
reported direction of the wind at that point, 
but conforms to the “general” direction of the 
winds over the surrounding area, or lies some- 
where between those two directions. Con- 


sciously or otherwise, the analyst evidently 
recognizes that wind reports are at least subject 
to instrumental and round- off errors, and may 
be “unrepresentative” in the respect that some 
of the observed irregularities in wind direction 
may be associated with disturbances whose 
scale is comparable with or less than the distance 
between adjacent observing stations. From this 
point of view, one of the most important 
functions of synoptic analysis is the construc- 
tion of a “representative” or “smoothed” state, 
which is relatively insensitive to nonsystematic 
instrumental and round-off errors and which 
does not reflect irregularities that cannot be 
reconstructed in detail. 
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The analyst’s subjective process of smoothing 
may be regarded as a rather complicated and 
ill-defined type of averaging. That is to say, 
he bases his analysis at any given reference 
point not only on observations in the imme- 
diate vicinity, but also on observations at a 
number of points in the surrounding area, 
attaching greatest weight to the data nearest 
to the reference point, less weight to data at 
points farther away, and no weight at all to 
observations at great distances from the ref- 
erence point. This suggests that the analyst’s 
“smoothed” state is roughly equivalent to a 
running area average, in which the raw data 
at a number of points in the surrounding area 
are weighted according to distance from the 
point at which the average applies. Viewed in 
this light, the desirability of “smoothing” is 
immediately clear. Owing to the nonsystematic 
character of round-off errors and of some types 
of instrumental error, the average of a large 
number of reported measurements, taken at points 
spaced regularly over an area surrounding a 
given reference point, generally differs from 
the corresponding true average by an amount 
less than the error in a single reported measure- 
ment at the reference point. Thus the “smooth- 
ed” state, constructed by averaging the original 
data, suffers less from outright error and more 
accurately characterizes a “true” state than does 
the unaveraged state revealed by the raw ob- 
servations. 

Any process of averaging, however, also 
tends to obliterate small scale fluctuations 
which, although their detailed structure is nei- 
ther known nor observable, are nevertheless 
quite real. At first glance, it might appear that 
the resultant loss of detail might more than 
overbalance the advantages of reducing non- 
systematic error. On the other hand, the density 
of reporting stations independently places a 
definite upper limit on the resolving power of 
our system of meteorological measurement. 
This, in fact, provides a clue as to how strongly 
‘the state should be averaged. In general, to 
minimize the errors of reported observations 
without sacrificing what little resolving power 
we do possess, the observed state must be 
averaged in such a way as to suppress fluctua- 
tions whose scale is comparable with or less 
than the distance between neighboring ob- 
servation stations, but leaving disturbances of 
much larger scale essentially intact. Even so, 
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any consistent system of averaging will result 
in some loss of detail; however, this is simply 
the price that must be paid for a more accurate 
and “representative” characterization of the 
state. 

A second, and more serious objection to the 
characterization of a state by averaged variables 
is that the Navier-Stokes equations apply “in 
the small” — that is, they apply to the unaverag- 
ed state, in which fluctuations above the molec- 
ular scale are not smoothed out, rather than 
to an averaged state. For this reason, it is not 
strictly correct to apply the classical equations 
of hydrodynamics to a “smoothed” analysis, 

owever artistically drawn, without modifying 
them for the effects of fluctuations that have 
been averaged or smoothed out. Granting, 
therefore, that it is desirable to carry out some 
kind of averaging process on the variables that 
characterize the state of the atmosphere, it is 
equally desirable that we derive from the 


hydrodynamical equations a new system of 


equations that is specially designed to apply to 
the averaged variables. 

At this point, it is relevant to note that the 
meteorologist’s situation is very much akin to 
that which motivated the development of the 
kinetic theory of gases and the Reynolds theory 
of turbulence. In the former case, owing to the 
impossibility of observing the position and 
velocity of every molecule of a gas, the physic- 
ist has been forced to characterize the state of 
a gas by certain statistics of its molecular state, 
and has derived from the Newtonian collision 
equations a closed system of hydrodynamical 
equations that apply to such quasi-continuous 
statistics as fluid pressure, fluid density (mean 
concentration of molecules),and fluid velocity 
(mean molecular velocity). For similar reasons, 
Reynolds proposed averaging the fluid state, 
in order to smooth out “turbulent” Huctuations 
whose scale or period is too small for observa- 
tion of their structural details, but whose scale 
is larger than that of molecular motion. As 
in the kinetic theory of gases, the central prob- 
lem of the Reynolds theory of turbulence is 
to derive a closed system of equations that 
apply to an averaged state. 

According to Reynold’s formal definition, 
what constitutes “turbulence” clearly depends 
on the mode of averaging-e. g., on the effec- 
tive size of the domain of averaging and, ac- 
cordingly, on the scale of the fluctuations that 
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are “smoothed out”. As far as this definition 
goes, “turbulence” is not a real and distinctive 
physical phenomenon, but is simply a formal 
and rather useful concept for dealing with fluid 
systems whose state is not completely known. 
From this same point of view, the turbulence 
problem, stated as an initial value problem, is 
essentially equivalent to the problem the mete- 
orologist faces in trying to predict the state of 
the atmosphere with incomplete knowledge 
of its initial state. 

At the same time, it must be recognized 
that the problem of deriving a system of 
equations for the synoptic analyst’s “smoothed” 
or “representative” state does differ from the 
turbulence problem in at least one real and 
fundamental respect,namely, in that fluctuations 
whose scale is intermediate between the scales 
of the longest “planetary waves” in the mean 
flow and of fluctuations at the threshold of 
detectability are not completely smoothed out 
by a systematic process of averaging, and are, 
therefore, reflected in the mean state, though 
with reduced amplitude. The theory of tur- 
bulence, on the other hand, aims at deducing 
the interaction between the mean motion and 
the fluctuations that are smoothed out by 
averaging, from only the knowledge of a mean 
state that does not directly reflect any of the 
structural details of the fluctuations nor, for 
that matter, their very existence. To point up 
this distinction, it is useful to return to the 
example of the kinetic theory of gases, whose 
accomplished objectives are close to the un- 
realized ideals of the theory of turbulence. In 
this case, the molecular motions or fluctuations 
cannot be observed at all, whence the details of 
the motion cannot enter directly into the de- 
termination of the mean state. Nevertheless, 
it is possible to deduce a complete system of 
equations in which certain statistics of the 
molecular motions are the dependent variables. 
Fortuitously enough, those variables may be 
identified with quantities that can be measured 
directly, without knowledge of the details of the 
molecular motion. In this respect, the derivation 
of equations for a weighted spatial average of 
the observed motions falls short of the ideal of the 
theory of turbulence, and, in that sense, does 
not represent a fundamental approach to the 
turbulence problem. From the standpoint of 
examining the interaction between the averag- 
ed motions and the fluctuations that are 
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smoothed out, however, it is still of some inter- 
est to derive the equations of mean motion 
corresponding to a given mode of averaging, 
whether the unaveraged motions are actually 
observed or not. 
To provide a concrete basis for further discus- 
sion, we shall briefly review certain points in 
the pattern of development of the classical 
theory of turbulence. For simplicity, we shall 
consider the motion of an incompressible auto- 
barotropic fluid relative to a nonrotating coor- 
dinate system, with a view to predicting its 
mean motion. The unaveraged state of motion 
is governed by the well-known equation for 
the conservation of the vertical component of 
relative vorticity, taken together with the con- 
tinuity equation. 

a 


al ve Vv =O 


V :V=0 


in which v is the projection of the total vector 
velocity on a horizontal plane, v is the hori- 
zontal vector gradient, and the vorticity & is 
the magnitude of the curl of v. Following 
Reynolds’ development, we next average the 
equations over a finite interval of time centered 
on time ¢, denoting an averaged quantity by a 
bar placed above it. Since the operation of 
averaging is commutative with differentiations 
with respect to time and the horizontal coor- 
dinates, 

at = 

Er + V -Êv=o 


V-V=0 


At this stage, it is customary to assume that a 
mean varies much more slowly than deviations 
from the mean. Thus, denoting a quantity’s 
deviation from its mean by a “primed” quan- 
tity, 


ee, 


Finally, by introducing the stream-function y, 
the system governing the mean motions can 
be put in the form 


Ho = - = 

Fe: VE+YV Da =O (x) 

v=K x vy C= v2y (2) 
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where K is a unit vector directed vertically 
upward. The first of these equations simply 
states that the rate at which the mean vorticity 
of an element, travelling with the mean veloc- 
ity, is created or destroyed is the net rate at 
which vorticity is transported in or out of a 
neighborhood of that element by the “eddies” 
that have been smoothed out. 


It has been realized for some time, of course, 
that Egs. (1) and (2) are not satisfactory from 
the standpoint of predicting the mean motions, 
beginning only with a given initial distribution 
of the mean vorticity £. Aside from the fact 
that the Reynolds’ mean is not well defined 
(being regarded as variable for some purposes 
and constant for others), Eqs. (1) and (2) do 
not form a closed system, for the eddy transport 


¢’v’ cannot obviously be expressed in terms of 


£. This problem that of deducing the form of 
the eddy transport is, in fact, the key problem 
in the Reynolds theory of turbulence. Simul- 
taneous with the solution of that problem 
would come the means for predicting the 
averaged state of certain simple fluid systems, 
whose behavior approximates that of the 
atmosphere. Up to the present, however, the 
most strenuous efforts have not succeeded in 
reducing Reynolds’ eddy transport to the desir- 
ed mathematical form. 


On the other hand, it is conceivable that a 
closed system of equations for an averaged state 
might be derived by applying an equally use- 
ful, but entirely different mode of averaging. 
For the sake of argument, let us consider a 
running area mean, in which the values of the 
unaveraged quantity at regularly spaced points 
over the infinite plane are weighted less and 
less with increasing distance from the point at 
which the average applies. Now, the mean 
value associated with any given point is simply 
a linear combination of the unaveraged values 
at all points in the plane. Thus, corresponding 
to each point in the plane, one may write 
(schematically) one linear algebraic equation, 
in which the unknowns are the unaveraged 
values at all points in the plane, the coefficients 
are the values of the weighting factor at corre- 
sponding points, and the non-homogeneous 
term is the known mean value at the point in 
question. Since there are as many points as 
unknowns, the totality of such equations forms 
a closed system of linear algebraic equations 
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which, in principle, can be solved by inverting a 
determinant of infinite order. The consequence 
of this argument is that the unaveraged distribu- 
tion is uniquely determined by the averaged 
distribution, and can be reconstructed from it. 
This is evidently a sufficient, if not necessary 
condition for being able to express the devia- 
tion from a weighted mean in terms of the 
mean itself, from which it follows that the 
“eddy transport” can be expressed in terms of 
the mean motion. 

The argument outlined above is, of course, 
a purely formal one, and is not intended to 
suggest that the unaveraged distribution should 
actually be recovered from an averaged distri- 
bution by solving a system of infinite order. 
As it turns out, it is possible to define an averag- 
ing process such that the unaveraged distribu- 
tion can be reconstructed by repeated applica- 
tion of an integral operator. 

Before setting out on a line of attack that 
differs radically from Reynolds’ classical ap- 
proach to the problem, it is appropriate to 
comment on the relative merits and disadvan- 
tages of alternative modes of averaging. To 
begin with, it is probably fair to say that any 
type of average that displays statistical stability 
is suitable for purposes of reducing nonsys- 
tematic error. Beyond that, the choice of a 
mode of averaging appears to depend primarily 
on experimental conditions or systems of meas- 
urement, and on the desired form of the results. 
It is quite probable, for example, that Reynolds 
chose to average time-wise simply because the 
most economical system of laboratory meas- 
urement is capable of greater time resolution 
than space resolution. The system of meteoro- 
logical measurement, on the other hand, re- 
solves about as well in space as in time. More- 
over, since the problem of weather prediction 
is essentially an initial value problem, it is most 
natural for the meteorologist to average space- 
wise. There is one other consideration that 
enters into the choice of a mode of averaging, 
namely, the distinctive form of the equations 
for the unaveraged motion. From this stand- 
point, an area mean is preferable to the Rey- 
nolds time mean, since the vorticity equation 
for nondivergent or quasi-geostrophic flow is 
distinguished by the way in which space de- 
rivatives enter the equation. 

Having dwelt at some length on background 
and motivation, we can now state very simply 
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the problem that takes up the remainder of 
this paper. The main purpose is to develop a 
closed system of equations for the averaged 
motions of certain idealized fluids, whose be- 
havior is, in several crucial respects, quite 
similar to that of the middle and upper troposp- 
here. For clarity and simplicity, we shall fırst 
consider the motions of an incompressible auto- 
barotropic fluid, later indicating how the meth- 
ods for handling that case can be extended to 
the case of divergent barotropic flow and simple 
types of baroclinic flow. By judicious choice of 
a mode of area-wise averaging, it will be found 
possible to reconstruct the unaveraged distribu- 
tion from a given averaged distribution. A 
further consequence of the fundamental opera- 
tional property of the mean is an algorithm 
for inverting operators of the Laplace and 
Helmholtz type. Taken together, these results 
provide the ingredients of a method for expres- 
sing all quantities in terms of mean values that 
are either initially known or predictable. 

Finally, for purposes of orientation, the re- 
sults of this study are compared with one of the 
principal conclusions of the classical mixing- 
length theory of eddy diffusion. By introducing 
the mixing-length hypothesis without the usual 
restrictions on the statistical properties of turbu- 
lent fluctuations, it is shown that the eddy 
transport is directed normal to the gradient 
of mean vorticity, rather than down the gra- 
dient, and that the coefficient of eddy transfer 
depends on the correlation between eddy dis- 
placements in one direction and eddy velocities 
in a direction normal to it. This clearly indicates 
that the process by which mean vorticity is 
generated or destroyed is not the classical pro- 
cess of eddy diffusion, and that assumptions 
about the statistics of turbulent fluctuations 
should not be introduced without regard to 
the nature of the quantity whose unaveraged 
value is conserved. 


2. The Weighted Area Mean as a 
“Smoothing” Operator 


Asa preliminary to introducing mean values 
into the hydro-dynamical equations, we shall 
first discuss the general mathematical properties 
of a weighted area average which, as will be 
seen later, is specially adapted to the form of 
those equations. 
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The unaveraged variables under considera- 
tion are, in general, continuous and continously 
differentiable functions of the rectangular coor- 
dinates x and y in an infinite “horizontal” 
plane, and time 1. Corresponding to the spatial 
distribution of a variable at a fixed time, we 
shall define a mean value of that variable, 
denoted by placing a bar over the unaveraged 
quantity. 


ais N= Xf {Kone mdr (3) 


where @ is a typical function of x and y, a 
is a real constant, Ky is the zero order Bessel 
function of the second kind with imaginary 
argument, € and 7 are the coordinates of the 
variable point (dummy variables of integration 
corresponding to x and y), and r is the distance 
between the fixed point (x, y) and the variable 
point (& 7) 1e, r?=(«-€)2 + (y—m)?, The 
indicated area integration extends over the 
entire plane. It must be kept in mind that x 
and y, the coordinates of the point at which 
the average applies, are constant with respect 
to the integration. 

With regard to the properties of this average, 
we first note that it is normalized, in the sense 
that the mean value of “one” is still “one”. It 
follows that, if y is finite over the entire plane, 
the average is always convergent, for the abso- 
lute magnitude of the average cannot exceed 
the maximum unaveraged value without re- 
gard to sign. The weighting function K)(«r) — 
or more simply, K-decreases monotonically 
with increasing distance from the point (x, y) 
and is never negative. It has a logarithmic 
singularity at (x, y) and approaches zero like 

1 


(ar) “e-% with increasingly great distance from 
(x, y). Thus K satisfies the essential conditions 
on a proper weighting function, namely, that 
it never gives negative weight and places less 
and less weight with increasing distance from 
the point at which the average applies. 

The constant «, which is unspecified in value 
and which we are free to fix, deserves special 
comment. From the foregoing discussion, it 
can be seen that, as & is increased, relatively less 
weight is attached to points at a large fixed 
distance from the point (x, y), relative to the 
weight attached to points at a somewhat smaller 
distance. Thus the value of & determines the 


“effective” size of the domain of averaging, 
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and thereby fixes the scale of transition between 
the small-scale fluctuations that are smoothed 
out by averaging and the large scale disturb- 
ances that are left essentially intact. In that 
sense, one might speak of this smoothing proc- 
ess as “averaging of modulus &”, and define 
“turbulence of modulus &” as the fluctuations 
that are smoothed out by that mode of averag- 
ing. In any case, “turbulence” in the formal 
sense must be defined with reference to a 
particular type of average. 

To illustrate the role played by the “scale” 
parameter «, we shall apply the average to a 
function that is independent of y and varies 
sinusoidally in the x-direction. By a straight- 
forward integration, it can be shown that, if 


p=cos px, 
9= (5) p (4) 


Two extreme cases may be distinguished, 
according as « < B or «> f. 


lim = lim w=0 
Ha 2.0 


This result may be interpreted from two dif- 
ferent points of view. Relative to a fixed 
value of x, disturbances of very small wave 
number (large scale) are virtually unaffected 
by the process of averaging, while disturbances 
of very large wave number (small scale) are 
completely obliterated. Similarly, if the un- 
averaged distribution is a continuous spectrum 
of wavelike components, the average distribu- 
tion will reflect only those components whose 
wave numbers are comparable with or less than 
the scale parameter «. Relative to disturbances 
of a fixed scale, on the other hand, averaging 
for a very small value of « tends to smooth 
out the disturbances, whereas averaging for a 
very large value of « leaves them almost intact. 
Accordingly, we shall later speak of the case 
where « is much greater than some character- 
istic wave number as “weak averaging” and, 
where it is much less, as “strong averaging”. 


3. The Fundamental Operational Property 
of the Mean 


Although the “‘bar’’-operator has been defin- 
ed as a weighted average, it also has many of 
the features of a transform or integral operator, 
and it will frequently be convenient to think 


Tellus VI (1954), 2 


155 


of it in that light. Our next concern is to 
investigate some of the operational properties 
of the mean, regarding it as a special type of 
double transform. 

The fundamental operational property of the 
mean, from which all others are deducible, 
stems from the peculiar character of the weight- 
ing function K. It may be verified by direct 
differentiation that K satisfies the Helmholtz 
equation 


V?K= 0K (5) 


The usefulness of this relation is most easily 
demonstrated by applying the bar-operator to 
the Laplacian of a typical variable @. According 
to the definition of Eq. (3), 


za 2 
ie Kv 2ode&dn 


By making use of well-known vector identities, 
the equation above can be rewritten in the 
form 


We next substitute into the equation above the 
value of v?K given by Eq. (s), and again 
introduce the definition of the bar-operator, 
with the result that 


vp = cp + 


2 
+ lim =f fe (KY p-gv K)-dédy. 
SR>0 S 
Loo 
where the domain of integration S is an annular 
region bounded by a circle J’ of radius R and 
a circle y of somewhat smaller radius o, both 
with center at (x, y). Thus, as R becomes 
infinite and © approaches zero, the domain S 
will cover the entire plane. This limiting proc- 
ess is necessitated by the singularity of K at 
my. oe HORS 
It remains to evaluate the limit indicated in 
the equation above. Since @ has no singularities 
at all and K has no singularities in S, the area 
integral over S may be transformed by Gauss’ 
theorem into a line integral taken around a 
path C which encloses S. The path of integra- 
tion is traversed by proceeding around J” in a 
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counter-clockwise direction, along the left hand 
edge of a “cut” adjoining I’ and y, around y 
in a clockwise direction, back along the right- 
hand edge of the “cut”, and so on around I’ 
to the starting point. Thus 


he a OR 
29 = u? = dC 
ie rn =f (xz oF) 
f R—00 c 
Loo 
- à 
in which ce and “as are the projections of 


on on 
vo and VK on a unit vector directed outward 
and normal to the curve C, and dC is an in- 
crement of length along C. It is evident, how- 
ever, that the contributions along one edge of 
the “cut” and along the other edge exactly 
cancel, so that 


27 

in 28 f (120 A 
lim an (X ee a 
o>0 0 


where © is the angle between the x-axis and a 
radius from (x, y). The second term on the 
righthand side of the equation above vanishes 
as R becomes infinite, since RKy(«R) and 
RK, («R) approach zero exponentially as R 
becomes infinite. The third term on the right- 
hand side may be evaluated by the techniques 
developed to prove Catichy’s residue theorem. 
Owing to the singularity of K at (x, y), it 
yields a finite contribution proportional to the 
value of pat the point (x, y) where the average 
applies. The final result, in fact, is simply that 


v?p = a(p — 9) (6) 
This is the fundamental operational property of 
the mean. Taken as it stands, it carries the 
rather remarkable implication that the devia- 
tion of a variable from its corresponding local 
mean is proportional to the mean of its Lapla- 
cian. 


4. The Commutative Properties of the 
Bar-Operator 


There are a number of other convenient 
operational properties of the bar-operator that 
follow directly from the fundamental property 


PHILIP DUNCAN THOMPSON 


given in Eq. (6). For example, applying the 
bar-operator to both sides of Eq. (6), 


v2 = a(p =) 
We also apply the bar-operator to 729, which 
is given by Eq. (6) as 

vip = 02 (p — g) 
Subtracting the second of the two equations 
above from the first, 

v29 — V29 =0 
Now, it can be shown that the weighted 
average of a function can be everywhere zero 


if and only if the function itself is identically 
zero. In particular, 


7p = V2 


Thus, the bar-operator may be commuted with 
the Laplacian operator. 

With this result, it is a relatively easy matter 
to show that the bar-operator can be commuted 
with all vector differentiations. That is to say, 


Vp = VE 
VERS 
7 OA = oon 


for all scalar variables 9 and all vectors A. 
Moreover, since the order of integration with 
respect to the space coordinates and differentia- 
tion with respect to time can be reversed, the 
bar-operator is commutative with time dif- 
ferentiation. 


ap _ Ap 
dt ot 
As will be found later, all these properties are 


useful in manipulating the equations of mean 
motion. 


5. Two Inversion Formulae 


By a simple rearrangement of Eq. (6), the 
fundamental operational property of the mean 
can be put in the form 


ee 
RL oe GO 
Denoting the bar-operator ® ) by L. ()fand 


successively applying the operator L ( ) to the 
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equation above (n—1) times, we generate the 
following array of equations. 
I 


gy — L(y) + 52 L(V) = 0 


Lip) - L(g) + GL(9%)  =o 


I 
© 


1 2 
L*-2(g) En 1(g) +, L*-*( 929) 


it 
Lr =1(o) =] # (œ) ae 5 L"(v2@) =o 


We next add all these equations, noting that 
the second term of the ith equation just cancels 
the first term of the (i+1)st. Thus, 


n 


I NDR 
ee > Li(v*9) 


i=1 


Finally, we shall consider the limit of the 
expression above as n becomes infinite. At the 
same time commuting the Laplacian and L ( ) 


Operators, 
Sk Î v2Li( ( ) 
eg eae v?L*(@) 7 
i=1 


where g* is the limit of L,(p) as n becomes 
infinite. 

Our next concern is to examine the be- 
havior of L"(g) for increasingly large values 
of n. Now, it is intuitively obvious that the 
result of applying a normalized averaging proc- 
ess to a function infinitely many times will 
tend to a definite limit. That is to say, q* 
exists. A necessary condition for the existence 
of g* is that one more application of the bar- 
operator leaves the limit unchanged -i. e., 


Sk 

PO 
But the fundamental operational property of 
the mean, coupled with its commutative prop- 
erties, implies that 


ik 


v2p* = a2(g* — p*) 


a 
from which it follows that 

v2p* = 
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This is simply Laplace’s equation, whose solu- 
tions are well known. By hypothesis, @* is 
nonsingular and becomes infinite at infini 
more slowly than an exponential function. The 
only function g* that satisfies both Laplace’s 
equate and the condition above is of the 
orm 


pr =ax+by+c 


where a is the arithmetic average of ay , taken 


Ox 

over the entire plane, b is the corresponding 
Op 
dy 
depending on the location of the origin. With 
reference to meteorological variables, the arith- 
metic mean over a very large area is highly 
persistent, so that @* may be regarded as 
known for purposes of short range prediction. 

A further consequence of the argument out- 
lined above is that 


average of —, and c is an arbitrary constant 


lim v?Li(g) =o 

i>oo 
For this reason, and because each application 
of the L ( ) operator involves multiplication 
by a factor of «?, the series on the righthand side 
of Eq. (7) converges for sufficiently small «. 
It is demonstrable, therefore, that the limits 
indicated on the righthand side of Eq. (7) 
exist. Writing Eq. (7) in slightly different form, 


I nA 
a er: > Li(v2q) (8) 
i=1 


This equation may be regarded as an inversion 
formula for solving equations of the Poisson 
type. That is to say, given the distribution of 
v2p over the entire plane, y can be computed 
from it by repeated application of the integral 
operator L ( ). Stated in this way, the inversion 
formula is reminiscent of a result obtained by 
Fyortorr (1952), through a formal expansion 
of the finite-difference form of the Laplacian 
operator. 

A more striking consequence of the opera- 
tional properties of the mean can be derived by 
applying the Laplacian operator to both sides 
of Eq. (8). Noting that v?p*=o, and letting 
VU D, 
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It must be borne in mind that the equation 
above is an identity, which holds for all func- 
tions ®. One may, therefore, define an “un- 
averaging” operator L-!( ) which, when ap- 


plied to the averaged distribution ®, yields the 
unaveraged distribution ®. That is, according 


to Eq. (9), 


Eee DE 


i=0 


In this respect, Eq. (9) is an inversion formula 
for computing the unaveraged distribution cor- 
responding to a given averaged distribution. 
By following a line of reasoning similar to that 
sketched out in the preceding paragraph, it can 
be shown that the summation on the righthand 
side of Eq. (9) converges to a definite limit. 

In order to illustrate how the unaveraged 
distribution is to be computed from the mean, 
and to examine the conditions under which the 
computational procedure converges rapidly, 
we shall return to the example of Eq. (4). In 
that case, Eq. (9) reduces to 


is SS OP ONE eset 
Looe a2 + B? foe” 
i=0 


By summing the geometric series exactly, it 
is easily verified that the equation above is an 
identity. It is also evident that the series con- 
verges for all values of & and ß, and that it 
converges rapidly for values of « much less 
than the wave number. The latter case, which 
corresponds to “strong averaging”, is the one 
of greatest interest in connection with the 
development of equations for the mean motion. 


6. The Equation of Mean Motion for an 
Incompressible, Auto-barotropic Fluid 


Up to this point, our discussions have dealt 
only with the mathematical properties of a 
particular type of average -i. e., with mathe- 
matical definitions and identities. Our next 
concern is to apply those mathematical methods 
and results to the problem of deriving a closed 
system of physical equations for the averaged 
motion of a fluid. 

In order to present general features of method 
in the simplest possible form, and to facilitate 
comparisons with the classical theory of atmos- 
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pheric turbulence, we shall consider the mo- 
tions of an incompressible auto-barotropic fluid 
relative to a non-rotating coordinate system. 
This, of course, is the model discussed earlier 
in the introduction, in connection with Rey- 
nolds’ formulation of the turbulence problem. 
The unaveraged motions of the fluid are govern- 
ed by the following closed system of equations. 


K4vıri=o 
v=Kx vy C= Vp 


where v is the horizontal component of veloc- 
ity, & is the vertical component of relative 
vorticity, @ is the streamfunction, and K is a 
unit vector directed vertically upward. To 
obtain the equations for the averaged motion, 
we simply apply the bar-operator to the equa- 
tions above. 


Ovo = (10) 
v=Kx vy (11) 


As was also pointed out earlier, the key problem 


is to express the mean vorticity advection v- VE 


in terms of the mean vorticity £. If it is possible 
to do so, then Eq. (10) provides a method for 


predicting the future distribution of ¢ from 
its known initial distribution. 

For later purposes, it will be convenient to 
decompose the mean vorticity advection by 
introducing the deviation from the mean. De- 
noting the deviation of a variable from its 
mean by a “prime”, we may state the follow- 
ing general result, which is a direct consequence 
of the fundamental operational property of the 
mean. [Eq. (6)]. 

/ I Fe 
Pie Sap 


x? 


Applying this result to the stream-function, 


y= Hess = 
BOOS 22° 

; I 
VE HH XG 

I — 
rt - 
Sul 

a2 
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The expressions above are now introduced into 
the vorticity advection. 


BD yvl=v.vitvV.vötv.vOrv.ve- 
re 2F ra F 
SAPS a A et) 


in which J is the Jacobian, or “functional” 
determinant. Finally, averaging the equation 
above, again making use of Eq. (6), and sub- 
stituting into Eq. (10), 


ea 27\ 
ET IE 


Bi ee == — 
et 0 (13) 
Formally, at least, this equation represents the 
solution of the problem that was originally set. 
Starting with the known initial distribution 
of ¢, we may compute y by making use of 
Eq. (12) and the inversion formula for solving 
equations of the Poisson type [Eq. (8)]. Having 
computed y, it is then possible to compute 


v from Eq. (11). In this way, since they involve 
no time derivatives, all terms of Eq. (13) except 
the first can be computed from the initial distri- 


bution of ¢. Thus Eq. (13) gives the instan- 


taneous rate of change of © at some arbitrarily 
chosen initial moment, in terms of the spatial 


distribution of £ at that time. This property 


enables us to predict the distribution of ¢ at a 
time slightly later than the initial moment, 
whence the process can be repeated indefinitely. 
This is sufficient to demonstrate that Eq. (13) 
is a true prognostic equation, and that Eqs (10), 
(11), and (12) constitute a closed system of 
equations for the mean motion. 


7. The Case of “Strong Averaging” 


Although Eq. (13) comprises a satisfactory 
solution to the problem that was originally 
set, it is difficult to interpret this result as it 
stands. For simplicity and ease of comparison, 
we shall next consider the case of “strong 
averaging’ -i. e., averaging for which the 
scale parameter « is much less than some charac- 
teristic wave number. In this case, the weighted 
average of the vorticity advection displays the 
same properties as an unweighted average, and 
takes on a more familiar form. 

By introducing the relationship between the 
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mean and deviations from the mean [Eq. (6)], 
the product of any two typical variables u and 
v may be written as 

uv = uv + (uv + vu’) + u'v’ 


1 SIE a Ne 
=uv —— (uv9v + vv2u) +— v2uV®v 
œ ci 


Le 2 -— — I _— — 
=uv—— V*(uv) + Vu. vv+—v2uV2V 
04 & he 


where, as before, a prime denotes the deviation 
from the mean. Applying the bar-operator to 
the equation above, and again making use of 
the fundamental operational property of the 
mean, 


DIT = I ni 
uv=uv+ — vu-vv+— v2uv®v (14) 
a at 


For purposes of estimate, we now ascribe to 
the disturbances in the mean flow a characteristic 
wave number »!. Thus the relative orders of 
magnitude of the second and third terms on 
the righthand side of Eq. (14) are 

T—- 2 

2 VAIO x 

== yt 

er Veuve — = 
In general, therefore, if « is much less than the 
characteristic wave number of the disturbances 
in the mean flow (the longest wave components 
in the unaveraged spectrum), the third term 
on the righthand side of Eq. (14) is much 
greater than the second. In other words, in the 
case of “strong averaging”, the weighted aver- 
age of a product reduces to 


uv = uv+ uv' (15) 
This property, it will be recognized, is peculiar 
to the usual unweighted mean or, more gener- 
ally, to a “weakly” weighted mean. 

Specializing the general result of Eq. (15), 
we may now put Eq. (13) in the alternative 
forms 


at FE = (al? 
2); v-ve= -v Tv 4 
un fi 
& 
Ir 
= ~ =a JK6; 0) = 
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As pointed out earlier, the conditions under 
which the equations above apply are also those 
for which the series in the inversion formulae 
(8) and (9) converge most rapidly. From this 
point onward, the case of “strong averaging” 
is the one with which we shall be primarily 
concerned. 

The righthand side of Eq. (16) represents 
the rate at which the mean vorticity of an 
imaginary element, traveling with the mean 
velocity, is created or destroyed or, from an- 
other point of view, reflects the interaction 
between the mean flow and the fluctuations 
that have been smoothed out by averaging. 
According to Eq. (16a), mean vorticity is 
generated by a net convergence of the “eddy” 
transport of vorticity —i. e., the transport of 
vorticity effected by the “eddies” or fluctua- 
tions that have been smoothed out. From the 
distinctive form of Eq. (16 c), it may be seen 
that no mean vorticity is created or destroyed 
if the unaveraged vorticity distribution is lined 
symmetrical, a fact that has been independently 
and adequately established by earlier writers. 


8. The Relation of the Equation of Mean 
Motion to the Theory of Turbulent 
Diffusion 


According to 16a, the equation of “strongly 
averaged” motion can be written as 


Pee ey =O 
Er 4 % 

This equation, it will be recognized, is identical 
in form to the Reynolds-averaged vorticity 
equation. Our present purpose Is to compare 


the form of the eddy transport ¢’v’ — which, 
in this case, can be deduced exactly - with the 
form derived in the classical theory of eddy 
diffusion. Following the usual method of devel- 
opment, we introduce the “mixing-length” 
hypothesis, imagining that turbulent fluctua- 
tions of vorticity are brought about by small 
displacements of fluid elements from their mean 
positions, each conserving its original mean 
vorticity. 


(ly! = - (r'- ve)’ 


wherer’ is to be interpreted as an eddy displace- 
ment. Subject to several restrictions, including 
the assumption that eddy displacements in one 
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direction and eddy velocities in the direction 
normal to it are uncorrelated, the expression 
above can be transformed into 


Be 


This result simply states that the eddy transport 
is directed down the gradient of mean vorticity, 


and that the coefficient of eddy diffusion r’ - v’ 
depends on the correlation between eddy dis- 
placements and eddy velocities in the same 
direction. If, as is usually supposed in the classi- 
cal theory of eddy diffusion, the latter correla- 
tion is always positive, then extreme values of 
mean vorticity would tend to disappear as time 
went on. 

In the present case, on the other hand, there 
is no real compulsion to introduce restrictions 
on the statistics of the eddy motions, for the 
eddy velocity v’ is expressible in terms of the 
mean vorticity gradient and vice versa. Thus, 
according to the results of section 6, 


= ! 


Cy’ = — (r-vé)v’ 


= —(K-r’ x v’) (Kx vi) 
or, approximately, 
Cv’ = -(K-r’ xv) (Kx vo) 


In contradistinction to the classical theory, this 
result states that the eddy transport is directed 
normal to the mean vorticity gradient, and that 
the coefficient of eddy transfer depends on 
the correlation between eddy displacements in 
one direction and eddy velocities in the direc- 
tion normal to it. 

It should be emphasized that it is not even 
necessary to introduce the mixing-length hy- 
pothesis in the present case, for the exact form 
of the eddy transport of vorticity is implicit 
in Eq. (13). The comparison above is, however, 
useful in demonstrating that the process by 
which mean vorticity is generated or destroyed 
is not one of pure eddy diffusion, and in point- 
ing out the dangers of imposing arbitrary re- 
strictions on the statistics of the eddy fluctua- 
tions, however plausible they might otherwise 


appear. 


9. Extensions to Simple Types of Baroclinie 

Flow 
The foregoing results, it will be realized, 
apply to a very special type of flow, and are 
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intended mainly to sketch out the pattern that 
a more general development might follow. 
Some generalizations-e. g., the inclusion of 
effects engendered by referring the motions to 
a rotating coordinate system require only 
slight modifications of the preceding theory 
and need no elaboration. The extension of this 
theory to systems whose behavior is controlled 
by basically different physical mechanisms is 
not, however, so obvious. Our present concern 
is to show how the methods developed in the 
foregoing sections can be used to derive equa- 
tions for the mean quasi-geostrophic motions 
of a simple baroclinic Auid. Divergent baro- 
tropic flow is included under the latter as a 
special case. 

The equations for the simple types of baro- 
clinic flow recently studied by Errasson (1952), 
Eapy (1952), THOMPSON (1953), BusHBy and 
SAWYER (1953), and CHARNEY and PnirLıps 
(1953) can, in general, be put into the form of 
a simultaneous system of two equations. 


IP - BP) = NT, var) ext AJ(®, v2) = 


(17 a) 


= (720 - 120) = - J(P, 720) - J(®, v2P) - 
- BJ(®, 928) + „AJCP, ®) BTE (176) 


Here the dependent variables Y and ® are the 
stream-functions for the vertically-integrated 
mean wind and thermal wind, respectively. 
The quantities A, B, u,, and u, are regarded 
as known constants, ß is the northward deriva- 
tive of the Coriolis parameter, and the coordi- 
nate x is directed toward the east. Since the 
right hand sides of Eqs (17a) and (17b) do not 
involve time derivatives, it is evident that they 
can be expressed in terms of Y and ®, by 
making use of Eq. (6). Accordingly, if the 
bar-operator is applied to the equations above, 
they may be regarded as a means for predicting 
(v2-u2)P and (v2 - w)®D. _ The remaining 
question is whether or not P and ® can be 
computed from the predicted values of (v? — 
u?) Wand (v?- w)®, to regenerate the initial 
conditions at a slightly later time. 

Stated in general terms, the crux of the 
matter is whether or not it is possible to com- 
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pute any typical variable » from a known 
distribution of (V?-u?)p. Although there are 
already well known methods for doing so, 
we shall present an integral operator method 
similar to the inversion procedure developed 
in section 5. As before, we generate an array of 
equations, obtained by successive applications 
of the L () operator to Eq. (6), and by 
multiplying the ith equation in the array by 


t (a2 —p2\i-1 
a? a2 ; 


(SS) [pow - (4 )ew]- 


I ao? = 2 tal 
= E ep) ( ae ) 


We next add all these equations, noting that 
the second term on the left hand side of the 
ith equation exactly cancels the first term of 
the (i+ 1)st. 


= a? — y? n 4 a Ni Er 
en (S34) 


L’(v’p—u?p) 
It will now be supposed that w is less than « 
—a restriction that is not very severe, since u 
is generally much less than the wave number 
of the large scale disturbances. In that case, 
the limiting form of the equation above is 


_ DENT 2 NE 
BES on Ben 


(18) 


This result shows that 9 can be computed from 
a given distribution of (V?-42)p. Since y is 
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less than «, the righthand side of the equation 
above converges more rapidly than the series 
in Eq. (8), which has already been found con- 
vergent. It can be concluded, therefore, that 
the methods developed in preceding sections 
will lead to a closed system of equations for 
the mean quasi-geostrophic motions of simple 
baroclinic fluids. 


10. A Weighted Time Mean as a Smoothing 
Operator 


In order to point up the advantages of the 
space mean, and to indicate how the foregoing 
methods can be extended to time means, we 
shall derive the fundamental operational prop- 
erties of a weighted time-average. In this 
connection, it should be noted that a time 
mean — if it is to be of any use for purposes of 
prediction — should not depend on unaveraged 
values at times later than a specified “moment 
of latest information.” Also, it is intuitively 
evident that the mean should place greatest 
weight on the state at the “moment of latest 
information”, with decreasing weight on the 
state in the increasingly remote past. This sug- 
gests an average of the form 

= t 

p(th=k f cp (rt) de (19) 

— CO 

where 7 is a dummy variable of integration 
corresponding to f, and k is an unspecified con- 
stant that we are free to fix. The average is 
normalized, in that the mean value of “one” 
is “one”. The exponential factor acts as a 
weighting function which attaches no weight 
on data beyond time f, greatest weight on data 
at time f, and decreasing weight to data at times 
in the increasingly distant past. It also satisfies 
the formal conditions on a proper weighting 
function, namely, that it is monotonic and 
never negative. 

Applying this averaging operator to a sinus- 
oidal time series, for example, we find that the 
average value of g=cos At is 


gy = sin © sin (At + ©) 


where 


O = sin” —— 

Vk? + 72 

so that, in this case, 
lim p= 


5 00 
4 


lim g=o0 
k 
——=-\0) 
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Thus the role of k is very similar to that of the 
“scale parameter” «, and what has previously 
been said of comparisons between x and the 
characteristic wave number is equally true of 
comparisons between k and the frequency. 
With regard to operational properties of the 
mean, it will be recognized that the “smoothing 
operator” defined by Eq. (19) is related to the 
Laplace transform, and that it has similar prop- 


erties. Applying the bar-operator to for 


example, and integrating by parts, 


em ke 2P = 

u ae 
B 

= klgekle 49 kt f ekle ddr 


— CO 


Whence, by definition, 


OP = 
— =k(p - 
= (9-9) 
This is the fundamental operational property of 
the time-mean. By direct differentiation of Eq. 
(19), it may be verified that 
dp _ Op 
dt ot 
Moreover, since the order of integration with 
respect to time and differentiations with respect 
to the space coordinates is reversible, the bar- 
operator may be commuted with all vector 
difterentiations. 

Eq. (20) enables us to derive a closed system 
of equations for the time-averaged motions, 
for deviations from the mean can be expressed 
in terms of the mean itself. That is, denoting 
deviations by primed quantities, 


(20) 


gene 
k ot 

We now introduce the expressions above into 
the time-averaged vorticity equation for in- 
compressible auto-barotropic flow. 


Tellus VI (1954), 2 


MEAN MOTIONS AND LARGE SCALE TURBULENCE 


at 

= EV Vo= 0 

Ot = 
After considerable manipulation, the equation 
above can be put in a form equivalent to that 
of the Reynolds-averaged equation. 


EME Er 
EN e+ gs (HF) =o a 


Since v=Kx vy and ¢ = vy, this equation 
involves only one unknown, -i. e., the mean 


stream-function y —and thus constitutes a clos- 
ed system of equations for the time-averaged 
motions. 

Although Eq. (21) states exactly and concisely 

Ow mean vorticity is generated by the fluctua- 
tions that are smoothed out in the averaging 
process. it must be noted that the time deriva- 
tive of the mean stream-function enters non- 
linearly. In general, therefore, the mean stream- 
function tendency cannot be obtained from 
Eq. (21) in closed form. On further reflection, 
it will be seen that this defect would be present 
in any system of time-averaging of this type. 
On the other hand, as was pointed out earlier, 
space-averaging leads to equations which are 
linear in the time derivatives and which can, 
accordingly, be solved by known methods. 


11. Summary and Conclusions 


The problem of deriving a closed system of 
equations for the mean motion of a fluid has 
been approached from the standpoint of re- 
constructing the unaveraged distribution of a 
variable from the distribution of its average. 
For purposes of mathematical simplicity, it has 
been found convenient to consider weighted 
averages, in which the weighting functions are 
chosen with regard to the peculiar form of the 
equations for the unaveraged motion. The 
averaging or “smoothing” operations associated 
with these weighted means are essentially trans- 
forms or integral operators, and have similar 
operational properties. 

By making use of the operational properties 
of a weighted area mean, we have derived an 
inversion formula |Eq. (9)| for computing the 
unaveraged distribution of a variable from its 
averaged distribution. 

An integral operator method [Eqs (8) and 
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(18) for inverting operators of the Laplace 
and Helmholtz types has been developed as a 
natural by-product. Finally, it has been shown 
that a certain type of weighted time mean has 
operational properties analogous to those of the 
area mean. 

To illustrate general features of method, the 
mathematical results of the first part have been 
used to derive a closed system of equations for 
the mean motion of an incompressible auto- 
barotropic fluid, relative to a nonrotating coor- 
dinate system [Eqs (13) and (21)]. This deriva- 
tion has been carried out for both area and 
time averages, either of which is satisfactory 
from the standpoint of reducing nonsystematic 
instrumental and round-off errors in reported 
observations. The equation for the area-averag- 
ed motion, however, is linear in the time deri- 
vatives and can be solved by known methods, 
whereas the equation for the time-averaged 
motion is nonlinear in the time derivatives and 
cannot be solved in closed form. This result, 
which stems from the distinctive form of the 
equation for the unaveraged motion, indicates 
that the area mean is preferable for purposes 
of applying the theory of mean motions. 

The case of “strong averaging’, correspond- 
ing to a suitably chosen value of an adjustable 
scale parameter, is of greatest intrinsic interest, 
and has been examined in some detail. In that 
case, the area-averaged vorticity equation is 
identical in form to the Reynolds-averaged 
equation. By introducing the “mixing-length” 
hypothesis without the usual restrictions on the 
statistics of eddy fluctuations, it is shown that 
the eddy transport of vorticity is directed nor- 
mal to the mean vorticity gradient, and that 
the coefficient of eddy transfer depends on the 
correlation between eddy displacements in one 
direction and eddy velocities in the direction 
normal to it. This result, of course, is at direct 
variance with the classical theory of eddy dif- 
fusion. 

Finally, we have outlined a general pattern 
of development by which these methods can 
be used to derive a closed system of equations 
for the area-averaged motions of a compressible 
barotropic fluid and of certain simple types of 
baroclinic fluids. It appears probable that more 
general methods for inverting the process of 
averaging will be capable of extension to a 
wide variety of physical problems. 
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The Exact Determination of the Effective Domain of Dependence 


of a One-dimensional Numerical Prediction Formula 


By LOUIS BERKOFSKY, Geophysics Research Directorate, Cambridge 


(Manuscript received April I, 1953) 


Abstract 


By assuming sinusoidal initial conditions, an exact expression for determining the 


€ 


“cut-off” 


point of a one-dimensional influence function is derived. From this result one may construct 
tables of percent error incurred by terminating the influence function at various distances. 


I. Introduction 


One-dimensional numerical prediction tech- 
niques usually require the integrated product 
of initial conditions and an “influence function” 
along a latitude circle. Fortunately, the nature 
of the influence function is usually such that 
one need not know the data all the way around 
a given circle of latitude. That is, in order to 
compute a sufficiently accurate forecast, one 
may terminate the influence function at some 
large, but finite distance from the origin of 
computations. This distance is referred to as 
the radius of the “effective domain of depend- 
ence’. As it turns out, the radius of the domain 
of dependence varies with the wavelength of 
the initial distribution, as well as with the phase 
difference between the initial distribution and 
the origin from which the forecast for a given 
point is computed. Thus, in making a forecast 
at a given latitude on a given day, it appears 
necessary to terminate the influence function 
at different distances for different forecast 
points in order to achieve the same level of 
accuracy everywhere. 

In a recent paper, THOMPSON (1952) has 
shown how to obtain the Fourier integra 
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solution of the Rossby vorticity equation in 
the form 
0 
z(rUt, t) =z(0, 0) + [z(E, 0) G(é)dé (x) 
where z(x, f) is the height of the soo-mb 
surface at time f, and 


GNS, <9) (à 


is the influence function. J, (2 = té) is the 
Bessel function of the first kind of order one, 
B is the N-S variation of Coriolis parameter, 
x is the distance, positive toward the east, 7 
is a constant, U is the undisturbed zonal wind 
speed. The essential difference between this 
solution and that of CHARNEY and ELIASSEN 
(1949) is that the latter solution is expressed as 
a Fourier series. Both solutions have the virtue 
that the “effective domain of dependence” is 
finite, so that one need not know the data 
everywhere to make a sufficiently accurate 
forecast. 

Lönngvıst (1949) has applied the Charney- 
Eliassen formula to obtain numerical predic- 
tions for sinusoidal initial conditions. Due to 
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the nature of the influence function in that 
formula, it appears that an attempt to carry 
out exact integrations for sinusoidal initial 
conditions would be an extremely difficult task. 
In the present study, however, it is possible, 
by using formula (r) above, to carry out the 
integration for sinusoidal initial conditions 
exactly, and thus compute the exact height 
change of the soo-mb surface. Actually, such 
a procedure is unnecessary for prediction of 
the propagation of sinusoidal initial conditions, 
since such waves will move with exactly the 
Rossby phase speed, suffering no change of 
amplitude or wavelength. Nevertheless, the 
performance of the exact integration yields 
information which is useful for determining 
quite accurate estimates of the “effective 
domain of dependence”. 

CHARNEY (1949) has derived an expression 
for determining the error incurred in cutting 
off the influence function at certain points. 
The integrals involved were determined 
empirically from an observed pressure distri- 
bution. In the present study, we evaluate the 
necessary integrals exactly by assuming sinu- 
soidal initial conditions. We assume that the 
results obtained in such a manner will apply 
sufficiently well to arbitrary initial conditions. 

We proceed in the following way: We first 
carry out the integration in (1) for sinusoidal 
initial conditions, and obtain an expression 
which is a function of the wavelength of the 
initial distribution and the phase difference 
between the sine wave and the origin from 
which the forecast is computed. For a given 
latitude and forecast period, we may then 
compute the height change as a function of 
wavelength and phase angle. We then perform 
the same integration, but this time terminate 
the integration at some finite distance from the 
origin. This integration yields an expression 
which is a function of distance out, as well as 
the previously mentioned parameters. Using 
the results of these two integrations, we are in 
a position to determine the finite distance out 
which will give a percent error equal to or 
less than some previously assigned maximum 
value. 

Thus, if we let 


2(6,0)=Asinml—(£-7)], (8x <o), (3) 


where A=amplitude, y= phase angle, m= 


? 


EOS, 18508) REI OURS IR NE 


L=wavelength, then we wish to satisfy the 
inequality 


J Asin m[ — (€ - y)] G(é)dé - 


J A sin m[- (&-)]G(&)dE 


t (4) 
- fAsinml-(£-7)]G(é)dE 
a AE, 
ee A sin ml - (€— 4) | G(é)dé | 


where E is the maximum percent error 
allowed. 


2. Evaluation of the Integrals 
Let 


B= % À sin m[ - (€- y)K - VBt) : 


— CO 


un (2V - BtË) 
V-E 


which, because (£, y < 0), may be written 


(5) 
dé, (E, x <o), 


B= - Ayßt f sin miele PE 
= vis 
(6) 


If we write 


y=2vßelel (7) 


we obtain 


B- AJ- sin (a2y2— K)J, (y)dy, (8) 


where 
m 
ee 
4ßt | (9) 
K=m| y] | 
Therefore 
B=Im[Af J, (y)eXe-P'v'dy], (ro) 
0 
where 
pe = ie’, (i= —1) (11) 


From WATSON (1944), p. 394 (5): 
co Let 2 
»(ay)e-P dy = Ne Br tee 
JJa(ape dy = ap) (2) 
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: x 
With v=-, a=1, we find: 


AV Ka i 
B=Im Fea n( - 345) |. (13) 


where fy, is the Bessel function of the first kind 
for imaginary argument. 
Since 


D 


-a<arg(-<5)<%, (14) 


we may use the formula (Watson, p. 77) 


i Br aie 
1-35) =: ul - ge | (15) 


Making use of (15) and the expression 


(16) 
and noting that 


pe (17) 


L 
where L is the wavelength of the initial 
distribution, (13) finally reduces to 


B= ~2Asin (PE) c cos (+) (18) 


The expression (18) may be used to compute 
the exact height change of a given sinusoidal 
wave for various phase angles ([7|) and 
various wavelengths L at given latitudes for a 
given forecast period. 

The next step is to evaluate the same integral 
as (5), but with the lower limit replaced by 
(- &,). We therefore consider the integral 


; see % 
Bu = f'Asnnt -(8-I(V- A + 


Fr (19) 
z AV = Bs) dé, (ET = 0), | 


= 


or 


Sal 
B, = AVBt/f sin m(| €] 
0 


(20) 
If we let 
IE] Zr (7) 
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(20) becomes 
Ba a" sn ("u mlzl).Ao)dy (22) 
With the following substitutions: 
b = 2/6: é, | | 
y=bz (22) 


and K=m|y| as in (9) above, (21) finally 


becomes 


B, = A|sin K/b cos ("= *) jb) ues ye 
— cos Kfb sin GR dz 


The integrals in (23) are expressible in terms 
of the Lommel functions of two variables, 


U, (=. ), 
Ww 


From Watson, p. 543 (1 & 2): 
Jun (bz) cos (= 5 Jde =sin (2) (Ex (=. ) = 
— cos (2) u,(= i) + u, o) 
2 w w 
22 ) b2 
fo (bz) sin (= ae (“u (re b) 
b2 


h2 af co x n+2r | 
ui b= ECG) oad je 
b: n+2r Uo 

UT. 0) = en Gn) lé 
w ei) | 


= sin (K- | LE (=. i) + sin KU, es 0) - 
b2 
— cos KU, ne | 


(26) 
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If we now substitute (22), (7), (17), (9) into 


(26), we obtain 
(= vB) + 


+sin (Al 1x) Us (= SE) . 
+ sin 2 (EE, o) = 


ll y (BL, 0) 


B, = Al cos el —|y|)U, 


= à 
cos L (27) 


If we consider what happens to B, for large 


1&,1, we find ; 
Lim B, = = A sin 2x xl u o) = | 
al L IT 


Zu 


if (Ee), | 


since the U, and U, functions go to zero for 
large | &, |. 

It may now readily be shown that the right 
side of (28) is indeed equal to the right side 


of (18). 
Thus 


(28) 


Lim B, = B 


|&;| 00 


(29) 


which, of course, is to be expected. Therefore, 
from (27), we have, 


BB - AI cos (1611 - 1x) 


Lt 
U, (= 
TT 


Finally, 


i 2 


Les 


2 FA + sin 2 ((8,|- lg) $30) 


u, (Evi) | 


| cs (lé: — Il) U (& a) + 
+sn7  l&l - x) u,( ©, vil) if 


(31) 


TLOLUITESEIBEEFRSKSOJEISIRENY: 


3. Results 


The result (31) is exact. It is unfortunate that 
available tables of Lommel functions (See E. 
LomMEL, 1886, pp. 229-239) make it im- 
possible to carry the calculations out to very 


5 10 ı5 20 25 30 35 40 


> 
(DEGREES. LONGITUDE) 


Fig. 1. Absolute value of error (E) in percent, vs. distanc 
out (| &,|) from origin. day latitudes SN, 
L = 90° longitude, | y | = o. 


large distances from the origin. Fig. 1 shows 
the graph of E vs | &, | for t=1 day, latitude 
45° N, L=90° longitude, |y|=o. From the 
graph it is apparent that one may be assured 
of an error 10 % if the numerical integration 
is carried out beyond 20° longitude from the 
ae In fact, at 35° away, the error is only 
2.6 %. The rest of the curve will oscillate with 
es fof decreasing amplitude, much like the 
Bessel functions, and finally go to zero. 

It is of interest to note that, in previous 
applications of equation (1), we have used 
APS longitude as a rough estimate of the 
cut-off point. It now turns out that this was 
a fortunate choice, particularly when the 
forecast is being made for a point of zero phase 
angle of initial conditions. 

in practice, since it is clear from (31) that 
E varies with | | as well as with ß, t, and L, 
one would construct a table of E for various 
, B, t, and L for quick reference in pre- 
paration oF forecasts. 
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Transport and Production of Vorticity in the 


Atmosphere 


By J. VAN MIEGHEM, University of Brussels 


(Manuscript received November s, 1953) 


Abstract 


The classical form of the vorticity equation is transformed into an equation of balance from 
which the vorticity transport and the rate of vorticity production are deduced. Transport and 
production are given for I) an arbitrary component of the absolute vorticity, 2) the absolute 
vorticity vector, 3) the absolute vorticity flux through a surface, 4) the total absolute vorticity 
within a volume. Special consideration is given to the absolute vertical vorticity, the absolute 
vorticity component normal to the isentropic surfaces and the total absolute vertical vorticity 
of the atmosphere. The influence of vertical motions and of heating and cooling on the vorticity 
transport is discussed. The sign of the meridional vorticity transport is tentatively determined. 


1. A convenient form of the vorticity equa- 
tion for the study of the general circulation is 
the equation of balance of the absolute vorticity 
giving explicitly the local rate of change of the 
absolute vorticity per unit volume in terms of 
a vorticity transport and a vorticity production. 

The classical form of the vorticity equation 
may be written, with the usual vector notation, 


er div v=Q-Vv,+(Vpx 


+ curl F, 


a) + 
(rer) 


where Q, represents the component of the 
absolute vorticity Q along an arbitrary axis z, 
fixed with respect to the earth, v the wind 
velocity, p the atmospheric pressure, « the 
specific volume of the air, F the frictional 


force per unit of mass, v the “del” operator 


eA | 
2 do ne 
dt dt 


along the air motion with respect to the earth. 
The suffix z designates the component of a 
vector along the z-axis. 


and individual derivative 


Substituting in (1.1) the well known iden- 
tities 

div Q = 0 and div(a A) =a div A+ A-vVa, 
where A is an arbitrary vector and a an arbi- 


trary scalar quantity, the vorticity equation 
assumes the form of an equation of balance: 


0 @ 
ot 


+ div(Q.v-v,Q) = (vp x va). + curl, F 


(1.2) 
a form which is appropriate for the evaluation 
of local vorticity changes. 

The equation of balance (1.2) is simply 
another expression for the circulation theorem 
of V. Bjerknes. Thus, the z-component Q. of 
the absolute vorticity Q in unit volume, at rest 


with respect to the earth, varies in consequence 
of: 


a) a vorticity transport 
C=Q.v-17,Q (1.3) 
across the boundary of the unit of volume; 
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b) a vorticity production within the unit 
volume at the rate 


(vp x Va). + curl, F (1.4) 


determined by the isobaric-isosteric solenoides 
and the vortex unit-tubes of the frictional 
force F. 

The term Q.v of the transport C of the 
z-component Q. of the absolute vorticity Q 
represents the convection with the wind 
velocity v of the instantaneous absolute 
rotation of the air about the z-axis, the non- 
convective term — v;Q, the contribution to 
the transport of the z-component of the 
absolute vorticity, of the instantaneous absolute 
rotation of the air moving in the 2-direction. 

The transport C, of the z-component Q; of 
the absolute vorticity Q in an arbitrary direction 
y is evidently given by 


C0, = 20 (5:34) 


where the suffix y designates the component of 
a vector in the y-direction. It should be noted 
that: 

1) the transport in the y-direction of the 
z-component of the absolute vorticity is 
numerically equal, but of opposite sign, to the 
transport in the z-direction of the y-component 
of the absolute vorticity; 

2) there is no transport in the z-direction of 
the z-component of the absolute vorticity. 

Hence, the three-dimensional divergence in 
(1.2) reduces to a two-dimensional one (see 


eq. (2.2)). 


2. Let us now suppose that the z-axis is no 
longer arbitrary but coincides with the upward 
pointing vertical z. The vertical component 
Q. of the absolute vorticity Q is the vortex 
quantity most commonly considered in dy- 
namic meteorology and is often called “the 
absolute vorticity’. 

The transport of the absolute vertical vorticity 
Q. is defined by the horizontal vector, (1.3), 


(2.1) 


where as usual the suffix h refers to a horizontal 
vector. Hence, the transport of the absolute 
vertical vorticity occurs along lines of the 
geopotential surfaces, (Cz = 0). 

Therefore the equation of balance (1.2) of 
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cC=Qv-,Q0=-Qv,-,Q,=C, 


DT 


the absolute vertical vorticity Q; assumes the 
simpler form: 


0246 
ee ati div, (Q: Vh = Vz Q;) = 


= (Vp x Va). + curl, F, (2.2) 
where div, represents the “horizontal” diver- 
gence (two-dimensional divergence in the 
geopotential surfaces). 

The convective part Q:v, of the vorticity 
transport € (or C,) results from the advection 
of the instantaneous absolute rotation of the 
air about the vertical axis z and the non- 
convective part -v.Q, from the instantaneous 
rotation of the air flowing vertically through 
the geopotential surfaces. The non convective 
transport -v.Q, represents the transformation 
or conversion of horizontal vorticity Q, into 
vertical Q.. This transport vanishes and hence, 
the absolute vertical vorticity Q: is simply 
advected in the geopotential surfaces, when the 
air moves horizontally (vz = 0). 

If N represents the external normal to the 
earth’s surface S (vn=0), the transport of 
absolute vertical vorticity across the earth’s 
surface is determined by -v;Qn. This transport 
is zero at all points where the earth’s surface 
S coincides with a geopotential surface (vz = 0), 
and therefore exists only in mountain regions 
(v: #0). The vorticity transport -v,Qx across 
the earth’s surface represents the mountain 
effect [see also eq. (3.1")| on the redistribution 
in the atmosphere of absolute vertical vorticity. 
It is readily seen from (2.1) that 


= v-QN = Q: (Vi : N) Là v(Qr à N), 


where N designates the unit-vector of the 
N-axis. Hence, the mountain-effect on the 
distribution of the absolute vertical vorticity 
consists of advection Q:(v,-N) = — Q.: 
. (v N) of absolute vertical vorticity and con- 
version -v;(Q, N) of horizontal into vertical 
absolute vorticity at the earth’s surface. With 
increasing slope of the mountains the advection 
term decreases whilst the conversion term 
becomes predominant. At the limiting case of 
a vertical wall the transport of vertical absolute 
vorticity across the wall is due only to con- 
version of horizontal into vertical vorticity. 

The rate of production of the absolute vertical 
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vorticity Q determined by the number of 
isobaric-isosteric solenoides and vortex unit- 
tubes of the frictional force per unit area in 
the geopotential surfaces, may be given the 
following expressions: 


(vp x Va)z+ curl F=curl,(F - avp) = 


= (9 x VO),+ curl, F= | 
“ove x u \ (2.3) 
= fu, ae + curl, F, | 


where © represents the potential temperature 
of the air, u, the geostrophic wind velocity, 
f the Coriolis parameter and V, the horizontal 
“del” operator. 

In the surface layer of the subtropical high 
pressure belt, curl. F>o and in the surface 
layer of the subpolar low pressure belt, curl, 
F<o. Hence, at the earth’s surface, there is a 
source region of (cyclonic) vorticity in the 
subtropical latitudes and a sink in the subpolar 
latitudes, (S. PETTERSSEN, 1950). The other 
term of the rate of production (2.3) depends 
upon the “horizontal” baroclinicity and is in 
fact rather small outside the frontal zones. 


3. Integrating the equation of balance (1.2) 
over a volume 7 at rest with respect to the 
earth (df = 0), one finds the equation of balance 
of the total absolute vertical vorticity ///Q.6t 

T 


within the volume 7, 
d 
JAN QT I (Qn 2; v;Q,)do = is (3 1) 


where n is the external normal to the boundary 
o of the volume t and where the rate of 
production & of the total absolute vertical 
vorticity within the volume r has the following 
expressions, (2.3), 


Cafff(vp x va). + curlFlör= | 
=/ff curl, (F - avp(ör = | 

T \ 
=/f [nx (F-avp)]o= | 


=/f 


[02 


[ fu, -n + (n x F).Jör, | 


n being the unit-vector of the n-axis. As the 
volume integral expressing the production & 
reduces, by virtue of Green’s theorem, to a 
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surface integral over the boundary of the 
volume, the sources of the absolute vertical 
vorticity at all the interior points of the volume 
t do not contribute to the production of the 
total absolute vertical vorticity within the 
volume r, and hence, the rate of change of the 
total absolute vertical vorticity within the 
volume t of air, at rest with respect to the 
earth, depends only on the state of atmospheric 
motion at the boundary o of the volume 7, 
(see also section 6). 

If the volume 7 is multiply connected, the 
boundary o is the “complete” boundary of r. 
Moreover we have implicitly assumed that the 
functions under the integration symbols in 
(3.1) and (3.2) satisfy the required analytical 
conditions. 

Considering the whole atmosphere instead 
of a volume 7 of air, we obtain the equation of 
balance of the total absolute vertical vorticity 
of the atmosphere 


[if de +ff IQ + | 


2 
Meee Ss 3-19) 
+[N x (F- aVp)].}öS= a 


or 


à 
3 [WQbr+ ff DeQx+ (NX Bet gun 
+ fur 0 N]öS = O, 


provided the horizontal components of the 
acting forces F and -«vp vanish at infinity 
with sufficient rapidity except when the ex- 
ternal boundary is a geopotential surface 
(see section 6). There are no reasons to believe 
that this requirement is not fulfilled in the 
atmosphere. The three terms in square brackets 
in (3.1”) are identically equal to zero at all 
points where the earth’s surface coincides with 
a geopotential surface. In the ideal case of an 
earth’s surface S coinciding everywhere with 
a geopotential surface (the mean sea-level, for 
instance), the total absolute vertical vorticity of 
the atmosphere is constant in time. This con- 
stant is zero, |see formula (3.3)]. 

Thus, the instantaneous rate of change of the 
total absolute vertical vorticity of the atmos- 
phere results from a mountain effect only. As 
stationary conditions prevail in the long run, 
the time average of the surface integral in 
(3.1”) for a sufficiently long time interval must 
be equal to zero. 
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The influence of the earth’s surface S on the 
total absolute vertical vorticity of the atmos- 
phere may be shown in another way. Re- 
calling that the absolute vorticity Q is equal 
to the sum of the earth vorticity 2@ and the 
telative vorticity curl v of the air, one finds 


JT Qzöt = [ff 2@ sin y dt + 


atm. atm. 


+ [ff curl, v ör, 


atm. 
or, applying Green’s theorem, 


JT Q6t = ff] 2m sin g or + ff (v x N).öS (3.3) 


atm atm. S 


where @ is the constant angular speed of the 
earth and the latitude. If the space occupied 
by the atmosphere is symmetrical with respect 
:o the earth’s axis and the plane of the equator, 
the volume integral in the second member of 
(3.3) is identically equal to zero. This double 
symmetry not being exactly realized, the value 


of //f 2m sin gör depends upon the shape of 
atm. 
the earth’s surface and is therefore a constant. 


It should be noted that 
SI 20 sinpör=/[f wrcos p(i x N)20S, 
S 


atm. 


where r is the distance from the earth’s centre 
of an arbitrary point of the earth’s surface S 
and i the unit-vector tangent to the latitude 
circle and pointing from west to east at the 
point considered. 

The surface integral in (3.3) is identically 
equal to zero if the earth’s surface is a geo- 
potential surface. This again is not truly so, 
and the total absolute vertical vorticity there- 
fore depends upon the slope of the earth’s 
surface and of the air motion at the earth’s 
surface in those regions where the earth’s 
surface differs from a geopotential surface 
(mountain regions). 


4. Integrating now the equation of balance 
(2.2) over a horizontal area 2 (an area of a 
given geopotential surface), at rest with respect 
to the earth, we obtain 


ker) +f (Q.v, - vzQ)ös= [(Fsös - xöp), 


s 
(4-1) 
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where » is the external normal to the boundary 
s of Xin the geopotential surface considered. 
The integration sense along the closed curve s 
is counterclockwise (right-handed screw rule). 
It is well known that / -xöp represents the 


Ss 
number of isobaric-isosteric solenoides em- 
braced by the horizontal closed curve s and 
that, in the northern hemisphere, [ F,ös is 


Si 
positive [production of vertical (cyclonic) 
vorticity] when the wind circulation is anti- 
cyclonic and, negative [destruction of vertical 
(cyclonic) vorticity] when the wind circulation 
is cyclonic along s. When stationary conditions 
are realized, a production of absolute vorticity 


flux //Q.6X through the horizontal area & 
Bs 


(or what amounts to the same thing, a produc- 
tion of absolute velocity circulation along the 
horizontal closed curve s, boundary of 2) 
must be compensated by a horizontal vorticity 
outflow /(Qzv,-v;Q,)ös>o across s and 


Ss 
a destruction of vorticity flux through 2, by 
a horizontal vorticity inflow across s. 


5. Assuming the geostrophic hypothesis and 
using the thermal wind equation, the absolute 
horizontal vorticity Q, may be given the 
following expression, valid as a first approxi- 
mation, 


g v9 | 
jad 


(5.1) 


Qi, = 209, 


where g is the acceleration of gravity. Substi- 
tuting (5.1) in (2.1) we find a useful approxi- 
mate expression for the transport of absolute 
vertical vorticity, namely: 

VO 


C=C,~Q.v, +7 ee (5.2) 


In an air current vy, vertical vorticity is 
advected with the air, but at the same time, 
vertical vorticity is also transported approxi- 
mately in the direction of the horizontal 
temperature ascendent (roughly to the south, 
on the scale of the general circulation) where 
the air ascents (vz > 0) or in the direction of the 
horizontal temperature gradient (roughly to 
the north, on the scale of the general circula- 
tion) where the air descends (vz <o). Hence, 
when in a horizontal layer the horizontal 
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temperature gradient (- v„©) ~ converges 
(cold pool), vertical (cyclonic) vorticity may 
accumulate as a consequence of convergence in 
the horizontal transport of vorticity, if the cold 
air descends. Similarly, when in a horizontal 
layer the horizontal temperature ascendent 
„9 converges (warm pool), vertical (cy- 
clonic) vorticity may accumulate as a con- 
sequence of convergence in the horizontal 
transport of vorticity, if the warm air ascents. 


6. More generally, the rate of production of the 
absolute vorticity Q per unit of volume is 


defined by the vector 


vp x Va+curl F, 


and the transport by the antisymmetrical tensor 
Qv Ex vQ, 


the component of which along the n-axis, 
Qv„— vQ,, determines the transport of the 
absolute vorticity Q in the direction n, (VAN 
MIEGHEM, 1951). 

Consequently, the vorticity budget of a 
given volume + of air, at rest with respect to 
the earth (öf=o), assumes the form 


2 Lf Qör + Lf (Qua -vn)ör=2, (6.1) 


where the production 2 of the total absolute 
vorticity //f Qör within the volume + may 
T 


assume the following expressions: 
2=//f[vpx va + curl Flôr = 
= {ff curl (F- «vp)ét = 


=/f nx (F-avp)do, 


use being made of Green’s theorem. The 
equation of balance (6.1)-with (6.2) -is 
similar to a vector equation obtained by C. 
TRUESDELL (1948). One obtains Truesdell’s 
equation by replacing the wind velocity v in 
(6.1) by the absolute wind velocity, the vector 
F — « vp by the absolute acceleration and the 


© z, 
rate of change > with respect to the earth by 
© 
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the rate of change with respect to the absolute 
frame. Truesdell’s conclusions may very easily 
be extended to the case of relative motion 
with respect to the earth, in particular the 
conclusion recalled in section 3. 

It should be noted that the gravitational 
force does not appear in (6.2), this force being 
expressed by the gradient of a potential func- 
tion (geopotential). 

Finally, considering the whole atmosphere 
instead of an arbitrary volume + of air, one 
finds the vorticity budget of the atmosphere 


2 IQ +f{[v;Qn+N x (F- a\/p)]6S = 0, 
atm. Ss 6.7’) 


where vs designates the wind velocity v at 
the earth’s surface S. Hence, the total absolute 
vorticity of the atmosphere is a function of the 
state of air motion at the earth’s surface only. 
It should be noted that the vector vs is tangent 
to S, (vn=0). 

The equation of balance (6.1’) is valid 
provided that one assumes that the normal 
component Qv, — vQ, of the tensor Qv - 
— vQ and the tangential component of the 
vector F — «vp upon a sphere concentric 
with the earth vanish more rapidely than the 
reciprocal of the square of the radius of the 
sphere, when the radius tends to infinity. If 
this requirement is not fulfilled, there exists a 
contribution from infinity to the vorticity 
budget of the atmosphere. In order to avoid 
this unrealistic circumstance we have assumed 
implicitly that the above condition is satisfied 
ipso facto, [see also eq. (3.1’)]. 


7. Using general coordinates in a frame 
moving arbitrarily with respect to the absolute 
frame (VAN MIEGHEM, 1951), it is possible to 
generalize the equation (6.1) to the case of a 
volume t moving with the new relative frame 
(0't=0). The same equation of balance is valid 
provided that one replaces the wind velocity 
v by the air velocity w with respect to the 


AV 


à : 1 à 
new frame and er by the time derivative — in 
[4 [4 


the same new frame. One obtains 
9’ 
al Qör + J(Qw, - wQ,)5'o = | 


(7-1) 
=/[nx (F-«vp)ö'o. | 
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Similarly, the equation of balance (4.1) may 
be extended. For instance, supposing that one 
of the new coordinates is the potential tempera- 
ture ©, we find the equation of balance 


a 
al Q,,0'2' + f (Quuw, — w,Q,) O's =f F,6's, 
(7.2) 


where » is an area of an isentropic surface 
(© =constant), s the complete boundary of 7, 
n the upward pointing normal to the isentropic 
surface, » the external normal to the curve s 
À : : .vO 

in the isentropic surface, Q, = Fe the 

Vv 

absolute vorticity component normal to the 
isentropic surfaces and w the air velocity in 
a relative frame moving with these surtaces 


(w, = O being the individual time 


O 
Fo]: 
derivative of ©). 

In (7.2) the term - x dp has dropped out, 
x Op being an exact differential on an isentropic 
surface. 

The transport 

r= Q, W —W,Q (73) 
of the absolute vorticity component normal 
to an isentropic surface is evidently isentropic 
(I„=0). 

Employing longitude, latitude, potential tem- 
perature, and time as independent eulerian 
variables, it may be shown that the zonal (x), 
meridional (y) and vertical (z) components of 
the vorticity flux I assume respectively the 
forms: 


if = Qn x z: WnQx| 
fie: = On, is wnQ) 


and 


Fr; = O7 Vh = WnQh; (Ge 


= (VO -T;,) : = 6); 


where v,, vy and Q,, Q, are the usual zonal 
and meridional components of the wind 
velocity v and the absolute vorticity Q and 


where ©, stands for —. 
dz 


The slope of an isentropic surface being small 
(10-8 to 107%), the isentropic transport T of 
the absolute vorticity component normal to 
the isentropic surfaces may be identified, for 
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all practical purposes, with the horizonta 
component I, of T. 

If sis a closed curve in an isentropic surface, 
along which the transport of Q, takes place, 
T',=0; hence the instantaneous rate of change 


of the absolute vorticity flux /f Q,62 through 
Z 


an isentropic area 2 (or what amounts to the 
same thing of the absolute velocity circulation 
along the isentropic closed curve s) is due to 
a frictional effect only. In an inviscid fluid, 
the absolute vorticity flux through an isen- 
tropic area is independent of time f provided 
that there is no isentropic absolute vorticity 
transport across the complete boundary of the 
isentropic area. When moreover the processes 
are adiabatic (w,=0), the absolute vorticity 
flux through an isentropic area is independent 
of time t when the boundary of the area is a 
streamline of the isentropic flow. In this 
simplified case, the Q,-vorticity is simply 
advected in the isentropic surfaces along the 
streamlines of the isentropic flow. 

Let us now consider an isentropic surface 
intersecting the mean sea level along a curve 
s entirely situated in the region of the westerlies, 
then /F,6’s<o, and if stationary conditions 


Ss 
prevail, the vorticity destruction along s must 
be compensated by an isentropic vorticity 
inflow across s, corresponding roughly to a 
horizontal vorticity transport to the north. If 
on the other hand s is entirely situated in the 
region of the tropical easterlies, [F,ö’s > 0, 


Ss 
and if stationary conditions again prevail, the 
vorticity production along s must be com- 
pensated by an isentropic vorticity outflow 
across s, corresponding roughly to a horizontal 
vorticity transport to the south. Hence the 
subtropical high pressure belt at the earth’s 
surface, which is a source region of (cyclonic) 
vorticity (see section 2), is also a region of 
horizontal divergence of (cyclonic) vorticity 
transport, stationary conditions being realized. 
Similarly the subpolar low pressure belt at the 
earth’s surface, which is a sink region of 
(cyclonic) vorticity, is also a region of hori- 
zontal convergence of (cyclonic) vorticity 
transport, stationary conditions being realized, 
(E. T. Eapy, 1950 and S. PETTERSSEN, 1950). 


8. From the point of view of the general 
circulation, it is of interest to consider the 
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transport C, of the absolute vertical vorticity 
Q. along the meridians (positive y-axis to the 
north), 


(Cr = Onn = 0.07 ORs 6% 5 (8.1) 


Let us first examine the advective meridional 
transport Q; v, in middle latitudes. Recalling 
that Q, is generally positive in the northern 
hemisphere, this transport is in a northerly 
direction to the east (v, > o) and southerly to 
the west (v,<o) of the upper air troughs in 
the westerlies. On the average along a latitude 
circle the advective meridional transport of 
vertical vorticity will be towards the pole if 
the vertical vorticity is larger to the east than 
to the west of the troughs, but towards the 
equator if, on the contrary, the vertical 
vorticity is larger to the west than to the east 
of the troughs in the westerlies. Apparently 
the first alternative occurs most frequently in 
the troposphere. 

It is much more difficult to state precisely 
the sign of the non advective mean meridional 
transport along a latitude circle. As a rule 


20 
ao, and presumably, on the average, v; > 0 
[4 


to the east of the troughs in the westerlies and 
v;<o to the west. Hence, the advective and 
the non-advective transports have generally 
opposite signs, but the non-advective transport 
is at least of one order of magnitude (very 
often even two orders of magnitude) less than 
the advective transport. 

The non-advective transport ee is ın a 
southerly direction to the east and northerly 
to the west of the upper air troughs. On the 
average, the non-advective meridional transport 
of absolute vertical vorticity across a latitude 
circle will be towards the pole when the 
meridional poleward gradient of potential 
temperature is larger to the west than to the 
east of the upper air troughs, but towards the 
equator if, on the contrary, the meridional 
poleward gradient of potential temperature is 
larger to the east than to the west of thetroughs. 
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Generally, in deepening upper air troughs, 
the baroclinicity is larger to the east than to 
the west of the troughs and hence, the zonal 
mean non-advective transport of vorticity 
across a latitude circle is presumably towards 
the equator. 


9. In the same way that formula (2.1) is 
approximately given by formula (5.2), we may 
also approximate the formula (7.3), 


I, = CO, Vn +? ve (9.1) 


where Q, ~Q., except in regions of high 
baroclinicity, and obtain for the non-advective 
term in (9.1) results similar to those of section 
5 regarding the non-advective term of (5.2). 

In an air current v,, absolute vorticity 
normal to the isentropic surfaces is advected 
with the air, but, at the same time, also 
transported approximately in the direction of 
the horizontal temperature ascendent (roughly 
to the south, on the scale of the general 
circulation) where the air is heated (w, > 0) or 
in the direction of the horizontal temperature 
gradient (roughly to the north, on the scale of 
the general circulation) where the air is cooled 
(w, < 0). Hence, when in a horizontal layer the 
horizontal temperature gradient (— V,@) con- 
verges (cold pool), (cyclonic) vorticity normal 
to the isentropic surfaces may accumulate as a 
consequence of convergence in the horizontal 
vorticity transport, if the cold air cools down. 
Similarly, when in a horizontal layer the 
horizontal temperature gradient diverges 
(warm pool), (cyclonic) vorticity normal to 
the isentropic surfaces may accumulate as a 
consequence of convergence in the horizontal 
vorticity transport, if the warm air warms up. 

Finally, it should be noted that the velocity 
components v, (meridional wind component) 
and w, (heating or cooling) have generally 
Opposite signs in the upper air troughs and 
ridges. Hence the non-advective meridional 
transport of absolute vorticity normal to the 
isentropic surfaces, due to heating and cooling, 
reinforces, as a rule, the advective meridional 
vorticity transport. 
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Introduction 


Several recent investigations, Kuo (1951), 
VAN MEIGHEM (1953), STARR (1953) and 
ARAKAWA (1953), have been concerned with 
examining the physical processes which account 
for the maintenance of the kinetic energy of 
zonal motion in the atmosphere. These studies 
have pointed up the important role played by 
the large scale horizontal eddies. The funda- 
mental point made by all these investigators 
centers on the existence within the free 
atmosphere of a systematic organization of the 
large scale eddy momentum exchange pro- 
cesses in such a manner that momentum flows 
from regions of low to regions of high mo- 
mentum. The reflection of this condition is 
most prominent equatorward of the jet stream 
where for example observational studies of 
Starr and WHITE (1952) show that the flux 
of momentum by the large scale horizontal 
eddy processes is poleward into the jet. STARR 
(1953) and ArakawA (1953) have actually 
evaluated the work done by such large scale 
eddy stresses and indicate that the rate of 
generation of kinetic energy of zonal motion 
by such processes seems to be sufficient to 
maintain the westerlies of middle and high 
latitudes. It is of interest, therefore, to inquire 
whether similar conditions prevail in the 
atmosphere with regard to the transfer of other 
atmospheric properties and in particular with 
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respect to the transfer of sensible heat. Sufficient 
data are now at hand to show that similar 
counter-gradient sensible heat flows do occur 
on a systematic basis. 


Data 


In connection with a more extensive study 
of the atmospheric energy balance, computa- 
tions of the large scale horizontal eddy flux of 
sensible heat were undertaken. Except for a 
constant factor, the large scale meridional 
horizontal eddy flux of sensible heat is measured 


by the term [v’ T’], where the brackets re- 
present a mean with respect to longitude, the 
bar a mean with respect to time, the primes 
represent deviations from the longitudinal 
means, v is the northward component of the 
wind velocity, and T the absolute temperature. 
The evaluation of this term at various latitudes 
and at the standard pressure levels in the 
atmosphere, was made using observed wind 
and radiosonde observations available in the 
data compilations of the Northern Hemisphere 
Historical Weather Map Series, prepared by 
the U.S. Weather Bureau, in cooperation 
with the Army, Navy and Air Force, for the 
full year 1950. The procedure was similar in 
many respects to that used in connection with 
investigations of the angular momentum 
balance reported by Starr and WHITE (1951), 
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Table 1. Numerical analysis of sensible heat flux data for the entire year 1950. Temperature in degrees 

absolute, velocities in m s2c-1 and N is the total number of observations. The 95 % confidence 
intervals defined as twice the standard error of the time means is also given. 


07 ° EN ıN 
Tr 707 N sen 42.5 3 
ao) I = — — = — = 
[7] A a | N | era | N 7 ae | NT ER) | ae N 
100 | 225} +0.4+2.9 | 288] 220] +5.6+2.7 | 763] 215] +4.8+3.8 | 546] 206] +3.2 +2.9 |1253 
200 | 223] +1.4+4.0 | 919] 220] +5.2+4.0 |1601] 2190] +8.3+4.5 |1684| 218) +7.442.7 |2859 
300 | 222| +0.0+3.2 |1302| 225} +2.5 +2.7 |2361| 231 +3.0+3.3 [2463| 236| +3.0+1.8 [4003 
500 | 244| +5.2+3.1 |2047| 249] +8.5 +2.4 13295] 256] +1.9+1.9 [3625| 262| +1.5+1.1 [5181 
700 | 259| +9.2 42.1 |2266| 264| +12.7 +2.0 |3599| 272] +6.6+1.3 |4319| 278| +2.3 +0.9 |5756 
850 | 266| +9.6+2.0 [2374| 271l +17.5 +2.0 |3684| 279| +10.8 +1.5 |4425| 285] +5.0+1.1 |5218 
1000 | 267| +7.6+2.3 |1716| 275| +11.8 +2.0 |2207| 284| +5.5 +1.2 !2494| 292| +2.840.7 [4507 


in which strings of stations in the vicinity of 
certain latitudes were selected from which 
observations were collected. A complete listing 
of these stations together with the frequency 
distribution of the observations will be forth- 
coming in a paper soon to be published by 
Starr and White. 


For purposes at hand it is sufficient to present 
the relevant data in Table 1, where the mean 
temperature, magnitude of the eddy transfer 
of sensible heat, and total number of observa- 
tions entering the evaluations, are shown. 
From this table it will be noticed that through- 
out the troposphere, at all latitudes up to the 
200 mb level the eddy-flux of sensible heat is 
poleward from regions of high to regions of 
low temperature as might normally be ex- 
pected. At and above this level the reverse is 
true. There is an eddy flux of sensible heat from 
the cold tropical lower stratosphere poleward 
to the warm polar regions at these levels. This 
condition was also recognized by PRIESTLEY 
(1949) and Mintz (1951) from data compila- 
tions which were restricted in space or time. 
The 95 % confidence limits of the time means 
of the eddy-transfer term are also indicated and 
in almost all cases exclude zero, indicating that 
at least a necessary statistical condition is 


satisfied by these data. 


Meteorological implications 


a. The existence of conditions illustrated by 
these data, and similar ones indicated by 
STARR (1953) point up the fact that meteorolo- 
gists should recognize not only the possible 
existence of these counter-gradient eddy 
transfer processes but also the actual prevalence 
of such conditions in certain regions of the 
atmosphere. 

b. The flux of sensible heat from cold to 
warm regions on the scale of the general 
circulation at and above the jet-stream- 
tropopause level indicates that the eddy 
processes are acting to build up rather than 
dissipate the existing temperature gradient. 

c. Certain fundamental questions concerning 
the location of the height of the tropopause, 
and the height of the jet are inseparable from 
the questions concerning the reversed merid- 
ional temperature gradient in the lower strato- 
sphere. Hence a process such as that indicated 
by these data which obviously acts to build up 
this reversed gradient may be of considerable 
importance. 

d. Since the full significance of these ob- 
servational findings is not clear, the data are 
presented primarily in the hope that they will 
encourage other interested investigators to 
consider the implications of such conditions. 
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The aiın of this note is to present more 
statistics concerning the meridional flux of 
angular momentum in the tropics similar to 
those already published (STArRR and WHITE 
1952 a). A further investigation for a period of 
one additional ycar has now been completed 
for 13° N, in order to verify the previous 
conclusions with independent data. A continu- 
ous two-year period has thus been analysed and 
represents a concluded phase of the investiga- 
tion for this latitude. 


Except as indicated, the form of the accom- 
panying tables and the computational pro- 
cedures used in obtaining them are the same 
as the corresponding ones in the paper cited, 
to which the reader is referred for details not 
repeated here. As an experiment the station 
network for the second year was altered con- 
siderably and enlarged, the added new key 
stations being listed in Table 1. Some former 
key stations were also dropped, the complete 
new list being that given in Table 2, which 


Table 1. List of added key stations (numbered) and their alternate stations (italics) 


Station Latitude (N.) Longitude Altitude (ft) Type 

ig (etre ISIE 5 35 Spo ao ado soe TE On OP AB! 18 262 pilot balloon 
SCUAOWOVe 2.6,00 Ho eterna: Ke) PRS) 94 21 EB 30 » 
BASSE eee wenn erene set 16° #46 94 46 E 13 » 
Mingaladon................ 16 54 OG wie IB 92 » 
TALONS ee coe: I4 06 98713 2 1012 » 
MICH EEG 2 BE 98 36 E 66 » 
Victoria ONE ER 09 58 SERRE 122 » 
TRES 6.0.6.9 2.0.0.0 0 on 08a one 05 18 Too 16 E 112 » 

OO Re Dacia Dites de I9 44 155 04 W 33 radio wind 

SH Vera Cruzin ro iste 19 12 96 08 W 10 pilot balloon 
TRACUOGV Gn ee uke) 27 99 12 W 7579 » 
IMICSLCO RG ULV yet Go 19 26 99 08 W 7340 » 
Cuidad Del Carmer......... 18 39 OI 49 W 5 » 
OUAIS ANT RTC 06 ED A 83 56 W 35 » 

4 Al brooks Ete dire EEE PRE 08 58 79 33 W Dit radio wind 
Plato Magdalend... see 09 48 74 48 W 190 » 

5. Cayenne-Rochambeau....... 04 50 Be Be M pilot balloon 
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Table 2. Percentage of total possible observations for second year 
u a a Te a DE IT BE TE a I VPE ee EEE 


Station 


Altitude in thousands of feet 


2 | 6 | Io | 14 | 20 | 25 | 30 | 35 | 40 | 45 | 50 | 55 
CAO PRE er Poe 64 65 62 44 II fo) o fo) o fo) fe) fe) 
BS HaSher ita cals oem 92 93 2 OI 84 33 27 I4 10 4 I I 
Aden. Sy SVE OSHS. EEE 75 85 73 13 7 21 64 10 58 3 3 AI 
PRTRCHINOPOLY sare ace mens 100 O | 100 fe) 88 84 75 Gi fo) fe) o o 
Port Blau Idee 97 86 OI 41 52 24 13 2 fo) fo) o fo) 
SAO TM tre cle RE Eat es 89 86 88 76 63 13 5 3 2 I fe) o 
lark: I SR RER 97 97 96 92 85 53 46 34 38 29 28 30 
PY ate eee eme me ne 75 84 78 58 64 43 46 30 32 20 Io 8 
Harmon ucld ER se eae 88 94 93 78 89 72 86 Gt 83 65 54. 48 
wa aline ee 99 | 100 99 98 99 95 98 90 96 80 62 57 
Johnson Island 2.7. 75 OI 90 70 79 57 63 45 49. 31 27 16 
EIS Oats Se eke esis sree ne 90 99 99 90 97 84 OI 76 82 57 Bo 30 
NOTAMCEUZ ER tee Seas 98 94 89 79 62 Io 2 fo) fo) fe) fe) fo} 
ENA TOO ke Et CLC ara 02 ce ce 82 93 92 Gr 83 61 65 47 49 24 6 it 
Wiallere Riel dsm = ass 100 | 100 | 100 96 97 84 OI TF 84 64 47 48 
Cayenne-Rochambeau..... 56 52 39 30 23 18 16 17 14 13 7 4 
Dakar RS SR 2 2 OI 78 41 35 18 28 16 2 12 8 8 


shows the frequency distribution of the availa- 
ble observations for the second year alone. Due 
to these changes in the stations used, it is not 
convenient to present a frequency table for the 
two years combined, as was done by the 
writers in the case of similar studies at 31° N 
(Starr and WHITE 1952 b). 


It is seen from Table 3 that the main results 
for the two year period are much the same as 
for the first year alone published previously. 
Again using parentheses to indicate vertical 
averaging with respect to pressure for the 


layer considered, ([u]) turns out to be — 0.35 
m sec”! for the two years as compared with 


- 0.50 m sec”! for the first year!. The quantity 


([v]) is -o.21 m sec"! for the two years as 
compared with -0.16 m sec! for the first 
year, indicating most probably that a slight bias 
in favor of northerly wind components is still 
present in the data. The integrals at the foot 
of columns 5 and 6 in Table 3 show that about 
II per cent of the momentum transport is 
accomplished through mean meridional cir- 
culations for the two years as compared with 
about 14 per cent for the first year. 


1 This last quantity was written erroneously as + 0.50 
m sec-! in STARR and WHITE (1952 a), although the 
associated discussion there given requires no correction. 
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Numerical Solutions of the Perturbation Equation for Linear Flow 


By ERIK ELIASEN, University of Copenhagen 


(Manuscript received August 20, 1953) 


Abstract 


Two-dimensional nondivergent perturbations of linear flow are considered. By means 
of a simple numerical treatment indications concerning the conditions of exponentilal instability 
are obtained. Numerical solutions are further used to represent the approximate development 
of initial perturbations up to a certain time. The correctness of this procedure is shown in the 
case of the Couette flow. Finally perturbations of a symmetrical harmonic velocity profile 
are considered in relation to the question about the instability in a barotropic atmosphere. 


1. The instability problem 


For a two-dimensional motion of a homo- 
geneous incompressible and inviscid fluid the 
vorticity is conserved. A flow in the xy-plane 
with the velocity components 


U= Uy) V=o 
is then stationary, and for a small perturbation 
of this flow expressed by a streamfunction 
D (x, y, t) we have the linearized equation 

I72D av7D DEU 
+ - (1) 
dt ox dx dy? 


As U only depends on y, we can develop ® in 
basic solutions of the form /(y)e*-“), and 
for (y) we get the following differential 
equation 

(U-c)p"-(U-c+U'ip=o (2) 
where ’ denotes a differentiation with regard 
to y. Together with the necessary boundary 


conditions this equation determines p(y). Here 
we shall consider the simple case 


p(o)=p(D)=0 (3) 
corresponding to rigid boundaries at y =o and 


y=D. 
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We have now given an eigen-value problem, 
as we shall find for which connected values 
of u and c there exist solutions of (2) and (3). 
If there exist solutions for complex values o 
c, this means that basic solutions will exist, 
whose amplitude will increase exponentially 
with time, and the considered linear flow will 
be unstable. The mathematical treatment of the 
problem is rather complicated, and generally 
we cannot find, even for simple functions U(y), 
analytical expressions for the solutions p(y). It 
therefore seems natural to treat the problem 
by a numerical method. 


In doing so we shall consider the function 
p(y) as given by its values @, in N points in 
the interval o < y < D in such a way that the 
interval is divided by these points in N+1 


D 
parts, each of the length d= The 


derivatives shall then be approximated by the 
corresponding quotients of finite differences. 
We shall here only use the simple case N=2 


D 2D 
with the values y, for y, =— and y, for y= — 


As p(O)=p(D)=0 we then have the following 
expressions for gy’ at yı and yp. 
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wm 


I Large 
PIE (92 - 29) P2 = 7 (P1 — 22) 
Accordingly we can write the equation (2) at 
these two points 


+ {(2 + u2d?)(U,—c) + AU} pa + (U, — pa = 0 
(U, 200, - (ana) (Ue) + EU, }p,= 0 
(4) 


where U,, U,, Uj’ and Uj’ are the values of U 
and U” at y, and y, respectively. In order that 
these two equations shall give non-zero solu- 
tions for @,, 9, the determinant must vanish, 
which gives 


= Ze Us” Pur up + — 


2 25 25 


- V{s(Uy - U,) + rd? (Uy! - Uy')}? + ABUS US 
(5) 
with 


2 


f= 2 bed, Ge = i 


From this it is seen that a necessary condition 
that c can get a complex value is 


Ue eae (6) 


in accordance with the well-known condition 
of exponentiel instability for a linear flow 
without friction, namely that the velocity 
profile shall have a point of inflection. That 
this point of inflection shall lie between the 
two points y, and y, is a result of the choice 
IN| S23), 

The condition (6) is, however, not a suffi- 
cient one. This can be demonstrated by con- 
sidering an anti-symmetrical profile 


Uy) = + U(D + y) 


for which (6) is satisfied. Equation (5) is then 
reduced to 


ge = V(sU, +rŒU; PU (7) 


Replacing U," with the quotient of finite 


differences à (Uy + 3U,), where U, = U(0), we 


get that c becomes complex if 
|U,u?d? + (U, + U,)u2d® + 2U, - 3U,| < | 


- (8) 
<|U= 3U,| J 
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With N=2 we can only describe a very 
smooth profile, and we assume that U, and U, 
are both positive. First we see that (8) can only 


be fulfilled if 
> CEA, (9) 


This is in accordance with the further condition 
of instability, which states that the numerical 
value of the vorticity | U’| shall have a maxi- 
mum at the point of inflection, cf. FJGRTOFT 
(1950). But even this is not a sufficient criterion. 
With 3U,>U,, (8) gives that c is complex 
when 


ped? <2 - D: (10) 
but this can only be the case if 


207 10, (11) 


As a particular example we consider the 
anti-symmetrical, harmonic profile 


Um Asin Gault Hz D (12) 


This profile has a point of inflection at y= 


> 


vb © 


where the numerical value of the vorticity has 
a maximum. (6) and (9) are then satisfied, but 
the third condition (11) is only fulfilled when 
Je! 
D> . (13) 
i.e. when the distance between the rigid 
boundaries is larger than half the profile 
wavelength. This result is proved by Homann 
(1952) for perturbations of infinite wavelength. 
The above simple considerations indicate that 
a linear flow with a velocity profile given by 
(12) and (13) is unstable for perturbations with 
wavelength L larger than a certain value, in 
our crude approach given by (10), ie 


Fig. 1 shows the imaginary part Ci, Of cas 
function of L for two such profiles with 


respectively H -> D and H=D. 
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— = 4 à Typ 


Fig. 1. The curves show the imaginary part of wc obtained 
by the approximation N = 2 for the anti-symmetrical 


0 2 ‘ é 


harmonic profile with respectively H = 3 Dand H = D. 
2 


2. The initiail-value problem 


From a total set of solutions of the eigen- 
value problem we can compose the general 
solution ®(x, y, t) of the perturbation equation 
(x), and this general solution can be adapted 
to a given initial perturbation. 

By the numerical method used in the fore- 
going approximate values of (y) can be 
obtained at some points and then by inter- 
polation for all values of y in the considered 
interval. In connexion with the initial-value 
problem, however, it seems more appropriate 
to find the approximate solutions p*(y) by 
a somewhat different method, based on a 
development of g*(y) in a finite sine series 


N 
œ*(y) = > a Sin ZT. n=1,2,3,---,N 
D 
n=1 
(14) 


The boundary conditions (3) are then satisfied 
and (2) can be written 


mi Se 
TT 11 
c= U) (wnt Ze) - U fe 
n=1 N (15) 
+ dy Sin si =O | 
We now make: the approximations 
= 
U sin = > Ems | 
HE { (16) 
U” sin = Cum Sin 
m=1 


Tellus VI (1954), 2 


185 
where 
D 
En fa 
D 
Gam == We sin si as dy 


Equation (15) then gives N linear, homogeneous 
equations in the N coefficients ap. 


N 


2 
(ern) a - > . 


| 
| 
m=1 ANT) 
| 
| 


| en E G [eee 

1 u rm D? mn +, = (En 

The condition that this system of equations 
shall possess solutions different from the zero- 
solution, is that its determinant is zero, which 
yields a relation between c and u, in form of 
an algebraic equation. For a fixed value of u 
we get generally N values cl"), real or complex. 
For each cl) we can find, apart from an arbi- 


trary factor, a set a") and the corresponding 
numerical solution 9*®. The N functions 
p*") (y) can be composed to a general expres- 
sion for the perturbation function given by 


N | 
DX (x, 7, => phn) eit el ) (18) 
n=1 


The N arbitrary constants entering this ex- 
pression can be fixed by adapting it to a given 
initial perturbation of the form 


N 
. NX) 
ES B, sin —-. 
fe 1D) 
n=1 


Since the perturbation function must be real, 
we must take either the real or the imaginary 
part of (18). 

It can now be expected, that the expression 
(18) will represent with some approximation 
the development of the initial perturbation up 
to a certain time. Denoting the error by 
p(x, y, t) we have by definition 


D (x, y, t) = O* (x, y, t) +P (x, y, t) (19) 
Writing formally 
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BEN ~ 5 | 
2 
ee = > sin | 
Far (20) 
20, ODE nr} 
U a oy U > Tn (x, t) sin D | 
n=1 
D* has been determined in such a way that 
ta) = Gale, 1 tN, 


By substitution from (19) and (20) in (1) we 
obtain 


OVE ere OT ee, 
er + a a U"’ + R(x, y, t) = 0 (21) 
where 
R (x, y, t) = fa (ef) sin 
n=N+1 


As W=o at t=0, Ÿ will at first grow only in 
consequence of the term R(x, y, f). The smaller 
R is, the longer time it lasts before the error 
becomes essential. 
We shall illustrate this method by con- 
sidering the simple Couette flow 
A 


Des 


D? osysD (22) 


for which the exact solution of the perturba- 
tion equation (1) can be obtained directly, as 
discussed by Fyorrorr (1951). Equation (2) 
can be written in this case 


(U-¢) (p" — wy) =0 (23) 


and the numerical treatment based on (14) 


yields 


N 
2 2 
c G + Ba) An — > (1 + m? 2.) . | en 


m=1 


* Funndm = © 
with 
A 
Eyes 
nn 2 
| O n+mevenn#m 
Bun Zr 8 
nm 
el 7 me a n+m odd 
707i (4? m) 


Choosing N=4 we obtain 


c1)=0.19 A, c® = 0.40 A, | 2 
c9) = 0.60 À, c4=0.81 À j 
and the corresponding solutions 
BT Non 
gr) = a) (in + 1.58 ky» an 
30 . AY 
+ 1.50 k, 3 sin a + 1.04 kı,, sin #2) 
ay amy _ 
pr = ql?) (sin Dt? 64 kı,, sin Fa 
er . 4ny 
— 0.62 k, sin a - 1.08 ky 4 sin 2) 
TE or 
pr) = al) (sin — 0.64 ky, 2 sin - 
AT 
~ 0,62 ky, sin an 08 k, asin 2) 
pr 4) = gif) (sin? - 1.58 k, 9 sin a. ~ 
+ 1.50 ky, sin Te 1.04 ky, 4 sin 42) 
(26) 
where 
TÈ 
HAE D? 
ye ag 
It 
EE RES 
Taking u =5 we obtain functions g*)(y) and 


g*)(y) as illustrated in fig. 2. They approxi- 
mate fe corresponding exact solutions of (23) 


al 
1.5} 


0 O2 0.4 0.6 0.8 1.0 y/D 

Fig. 2. The numerical solutions @*(n)(y) with N = 4 for 

the Couette flow corresponding to c(1) = 0.19A and 
c(2) = 0.40A. 
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Table ı 
m ne en sr Es u Tr u En ET u A Sn ee 
= a je * (1 2 1 2 1 2 
o 1.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00 0.00 
Ar 0.84 0.52 — 0.16 0.02 0.85 0.50 — 0.17 0.00 — 0.02 0.00 
zZ 0.47 0.71 — 0.50 — 0.14 0.50 0.68 — 0.50 — 0.19 0.00 — 0.03 
SUR 0.13 0.50 — 0.69 0,52 0.17 0.50 — 0.65 — 0.50 0.20 0.00 
4T — 0.03 0.14 — 0.50 — 0.91 0.00 0.19 -0.50 | -0.65 0.50 0.21 
N : Pr. D IV?@ Ov? 
with singularities at respectively y,=— and ee de 
S x 


ne, where U(y,) =c and U(y:) =c®). 


By means of the functions (26) we can 
compose the general perturbation function 
(18). If we take as the initial perturbation 


Considering again the initial perturbation (27) 
we have for V?®{x, y, t) the exact solution 


v2@= - a B sin {u(x — UN} sin 22 
D2 D 


@ (x, y, 0) =B sin ux- sine (27) (29) 
D This expression can be written as the following 
eR AS series 
DR, y; 0) = B} Qt sin (s- 24) sin + nr? = f À | 
YE ps oe Mee ee > : ae 5 
v2 (1 _ =) B E wn(» z 9} 


a. 0 
x a al 
+ ki 2Q3 cos a (x =) sin + 


à ANS Any 
* ae Bees 
+ ky, 3Q¥ sin u (= 5 ) sin 7 + 


AN. 
ee), 


where 
QŸ = 0.30 cos (0.3 uAf) — 0.70 cos (0.1 uAt), 
QŸ = 0.48 sin (0.3 uAt) + 0.44 sin (0.1 At) 
* = 0.45 {cos (0.3 u.At) — cos (0.1 uAf)}, 
QF = 0.32 sin (0.3 uAt) — 0.76 sin (0.1 uAt) 


The values of Q* at different times are given 


ONE) + cos u (x se ) | Qu) o | sin nr 


wb 


which gives 


m 


(30) 


Here Q”) =o for n even and Q®=o for n 
odd, exactly as in the numerical solution (28), 
while 


I 

IB U'= 2 uAf J = | 

in table 1, where uAT= or 17 Q, / (ua = (4 — a | 
: ein n odd er 

For t up to about 2T Q% is still small compared (wAt)? — (n+ 1)? n2| 

with the other ones, which indicates that the ae 

approximation N= 4 is sufficient for describing A { I 

the development up to that time. That this 1s QU = — 2uAt | AD? = (n = 1272 = 

really the case can be seen by considering the i 

exact solution. For the Couette flow equation _ I (| A 4 ER 

(1) reduces to (uAr)? — (n+ 1)?2| 
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The values of Q and Q at different times 
are given in table 1. Up to 2T there is a good 
agreeraent with the numerical solution. 


3. À symmetrical, harmonic velocity profile 
In this section we shall consider a velocity 
given by 


(= st 


n (32) 


Os yeu 
The instability of such a linear flow is possibly 
of some interest for the study of the large-scale 
atmospheric motion. 

For this profile the developments (16) are 
very simple and the equations (17) become 


ea 


(rn: =) CAE | 
Sr | D? n 
\ 


(33) 

el, O 

5 \' D> n+2 

1 | k zu? 

3 Li IN 3) Sudan. >> 

2 
HN wer 
in 24 T(N : Te 
2° LUN) 415 Aan-2+ 


2 
+ (wen) CAN = © | 


This system of equations consists of two in- 
dependent parts, the one containing a, with 
n odd, the other 4, with n even. Accordingly, 
we can consider those perturbations separately 
which are either symmetrical or anti-symmet- 
rical in y 

Taking N=4 the following two equations 
in a, and a, are obtained 
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Cim 


0.3 


0 1 2 3 4 L/D 


Fig. 3. The figure illustrates the stabilizing effect of the 
varying Coriolis parameter showing the imaginary part 
of c given by (35), (B = 0), and of c given by (39) with 


I 70 TT? 
= G GE zi) Aa, = (1 Te oF) Cd; = =) 
They possess non-zero solutions when 
ru Da I A 
4 ED? + 7 ~ 42D? + 72 
/u2D? it 
He DE At on? 


(sutD + 3072u2D2 — 77%) 
(35) 
c becomes complex for wavelength 
12 D=1,234.2°D 


The imaginary part &, of c as function of L 
is shown in fig. 3. It has its maximum value 


about L=2D (=). The two equations in 


a, and a, give 


ae my D? + 12 7° 
2 (u2D? + 4 x?) (u2D? + 16.72) 
which are real for all u. 
In table 2 are given the values of cm) corre- 


sponding to pas for the different approxi- 


mations N=1, 2, — — -, 8. It is seen that 
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Table 2 
ee EN A DET RE BIETET me te BEN RN 
A A A A A A A A 

iE — 0.50 

2) — 0.50 fe) 

3 — 0.25 +0.30i o — 0.25 — 0.30i 

4 —0.25 +0.30i 0.20 — 0.25 —0.30i — 0.20 

5 — 0.38 +0.301 0.20 — 0.38 —0.30i — 0.20 0.27 

6 — 0.38 +0.30i o — 0.38 —0o.30i 0.46 O7 — 0.46 

7 — 0.34 +0.18i fe) — 0.34 -0.18i 0.46 — 0.36 — 0.46 0.54 

8 — 0.34 +0.18i 0.14 — 0.34 -0.18i — 0.14 — 0.36 0.63 0.54 — 0.63 


the flow is unstable only for the symmetrical 
part of the perturbations. 


Considering the initial symmetrical pertur- 
bation 
TX. TH 
= B sin — sin - te 


D D 


® (x, y, 0) 


we get with the approximation N=8 the 
solution 


at wea Fon 
Or — BIS, sin H+ 5, sin 32? +, sin 52? 


7zy\ 


Tos an. 


S, = 1.65 M sin = (x + 0.34 At) + 


+0.55 P cos = (x + 0.34 At) — 2.28 sin = - 


+ (x + 0.36 At) — 0.02 sin > (x- 0.54 At) 


= O LT M sin = (x + 0.34 Al = 


- 0.26 P cos 5 (x + 0.34 At) + 0.21 sin > 


- (x + 0.36 At) + 0.01 sin (x — 0.54 At) 


S5= 0.07 Msin 5 (x+ 0.34 At) - 


+ 0.01 P cos = (x + 0.34 At) - 0.14 sin 5: 


+ (x + 0.36 At) — 0.01 sin — 5 (ie — 0. 54 At) 
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S, = 0.04 Msin (x + 0.34 Af) + 


+ 0.01 Don (x + 0.34 At) - 0.08 sin 


-(x + 0.36 At) + 0.04 sin (x — 0.54 At) 
nAt nAt At 
M 0,18 D = „18 P QE 0,18 D 
(36) 
nAt 


0,18, : 5 
The factor e’ P is doubled in the time 


T=1.3 =. At the times T and 2T the solutions 


el he 


(gga! el nt 
@O*(T)=B 1.98 sin (5 x + 1.64) sinn + 


| 


+ 0,38 sin (58- 0.46) sin “ae 


+ 0.06 sin (5% + 2) sin nr 


+ 
+ 0.02 sin Tx+y® cuil 
D 7 D | 


O* (27) = ls O4 sin (Se 2.97) s sin oe 


8 si pe ln 
+ 0.98 sin | x + 0.79 D 


ad ah Oy on oe 
+ 0.18 sin D*t Par sin 7 


. {2% NZ 
+ 0.10 sın DAT sin 


= 


— 


Q D/2 D 
k (2T 
Fig. 4. The isolines wee 2 ) according to (36), drawn 


for every whole number. 


y® and y® denote certain phase-angles. The 
isolines DX(x, y, 2T)=constant are shown in 
fig. 4. It is seen that the amplitude of the 
perturbation at the time 2T is increased to 
more than sB. 

The function R, which according to (21) 
acts as a source for the error y, has in this case 
the simple form 


R=275AB| -o 6 M cos (x + 0.34 At) + 
=2 53 4 cos (x 34 
SEIT 
+ 0.16 Psin + (x + 0.34 At) + 


+ 0.96 cos (x + 0.36 At) - 


— 0.04 sin > (x — 0.54 Ad) | sin 21 
(37) 
For the integrated effect we have 
D 
TO TC . omy 
Rdt = 0.26 2 = m 
if dt = 0.26 2 53 B sin (5x+ör) sin 7 
2T 
EN TE . omy 
[ Rit=201 ‚2 D? B sin (Hart) SMS, 


0 


where 6 is a certain phase-angle. This induced 
vorticity is for £ up to nearly 2T small com- 
pared with the vorticity v2®* of the approxi- 
mate field. But certainly the error can be 
further increased by the interaction with the 
basic flow U. 

As essential features of the large-scale 
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atmospheric motion can be described by a 
two-dimensional barotropic model, the ques- 
tion arises whether the generation of the large- 
scale disturbances can be explained, fully or 
partly, as an instability of the kind considered 
above. This problem has been investigated 
most thoroughly by Kuo (1949). The numeri- 
cal method used in this paper can give in a 
very simple manner rather detailed results. The 
limitation that the numerical solutions are only 
approximations in a certain interval of time 
does not seem serious as physical factors have 
been excluded, which will make any solution 
of little value after a certain time. 

The zonal flow at middle latitudes is often 
nearly symmetrical with a pronounced maxi- 
mum. As a typical case we can consider a 
linear flow given by (32) superimposed on a 
constant zonal flow. It has been shown above 
that this flow will be unstable for perturbations 
of certain wavelengths, thereby doubling their 


amplitudes in the time T=1.3 7 for the wave- 


length of maximum instability L=2 D. For 
the large-scale atmospheric motion we may 


quite well have conditions making a of the 


order of magnitude 1 day. 


The above results cannot be applied directly 
to the atmospheric motion as the curvature 
and rotation of the earth has not been taken 
into account. Further boundary conditions like 
(3) are not correct in a strict sense for the 
atmosphere. Strictly we should therefore base 
the numerical method on an harmonic analysis 
on a spherical surface. 


A well-known effect from the including of 
the variation with latitude of the Coriolis 
parameter is that of adding a velocity towards 
the west to the velocity of propagation c. 
Also the instability of the linear flow is in- 
fluenced by the varying Coriolis parameter. If 
we include the Coriolis parameter fin equation 
N in the usual way, we get instead of equation 
3 
(U-c)p"-{U-ou+U"-B}p=0o (38) 
Considering pai as a constant in the con- 


sidered interval and using the approximation 

(14) with N=4 we get for a symmetrical 

perturbation the following equation for c 
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PERTURBATION EQUATION FOR LINEAR FLOW 


Ss Ty? D?-: x 300 
wid D? + +702 
toe ets D? - 
(u2D® + 72?) (u2 D? + 9x?) = 
A 
4 (u2D? + 22) (u2D® + on?) 


-\/ (uD? 321%) (u2D° + om?) 
2 
[ (setD4+30 n2u2D’— 7 14 — 32 7° 72) + 
Vo 
+ 256.08 pr. | . (39) 


(39) compared with (35) shows clearly the 
stabilizing effect of the variation of the Coriolis 
parameter. It is seen that the ne of B 


depends upon the magnitude of À ;. The larger 
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is, the smaller is the stabilizing effect of ß. 


pi is 


A 
D? 
The imaginary part of c in the case — 


shown as function of L in fig. 3. 
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A Quantitative Prognostic Method for Thickness Charts 


using 


Advective Tendencies 
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Abstract 


A method is presented for the use of upper winds to construct a chart for the hori- 
zontal distribution of advective tendencies in the layer 1000—$00 mb. A 24-hour prognostic 
thicknesschart can then be constructed from these advective tendencies. The method has been 
tested in a few weather situations in Europe, especially the one of Jan. 30—Feb. 2, 1953. 


In the light of recent developments in quanti- 
tative forecasting techniques it is of some 
interest to explore a little more in detail the 
simple graphical method that may be used in 
the daily forecasting with the present set up of 
an ordinary weather service. We shall here 
discuss a purely advective method by which 
we can obtain a quantitative forecast of the 
thickness pattern. 

Such a prognostic chart for the thickness 
(relative topography), could be used together 
with the prognostic sea-level chart to con- 
struct a forecast for the absolute topography. 
It would also be useful in showing the new 
positions of fronts and possible new develop- 
ments to be introduced on the prognostic sea- 
level chart. 

The necessary analysis of the advective ten- 
dencies for the layer in question further shows, 
in a more distinct manner than the charts for 
fixed pressure levels with contour-lines and 
isotherms, the immediate magnitude of the 
advection, which is useful for any more detailed 
analysis of the mechanism of pressure changes. 


The method by which the thickness changes 


in any particular layer may be computed, is 
derived in analogy with the tendency equation. 
Let A denote thickness (relative topography), 
Ra gas constant for dry air, p pressure, v 
horizontal wind, T* virtual temperature and 
vz vertical velocity. Thus we have 


Po 
= Rf -vvuT+ + 
dt 2 p 
Po N ec oe (1) 
+ Raf ( u men öp 
dt dz} p 
p 


The last integral is called S. 
In most cases we may neglect S, and equa- 
tion (1) is then reduced to 


oA ot represents the change in the amount of 
geopotential unit layers contained in the col- 
umn of air between the pressure levels p, and 


Tellus VI (1954), 2 


PROGNOSTIC METHOD FOR THICKNESS CHARTS 


p, and from equation (2) we see that this change 
depends exclusively on horizontal advection of 
air. By neglecting S we have excluded the 
effects of vertical motion as well as non- 
adiabatic effects. 

Assuming that the thermal wind relation is 
on (1) may be transformed into (cf. Rosspy, 
1942 


where 


The integral has a simple geometrical inter- 
pretation and may be evaluated graphically 


from the hodograph for the layer between pu 


and pp. 

For practical purposes we may assume that 
the triangle between v, and v, represents the 
area swept out by the velocity vector, and 
then equation (3) may be written. 


IA | : 
—=ksin y-v,-v, sine (4) 
It 

where 

oa/dt = the “advective tendency” 

Vi = wind velocity at level 1 

Vo es LE] LE] LE LE] 2 

k = a constant depending on the units. If 


we express v, and v, in knots and 
want to get in gpm/3 hours, the value 
of k = 4.25 -10 2: 


Y = the latitud, counted positive in the 
northern hemisphere 
x = the angle between the wind direc- 


tions in levels 1 and 2 counted posi- 
tive when the wind turns to the right 
with increasing height. 


Some remarks on the construction of a 24 


‘hour prognostic chart for the layer between 


1000 and 500 mb. 


Although the computation accounts for the 
layer 1000 - soo mb, we cannot use the wind 
at the 1000 mb level to represent v,, as it is too 
much influenced by friction. The best alterna- 
tive is to take a wind value as close to the 
upper boundary of the friction layer as possible 
Tellus VI (1954), 2 
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(950 mb), or to interpolate between two wind 
values, one within and one somewhat above 
this layer (e. g. the surface wind and the wind 
at 850 mb). v, is taken at soo mb. 

We assume that the contour lines of the 700 
mb level give a reasonably good picture of 
the velocity field for the layer between 1000 
and s00 mb. Arrows showing the 24-hour 
displacement according to the 700 mb velocity 
field in certain points are placed on the OA /dt- 


Fig. 1. Estimation of 24-hour thickness changes using the 
chart of 3-hour advective tendencies. Full lines represent 
lines for advective 3-hour tendencies. The broken lines 
represent isohypses of 700 mb. The arrows demonstrate 
the estimation of magnitude and position of a 24-hour 
maximum and the location of a point on a zero line. 


chart with their mid-points in the points in 
question. We assume that the 0A/dt pattern 
moves with the 700-mb wind without chang- 
ing shape. The time integration may then be 
replaced by a space integration along the 24- 
hour displacement arrow. Since the advection 
tendencies are expressed in gpm/3 hours and 
we want to forecast the 24-hour change, each 
arrow is divided into eight parts and the ad- 
vective tendencies at each point of the arrow 
(except the terminal point), are added to obtain 
an approximate value of the 24-hour thickness 
change at the terminal point (fig. 1). 

Special care is taken to evaluate the motion 
of maxima and minima as well as zero lines. 
Before the final construction of this chart (prog 
As, 1000 — 500 mb) charts of the past changes 
(Az, 1000-500 mb) should be studied. A 
series of such charts will show not only the 
successive displacement of areas of major 
changes but also variations of intensity. 
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The utilization of past changes has been done 
in analogy with the method described by 
SCHERHAG (1948), while the displacements of 
maximum, minimum, and zero lines are made 
in analogy with the method for construction of 
the sea-level pressure field by the aid of 3-hour 
sea-level pressure tendencies, introduced at the 
Swedish Weather Bureau by Sima. 

The assumptions and approximations on 
which the prognostic method described here 
is based are: 


(a) the last integral in equation (1) is neglected. 
This implies a lapse-rate close to the adiab- 
atic or the absence of strong vertical wind 
components. Besides, we have disregarded 
all non-adiabatic processes. The assumption 
of no vertical motion is in same cases a too 
rough approximation. Since we know, 
however, that there is a close correlation 
between warm-air advection and upward 
motion and between cold-air advection and 
downward motion, the assumption essenti- 
ally influences the magnitude of the 24- 
hour-maximum and minimum; 

(b) the wind is geostrophic; 

(c) the area of the hodograph computed by 
vertical integration is equal to the area 
between the 950 mb wind and the 500 mb 
wind; 

(d) the geostrophic wind field in the 700 mb 
level gives a reasonably good picture of the 
velocity field in the layer between 1000 
and soo mb; 

(e) a linear extrapolation scheme is used. 


Test of the method 


The prognostic method described above has 
been applied to some actual situations. 

The resultant prognostic charts have been 
tested by computing a correlation coefficient 
between forecast and observed thickness 
changes. The tested area, containing 26 grid 
points, was the same as the area by Lönngvist 
(1952). For comparison, a correlation coefti- 
cient was also computed for the ordinary prog- 
nostic thickness charts, constructed by pure 
extrapolation as diseussed by Lönngvist (1952). 

Judging from the results of the checking, 
the proposed method can be regarded as 
giving surprisingly good results in view of 
the approximations we have done. The 
most Important approximation is the assump- 
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tion that the vertical motions can be neg- 
lected. This leads, as a rule, to too large 
absolute values of the computed AA/dt. In 
spite of this we have not found that the prog- 
nostic 24-hour changes are too large. This may 
depend on the fact that the availabe observa- 
tions used for the analysis of the OA /df-chart 
do not as a rule include absolute maximum 
or minimum values. As a consequense this 
means that the analysed JA /dt-chart shows a 
smoothed horizontal field of “advective tend- 
encies”. Furthermore past changes have been 
used as a guide in the final construction of the 
prognostic thickness chart. 


Table I. Correlation between predicted and ob- 
served thickness changes obtained by different 


methods 
| Correlation | 
coefficient 
Dante | The 
Ordinary] « ver 
method | AO 
| method”’ 
22 23, OCts TOO 0.45 0.68 
14. D NOV LOSOm eee 0.18 0.63 
1-2 NOV LISTE 0.52 0.80 
30—31 Jan. 1953 ...... 0.68 0.92 
31 Jan.—ı Feb. 1953 0.79 0.92 
2 eo, ROSE Var 0.81 0.96 


I ISCH SEES I 
a Sn J 7 


8% 


3 


\ 0 
er 
> N 
N eG 
5 ad iis 
RT Lu 
< 7 
7 D > 
GE 
A x KL 
AE) 
FA À ASE 
% 
g 4 Ze 
3 7 
7 SI 
7 SN 
0 ON 


Fig. 2. Chart showing advective 3-hour tendencies of Jan. 

30, 1953, 03 GMT. The lines are labelled in tens of gpm/3 

hour. The arrows show the predicted direction of the 

displacement of each area of advection. The terminal 

point of each arrow shows the location of a 24-hour 

maximum or minimum, obtained by the movement of 
the area in question. 
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Fig. 3. 24-hour thickness changes 1000—s00 mb, Jan’ 


30—Jan. 31, 1953, 03 GMT. 
a) Predicted changes. b) Observed changes. 
The lines are labelled in tens of gpm. 


The author has tried to be as objective as 
possible when constructing the prognostic 
charts. However, in the three cases from 1953, 
his knowledge of the subsequent development 
may have influenced the results somewhat. 


Application to the situation Jan. 30—Feb. 2, 
1953 


The weather situation over the Atlantic and 
Europe during the last week in January was 
characterized by a trough extending from the 
Lofoten Islands to southern Greenland and a 
high over France and the Bay of Biscay. In 


Tellus VI (1954), 2 


Fig. 4. 24-hour thickness changes 1000—500 mb, Jan. 
31—Feb. I, 1953, 03 GMT. 
a) Predicted changes. b) Observed changes. 


the south-westerly current south of Iceland and 
west of Scotland waves developed, which at 
first moved quite slowly north-east towards 
Scotland where they in most cases accelerated 
considerably and their amplitude increased. Th: 
resulting disturbances subsequently followed 
the WNW stream across southern Scandinavia 
to Poland where they gradually weakened. 


Jan. 30, 1953 

At 03 GMT a wave between Iceland and 
Scotland was moving ENE, i. e. parallel to the 
flow at about 700 mb. However, at the same 
time an intense cyclogenesis had occurred at 
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Newfoundland. In connection herewith astrong 
advection of warm air took place towards 
Greenland, and the cold air, which had long 
been stationary north of Iceland, began to flow 
south in the rear of the wave between Iceland 
and Scotland. The very strong advections of 
cold and warm air are clearly shown ın fig. 2. 
This chart was taken as the base of a forecast for 
the thickness 1000 - soo mb, until Jan. 31 03 
GMT. The resulting forecast map shows quite 
clearly that the advection of cold air initiated 
south of Iceland is a factor of outstanding 
importance, and indicates a cyclogenesis (fig. 
3 a) over the British Isles in the soo mb level. 
Fig. 3 a indicates the following development 
in the upper troposphere: a rapid displacement 
castwards (to the British Isles), and also an 
intensification, of the trough now extending 
SSW from Iceland. (cf. fig. 3 aand 3 b). 


Jan. 31, 1953 

The maps for sea-level and 700 ınb at 03 
GMT show an intense cyclone just north of 
Scotland; exceptionally high surface wind ve- 
locities were prevailing mainly over Scotland 
and Iceland. The strongest advection of cold 
air is now to be found over the British Isles, 
and the intensity of the advection of warm air 
towards southern Scandinavia has decreased. 
The strong advection of warm air which the 
day before had its centre south of Greenland 
has moved east to the sea between Greenland 
and Iceland, and has maintained its intensity. 
Displacing the area of strongest cold-air advec- 
tion now over the British Isles, we find that 
the maximum 24-hour thickness fall will occur 
over southern Germany. The strong-warm air 
advection south west of Iceland will move 
north-eastwards. The resulting prog A 5,1000 — 
soo mb (fig. 4 a) again shows a high correla- 
tion with A, 1000 - 500 mb (fig. 4 b). 


Feb. 1, 1953 


The deep low has now its centre near Ham- 
burg, and a full north-westerly gale prevails 
over the North Sea, the eastern parts of the 
British Isles, the Netherlands, and Belgium. 
The advection of warm air originally south of 
Greenland has now reached Iceland with 
slightly decreased intensity. As this flow of 
warm air continues in a north-easterly direc- 
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Fig. 5. 24-hour thickness changes 1000—500 mb, Feb. 
1—Feb. 2, 1953, 03 GMT. 
a) Predicted changes. b) Observed changes. 


tion, it cuts off any further advection of cold air 
to the European cyclone, and at the same time 
as the low becomes thermally symmetric, it 
fills slowly. The high correlation between prog 
A 1000 — soo mb and Az, 1000 — 500 mb is 


illustrated by fig. 5 a and fig. 5 b. 
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Letters to the Editor 


On Syönos Initial stages of Typhogenesis 


Dear Sır, 


In a recent paper in Tellus (Syöno, 1953) a 
physical picture of a typhoon model is presented. 
First, it is questioned whether the stages described 
by Syöno should be referred to as “typhogenesis’”. 
He states that: “It is well-known that at the initial 
typhogenetical stage, there occur numerous scattered 
local convective showers and squalls, but no system- 
atic wind. Then we may assume that there is no 
mean systematic converging current within the 
region, where the typhoon is about to form.” Such 
a picture of the formation of tropical storms is in 
sharp disagreement with the views of many authors 
(c.g. RIEHL, 1951; PALMER, 1952) who visualize the 
formation process as the intensification of an existing 
disturbance which often has been evident as a 
synoptic entity for days. The fact that many of the 
travelling disturbances in the easterlies and along 
the equatorial trough—some of which may deepen 
to typhoon intensity — produce rainfall amounts of 
the same magnitude as found in typhoons suggests 
that very definite convergence-divergence patterns 
exist prior to the typhoon stage. The divergence of 
opinions cited may partly be a matter of terminology, 
but it would seem that the Syöno model for the 
initial stage would be satisfied by innumerable small 
scale systems which never develop beyond the 
convective cloud cell stage. It would therefore be 
misleading to refer to these minor developments as 
“typhogenesis”’. According to the other viewpoint 
“typhogenesis” should properly refer to the mech- 
anisms whereby an existing tropical disturbance, 
having cold core characteristics in most cases, is 
able to take on the thermal and circulation charac- 
teristics of the typhoon. 

Secondly, it would seem that the mechanism 
proposed by Syöno for the pressure fall in the initial 
stage of “typhogenesis”’ is open to some doubt. He 
assumed that non-divergent flow is found about the 
initial convective cell so that the local pressure 
tendency is dependent only on the vertical motion 


: de 
term of the tendency equation; i. e. = gow 


The vertical mass transport term was separated into 


that due to moist air (ow), and to liquid water (ow,)1. 


d — + 
En — g( ow), + g(ow), 
C 


After making the statement that: “At the present 
(initial) stage, since numerous local convective cells 
are surrounded by downward currents, the mean 
upward mass transport in the form of air is expected 
to be very small”, the local pressure change is 
discussed by considering only the liquid water 
transport. From this reasoning the conclusion is 
reached that: “In fact the pressure ‘tendency is 
ı mm/hour corresponding to an intensity of the 
rainfall of 10 mm/hour.” Static reasoning shows 
that this statement can be true only if the precipi- 
table water in a vertical column has been decreased 
by this amount — which is a sizeable percentage of 
the total precipitable water even in a tropical 
atmosphere-since it has been assumed that no 
divergence is present. Clearly, from this argument 
the most effective means of obtaining a pressure fall 
would be rapid dissipation of the cloud system and 
a complete precipitation. Since the intensification 
of a convective cloud in the usual sense is associated 
with an increase in precipitable water in vertical 
columns through the cloud, Syöno’s mechanism for 
the initial pressure fall does not appear physically 
realistic. It could, of course, be argued that a net 
mass convergence of liquid water could take place 
in a non-divergent system provided that it is exactly 
balanced by mass divergence of air. However, the 
mechanics of such a system would appear to be as 
obscure as that of the pressure fall which it purports 
to explain. 


! The division of the transport term in this manner 
presumably implies that the liquid water is in hydrostatic 
equilibrium. If the water is falling such may not be the case. 
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LETTERS TO THE EDITOR 


Reply 


On my work Mr Jordan raised some questions, 
which are caused partly by the difference of defi- 
nitions. For the formation of typhoon, we should 
consider two mechanisms, one is concerned with 
the initiation of the process, in other words the 
mechanism of trigger action, and another is the 
mechanism of the continuous development to 
revolving storms. In mathematical language, those 
corresponding to initial conditions and differential 
equations respectively. In my article, I solved the 
differential equations under certain initial conditions. 
Mr Jordan’s questions are concerned with the initial 
conditions. The term “typhogenesis” is new, so far 
as I know. I used it just for the purpose to distinguish 
the process of formation of a typhoon from that ofan 
extratropical cyclones. So far as a definite meaning, 
which is universally accepted, has not been achieved, 
it may be possible that different authors use a term 
in more or less different meaning according to their 
opinions. I don’t think the first stage of my picture, 
strictly speaking not my own but one of many authors, 
is in sharp disagreement with what he describes. 
The concepts of “no systematic wind” in my article 
and of his “‘very definite convergence-divergence 
pattern” must be different, so far as the scale of 
pattern, which he considers, is not the same with 
that of typhoon in her incipient stage. I would like 
to ask him, if it is possible to get a definite pattern 
with the scale of typhoon by using sparce data 
obtained in the original region of typhoon in the 
ocean. Mr Jordan rewrites my sentence “‘numerous 
scattered local convective showers and squalls” as 
“innumerable small scale systems” and then he says 
conclusively “which never develop beyond the 
convective cloud cell stage”. Therefore I think that 
it is not necessary to answer this point. But here I 
would like to add that some mechanism should 
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exist as the cause of initial rainfalls. This mechanism 
must be studied by synoptic aerological analysis by 
using sufficiently abundant data. On this mechanism, 
Dr Riehl’s hypothesis (Rıekt, 1951) might be 
probable, but on this hypothesis I have the same 
opinion as Dr Rossby (RossBy, 1948). 

Secondly, Mr Jordan raised a question to the 
tendency equation. But he makes some mistakes. 
My estimation of pressure tendency by rainfall is 
not a conclusion but a remark, because the effect 
is currently neglected. Next, his footnote is unac- 
ceptable, because raindrops fall with final velocities 
by the frictional effect of air. Next, I explained the 
process by dividing it into three stages, in order 
to get clear physical pictures. But actually, these 
stages cannot be cut distinctly. The transition from 
one stage to the next is continuous and two stages 
are superposed. The amount of precipitable water 
in a tropical airmass is more than 30 mm (equiva- 
lent pressure is 3 mb), then far before exausting all 
precitable water, the second stage begins. Therefore, 
Mr Jordans worry about clearing up of clouds is 
unnecessary. But I believe that a discussion like 
this should be done quantitatively, or it will be an 
endless dispute. I shall not hesitate to revise my 
physical picture, if some one will show acceptable 
proof based on proper data and reasonable dis- 
cussions, but not objections in words, that my 
picture contains something wrong. 
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